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PREFACE
	The Electric Punched Card Method, known also as the
Hollerith, or Electric Accounting Machine Method, is now
generally known to the business and statistical worlds. Its
use in physical science has, however, been confined largely to
astronomy.
	The purpose of this book is to show the possibilities of the
method in scientific computation. It is not a manual of
operation for the various machines but is designed to enable
a scientist so to formulate his problem that any skilled operator
of the machines could carry it out.
	No account is given of punched card tabulating machines
in which the cards operate the machines mechanically instead
of electrically. The great advantage of the electric
method for scientific computation is its flexibility. In the
electric method the necessary machine changes in going from
one problem to another can be accomplished by the operator
in a few minutes.
	Nearly all astronomical calculations now being done by the
Punched Card Method are being performed in The Thomas J.
Watson Astronomical Computing Bureau, a scientific non-
profit-making institution which was organized for the use of
astronomers. Most of the methods in the book will therefore
be described in terms of the equipment of the Bureau. This
procedure will assist astronomers in using the Bureau, and will
furnish for others a working example. A brief description of
the organization of the Bureau is given in appendix C.
	The computing laboratory which is now operated by the
Bureau was developed by the writer, first in the Columbia
University Statistical Bureau and later in the Department of
Astronomy of Columbia University. The directors of these
institutions, Dr. Ben Wood and Dr. Jan Schilt, have given
every encouragement.
	The laboratory would not have been possible without the
interest and generosity of Dr. Thomas J. Watson, President
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of the International Business MachineS Corporation. The
undertaking has had his cooperation and encouragement from
the very beginning.
	It is a pleasure to acknowledge the technical and other
assistance of many officials of the International Business Machines
Corporation. Mr. G. W. Baehne made every eITort to
secure equipment that would effectively meet all of the requirements.
Mr. William Fitzpatrick, Customers' Service
Supervisor of the New York District, assisted by Mr. W. L.
Seasholtz of the Customers' Service Division, arranged the
special circuits for the process of numerical integrat ion described
in chapter VIII. The Bureau itself has grown out of
a suggestion by Mr. C. H. Tomkinson, then Manager of the
Commercial Research Department.
	Miss Lillian Feinstein has contributed largely to the success
of several of the applications in part III. In connection
with the work on the Yale Zone Catalogues and the Columbia
Photometric Program she has made valuable suggestions as to
the arrangement of the work and has personally performed
most of the computations.
	The author has been able to incorporate in the manuscript
a number of corrections and improvements suggested by Dr.
Frank Schlesinger, Dr. Dirk Brouwer, Mr. Hans Hertz, Mr.
George Essenwine, Miss Lillian Feinstein, Mrs. Grace O.
Savage and Dr. Robert Wolf.
				W. J. Eckert
Rutherfurd Observatory
Columbia University
New York City
December 15, 1939.
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PART I. THE PUNCHED CARD METHOD

CHAPTER I

INTRODUCTION

	There have been three major advances in computing
methods brought about respectively by the introduction of
logarithms, of manually operated computing machines, and
of automatic tabulators using punched cards. The first table
of logarithms was publishced by John Xapier in 16114, and for
nearly three centuries all extensive calculations were performed
by this method. The present generation has seen the
logarithm replaced by hand and electric computing machines.
These arc all similar in that the simple operations of arithmetic
are performed more or less automatically upon numbers
which are put into the machine by hand. For the automatic
tabulators the initial numbers arc punched on cards and
checked. From then on they receive no personal attention,
as they are automatically transferred through flexible wiring
to the proper place in the machine for the various operations.
	The complete line of equipment using the punched card
system is too extensive to describe here; suffice it to say that
there are standard machines for practically all purposes.
They have been developed for commercial uses but with
slight modifications become ideal for scientific work. The
main question in any case is not "can the problem be solved
by these machines?" but rather "have l enough operations
of this type or that to justify such powerful equipment?"
It might be mentioned also that, with standard parts, special
machines for almost any imaginable purpose could be constructed,
but except in very rare cases the expense and delay
involved make the standard machines preferable.
	The first extensive use of the early Hollerith Tabulator in
astronomy was made by Comrie. He used it for building a
[Footer: 1]
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table from succcssive difkrcnces, and for adding large numbers
of harmonic terms.*
	The laboratory at Columbia Univcrsity is the first capable
of performing intricate scientific calculations by machinery
without anv reading or writing of figures. The individual
machines have highly specialized functions, but since each is
capable of reading numbers from the punched cards and recording
the results in the same manner it is possible to complete
a long sequence of operations one after the other without
writing down a single number by hand.
	A few of the specific astronomical problems which have
been solved in this laboratory will be described in detail in
later chapters, but to give some idea of their nature brief
mention will be made of them here in general terms.
	One of these problems is that of planetary motion. It
consists of the numerical integration of a set of three simultaneous
differential equations each of the second order. The
solution is carried out with ten significant figures and at the
rate of twelve epochs per hour.
	Tables of functions are constructed from their differences
with great efficiency, either as printed tables or as a file of
punched cards. From the latter type the values of the
function may be taken automatically at the rate of many
thousands per day. Calculations involving such tabulated
functions are efficiently performed.
	The manipulation of harmonic series is useful in many
kinds of theoretical work. Such operations consist of numerical
harmonic analysis and synthesis and the multiplication
of harmonic series.
	The formation and the solution of normal equations, the
computation of correlation coefficients, of averages, etc., are
natural applications of the punched card method.

	*On the Construction of Tables by Interpolation. L. J. Comrie. Monthly
Notices, R.A.S., Apr. 1928. The application of the Hollerith Tabulating Machine
to Brown's Table of the Moon. L. J. Comrie. Monthly Notices, R.A.S., 
May 1932.
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	In one program the punched card has been found to be t he
best means of recording astronomical observations; observations
recorded in this manner may be automatically reduced.
	The results of punched card calculations are easily listed
for reference or for publication. The tabulator record may
be used not only for printer's copy but may in many cases be
reproduced directly by photography.

[Page Break]

CHAPTER II

THE PUNCHED CARD AND THE MACHINES

1. THE OPERARING PRINCIPLE OF THE HOLLERITH METHOD

	In the Hollerith method of computation the numbers are
recorded by means of holes punched in cards, and are read
into the machines by means of electrical contacts made
through the holes. In all the machines the principle is the
same. The tabulating card, acting as an insulator, passes
between a wire brush an<l a brass roller.

[
DIAGRAM of READING BRUSH

Text includes following:
"BRUSH"
"BRUSH 
HOLDER"
"TABULATING
CARO"
"CONTACT
ROLLER"
"FIG. 1. Reading Brush of Punched Card Machine"]

A hole punched in the card permits the brush and the roller
to make contact, thus completing an electric circuit which
contains an electromagnet. The position of the hole in the
card determines the time when the magnet is operated, and
all parts of the machine which are synchronized with the
passage of the card will therefore perform different functions
depending upon the position of the hole in the card.

[FOOTER: 		4		]

[Full Page Image of 80 column punched card, in addition to the numbers on the card the following caption is written vertically on the right side of the image: 
"FIG. 2 Tabulating Card (Actual Size)"   ]
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	In the sorting machine this magnet opens a chute along
which the card slides until 1t reaches the proper receptacle.
In the adding machines, the digits on the card are added.
In the multiplying machine the multiplicand and multiplier
are read from the card.

		2. THE TABULATING CARD
	The tabulating card (fig. 2) has 80 vertical columns with
12 punching positions in each column. Of these, ten are
indicated by the printed digits 0 to 9 which correspond to the
digits of the numerical data to be punched. The 11th and
12th punching positions at the top of the card (not indicated
in fig. 2) are commonly referred to as the "x" and "y"
positions. The "x" punch is usually used to control the
operation of the various machines. For example, negative
numbers are differentiated from positive ones by means of
the "x".*

[Small image of a partially punched tabulating card captions reads:
"FIG. 3. Tabulating Card Partially Punched" ]

	Figure 3 shows a card with data punched in the first six
columns. The punching reads : 1 3 8 y 6 x 9 0. The vertical 
lines have nothing to do with the punching, and are described
in the next paragraph.


* For the alphabetic machines each letter is indicated by two punched hole,
in one column the "x" or "y" being combined with one of the digits 0 to 9.
See sec.4.
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	The first step in the use of a card for a particular record
is the designation of groups of columns as "fields." Each
field defines a section of the card in which one particular type
of information will always appear. If a given problem is
large enough to warrant the expense, the cards may be ruled
at the factory to indicate and label the fields.* Most of the
computation of the Thomas J. Watson Astronomical Computing
Bureau is done on a card with standard ruling designed
only to assist the eye in reading from desired columns.
Figure 3 shows the standard tabulating card used by the
Bureau. The 80 columns are divided into fields of ten, with
further subdivisions of three and one.

[Image of a punched card, punched wiht the punches interpreted at the top and columns separated visually by ruled lines with groups of columns given abreviated text headings, captions read:
"FIG. 4. Tabulating Card with Interpretation and Special Ruling" ]

	Figure 4 shows a card which was ruled for an extensive
project. There are over 30,000 cards in the set and ten
duplicate sets are to be used at different places so the extra
expense of ruling was justified. The numbers at the top are
not part of the ruling but were added later by the "interpreter" 
(see sec. 3).

3. THE MACHINES OF THE THOMAS J. WATSON
ASTRONOMICAL COMPUTING BUREAU

	The machines of the Bureau are, with one exception,
standard accounting machines. Where necessary, circuit

	* There is no charge for the ruling but the preparation of a new electroplate
for a change in ruling costs about $25.
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changes were introduced, but the original functions of the
machines have not been disturbed. The exception is the
"Calculation Control Switch" which controls the operation
of several other machines in the complicatcd process of numer-
ical integration. In describing these machines we shall there-
fore give the trade name of each machine and note the
changes that have been made.

The Card Punch
	The card punch is used to transfer the given data of a
problem from the manuscript to the punched cards. The one
used for most of the work of the Bureau is a "motor-driven
electric duplicating key punch." In this model the cards are

[Image of a duplicating key punch, text is as follows:
"Duplicating
Rack"
"Keyboard" 
"New
Cards"
"FIG. 5.  Duplicating Key Punch" ]

fed and ejected by an electric motor, the punching is done
magnetically, and information from one card can be automat1cally
transferred to another by the duplicator.

[HEADER: TnE PUNCHED CARD AND THE MACHINES		9]

	The operation of the punch is much like that of an ordinary
typewriter. The card is fed into the machine with the
first column under the punclws. The depression of one of
the twelve punching keys causes the corresponding hole to
be punched in that column, and the card to move to the next.
Columns may be skipped by pressing the spacer key or the
card entirely ejected by pressing the eject key. Provision is
also made for automatic spacing and ejecting on any predetermined
column.
	If a prepunchcd "master card" is placed in the duplicating
rack, the information from the master card will be copied
onto the card in the punch. Thus, for example, we might wish
to punch common information in columns 1-5, 13, 14 of a
large number of cards, columns 6-12 being different for each
card. The common information is punched in the "master
card" and it is placed in the duplicating rack. As soon as a
card is fed into the punch, the duplicator copies columns 1-5,
the operator punches columns 6-12, and the duplicator copies
columns 13, 14. When part of a card is spoiled the duplicator
is used to copy the good portion.
	The speed of the punching operation depends not only
upon the operator but upon the legibility of the manuscript,
and the amount of coding required. The duplicator operates
at the rate of ten columns per second, and about one and a half
seconds are required to feed a card. A general idea of punching
speed may be obtained from a recent job where about
half of each card was filled, part by the duplicator, and with
a small part of the information in code. On 30,000 cards·an
average of 200 per hour was maintained. On another extensive
job, where about half the card was punched but considerable
coding and rounding of numbers was required, the
operators averaged less than 150. For straight numerical
punching from good typed manuscript 125 eighty-column
cards should be punched per hour.

The Verifier

	The mechanical verifier is used to check the punching of
the initial data on the cards. It is operated in the same
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manner as the key punch. When a key is depressed the
card is spaced provided the correct hole has been punched;
if not, the card stops. For most purposes this method is 
more accurate than proofreading. It does not, however,
check double punching, i.e., two holes in the same column.

The Counting Sorter
	The sorter automatically sorts the cards into groups ac-
cording to the information punched in any chosen column
of the cards. The cards to be sorterd are placed in the hopper,

[Image of a counting sorter, text includes "Counting Attachment" "Card Hopper" and
"FIG 6. Counting Sorter" ]

the movable sorting brush is set on the desired column, and
the sorted cards fall into the various receptacles. There
are thirteen receptacles, twelve corresponding to the twelve
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punching positions on the card and the thirteenth called the
"reject" for cards which had not been punched in the column
of the sorting.
	Provision is made for closing off any of these receptacles.
For instance, if we wish to pick out only the cards with an "x"
punch in the required column, we can shut off all the others
so the cards without an "x" will go into the reject.

[Image of sorting brush in operation arrow points to the brush caption reads:
"FIG. 7. The Sorting Brush in Operation" ]

The counting attachment gives the number of cards which
goes into each compartment, the total for the group, and a
grand total of all the groups. Switches are provided so that

[HEADER: 12		PUNCHED CARD METHODS	]
counting may be done without sorting, sorting without counting,
or sorting and counting.
	Figure 7 shows the movable sorting brush which may be
quickly set on any desired column. It also shows the card
in the position where the actual sorting takes place. At this
point the card passes into a chute leading to the proper
receptacle. The arrow shows the brush making contact
through the hole in the card.
	Sorting is done at the rate of 20,000 cards per hour.

The Tabulator

	The tabulator is essentially a high-speed adding machine
of large capacity which reads the numbers to be added from
holes in punched cards. When equipped with a printing
mechanism it is called an accounting machine. This printing
mechanism will print information from each card or the accumulated
totals only.
	The feeding mechanism of the tabulator has eighty brushes
of the type shown in figure 1, one for each column of the card.
As the card passes these "adding" brushes the adding wheels
in the "counters" revolve and the electromagnets controlling
the addition are actuated by the electric contacts made
through the holes in the cards.
	The electric circuits from the adding brushes to the
electromagnets of the adding mechanism pass through a plugboard
similar to that of a telephone exchange so that it is
possible, by means of plugwires, to select the positions on the
cards from which the numbers are to be taken and the positions
in the adding mechanism where they are to be added.
	The machine used by the Bureau is the so-called "type
3S Direct Subtraction Accounting Machine." It has been
modified to some extent. The standard machine of this type
has five adding units or "counters" each with eight-digit
capacity, whereas this one has four ten-digit ones. A few
special circuits have been added for the application described
in chapter XII, but these do not interfere with the normal
operation of the machine. It has the following features which
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are optional on all standard machines: "automatic plugboard,"
" major and minor control," and two "class selectors."
	The automatic plugboard is an arrangement whereby the
plugboard may be removed from the machine for convenience
in wiring. Several extra plugboards are provided so that a
number of wiring arrangements may be quickly available.

[Image of Subtraction Tabulator text includes "Card Feeding Mechanism"
"Printing Mechanism" and the caption at the bottom "FIG. 8. Subtraction Tabulator"]

	Automatic "control" is a feature of most tabulators. For
control purposes the feeding mechanism contains, in addition
to the "adding brushes" mentioned above, a second set of
eighty brushes known as the "control brushes." While one
card passes the adding brushes the following card passes the
control brushes. During the feeding process it is therefore
possible to compare the desired part of each card column by
column with that of the preceding. The principal use of the
control is to enable the machine to sense when one group of
cards has been finished and the next is about to begin. For
example we might have cards numbered consecutively from
one to one thousand and wish to tabulate them in groups ot ten.

[Full page image of an automatic plugboard with text: "Fig. 9. Automatic Plugboard"
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We would wire the controlling mechanism for the columns
on the card containing the tens and hundreds digits. In the
two columnc; selected for the control the first nine cards would
match and would "pass control." The tenth card however
would differ in the tens position and would therefore "break
control" and start a new "control group." This break in control
can be used to cause the tabulator to perform various
functions such as printing totals from the counters and clearing
them for the next group. By the use of major and minor
control it is possible in the above example to tabulate by tens
and also by hundreds. Any number of columns up to sixteen
may be used for control.

[  Diagram illustrating how feeding mechanism of a tabulator works. Text includes "CARDS TO BE TABULATED", "CONTROL BRUSHES", "ADDING BRUSHES" and "TABULATED CARDS".
The caption reads: "F1G. 10. Feeding Mechanism of Tabulator" ]

	Figure 10 shows the card feeding mechanism of the tabulator.
The cards to be tabulated are placed in the hopper at
the top. They arc fed one at a time from the front of the
pack and moved downward past the control and adding
brushes, and are finally stacked in a second hopper.
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	The class selector makes it possible to add in a given
counter from one of two different fields on, the card, depending
upon whether or not the card has an "x" in a certain column.
The use of two class selectors in series makes it possible to
select from three fields.
	The individual counters of the tabulator can be made to
add or subtract the numbers read from a card depending upon
whether or not the card contains an "x" in a given column.
The "x" controls are wired through the plugboard and hence
give great flexibility. For instance we can add the "x" cards 
and subtract the "no x" ones, or add the "x" cards and omit
the "no x." By the use of "x's" in two columns it is possble
to subtract algebraically the numbers on one set of cards from
those on another (see app. B). In addition to reading the
numbers from the cards a special circuit is provided for reading
a "l" with the passage of each card. When this circuit is
wired to a counter it gives a count of the cards tabulated.
	The printing mechanism consists of five printing "banks":
one connected to each of the four counters and one, known as a
"list bank," which prints directly from the cards without
adding. It is thus possible to print fifty digits from each
card and to add forty of them in the counters simultaneously.
When each card is listed separately the process is known as
"listing" but when only major or minor totals are printed
the term "tabulating" is employed.
	In connection with the use of control groups there is a
process known as "group indication." This means that for
each control group the numbers selected for group indication
will be printed or added only once. Information from
fields other than those used for the control may also be wired
through the "group indicator."
	The tabulator operates at the rate of 9,000 cards per hour
while tabulating, and 4,500 while listing. The time required
to break control and reset is 1.1 seconds.

The Summary Punch
	The summary punch is a motor-driven duplicating key
punch llttached by an electric cable to the tabulator. Its
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purpose is to punch on a new card the totals accumulated in
the counters of the tabulator. It may be operated by a
major or a minor break in control. A plugboard on the
summary punch makes it possible to select the positions from
the counters to be punched in any column of the card. By
throwing a switch on the summary punch it is possible to
change the signs of all the numbers punched from a given
counter.
	Information common to many cards may be punched from
a master card in the duplicating rack as with the key punch,
so that part of the information on a given "summary card"
may come through the tabulator and the rest from the master
card.
	The summary punch operates at the rate of ten columns
per second, and 1.5 seconds are required to feed a new card.

The Multiplier

	The multiplying punch reads two numbers from a card,
forms their product and punches it on the same card. The
"cross footing"* feature makes it possible to add or subtract
as well as to multiply two numbers from a card. The two
machines used by the Bureau have the crossfooting feature,
and one of them has an automatic plugboard.
	The following operations are possible:
	A ⅹ B		A ⅹ B + C + D
	A ⅹ B + C	A ⅹ B - C - D
	A + B - C

	The possibilities of the multiplier are best understood
from a description of its construction and operation.
	The feeding mechanism is similar to that of the tabulator.
There are 80 reading brushes which correspond to the adding
brushes of the tabulator and which are wired to the plugboard.
The eighty control brushes of the tabulator are here
replaced by three or four movable brushes which may be set

* Crossfoot is an accounting term which means to add numbers across a line
rather than down a column.
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in the desired columns Their function is limited to the de-
tection of "x" punches and they are commonly referred to as
"x" brushes. When a card has passed the brushes it goes
into a card punch where the answer is recorded on it.
	The multiplying machine contains seven counters of the
type used in the tabulator: six of eight-digit capacity and
one of ten. They are referred to as (1) multiplier; (2) multi-
plcand; (3)-(4) right-hand components; (5)-(6) left-hand
components; (7) summary products. Counters (3)-(4) are
connectcd in series to form one sixteen-digit counter, and the
same is true of counters (5)-(6). Counter (7) is the ten-digit
one.
	In straight multiplication the operation is as follows: As
the card is fed past the adding brushes, the multiplier and
multiplicand are read from it into counters (1) and (2) respectively.
The multiplication then starts. The first d1g1t of
counter (1) multiplies each digit of counter (2) electrically
and at the same time. Thus if the first digit of the multiplier
is a 9 and the number in the multiplicand is 7638, the
partial products 9 ⅹ 7 = 63, 9 ⅹ 6 = 54, 9 ⅹ 3 = 27, and
9 ⅹ 8 = 72 are formed from the "electrical multiplication
table" all at one time. These partial products together with
those which will come from the remaining digits of the multiplier
must be accumulated in the 16-digit counters of the right
and left-hand components. The reason for two sets of counters
is that only one digit at a time can be added on a given
adding wheel. The right-hand components of the above
operation are of course 3, 4, 7, 2 and the left-hand components
are 6, 5, 2, 7.
	After the multiplication by all the digits of the multiplier
is completed, the contents of the right-hand components
counters are transferred to the left-hand components counters.
The results are then punched on the card and also transferred
to the summary products counter. This counter accumulates
the sum of all the products of the group.
	In punching a given product on the card, the entire
product is usually not required. For instance in forming the
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product of two eight-digit numbers we usually want the product
to eight figures. In dropping the extra digits the last
figure is usually "rounded" to the nearest unit. This rounding
is automatically done by the machine with the use of the
"half pick up" whereby a "five" is automatically added in
any desired column of counter (4) immediately after the card
passes the reading brushes.
	The crossfooting operations are accomplished by reading
into counters (4) and (6) at the same time as (1) and (2).
When the numbers arc being transferred from (5)-(6) to
(3)-(4) they may be transferred as true figures or as complements,
thus giving the subtraction crossfooting. On some
of the later models it is possible to read also into counter (7).
	For purposes of checking, a switch is provided which
interchanges the multiplier and multiplicand, and spaces the
card without punching. The sum of the products is then
found for comparison with the value obtained during the
previous run.
	In ordinary multiplication or crossfooting all the counters
except the summary products counter are cleared after each
card. When a common multiplier is to be used for several
cards it is possible to retain the multiplier in the counter.
This may be done in two ways: "fixed multiplier," or "rate
card." In the case of the fixed multiplier the multiplier
counter must be reset by the operator. In the "rate card"
method each card containing a new multiplier is identified by
an "x" punch. This "x" punch causes the old number to
be cleared out and a new one to be put into its place. The
multiplicand from each card is thus operated upon by the
multiplier from the last preceding rate card.
	Special features of our machine, aside from those mentioned
in chapter XII, are:
	(1) device for reading the algebraic signs of multiplier
and multiplicand from "x" punches on the card, and indicating
the sign of the product in the same manner.
	(2) A special circuit for forming the quantity A ± (B ⅹ C)
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with the sign depending upon whether or not the card carries
an "x" punch in a given column.
	(3) Provision for additional crossfooting. The standard
machine docs not provide for the use of all sixteen possible
positions in counters (3)-(4) and in (5)-(6).

	The machine operates at 150 cycles per minute. Two
cycles arc required to feed the card, one for each digit of the
multiplier (zeros excepted), and one for the transfer from
counters (3)-(4) to (5)-(6). The speed is independent of the
number of digits in the multiplicand or in the crossfooting.
Punching requires one second for ten digits, but part of the
punching is done on the feed cycle of the next card. Thus
eight by eight digit multiplication is done at the rate of
730 cards per hour when ten digits of the answer are punched,
and at the rate of 870 for "non-punch." Three by eight
multiplication, "non-punch," is performed at the rate of 1,350.

The High-Speed Reproducer
	This machine is used to transfer information from one
card to another. It is much faster and more flexible than
the duplicator of the key punch. An automatic plugboard
allows conyenient changes of wiring. Comparing relays check
the new card against the original during the same run of the
cards.

[ Schematic representation of function of high speed reproducer
text includes: "A" "READING BRUSH" "COMPARING BRUSH" "B" "PUNCHES" "PUNCH BRUSHES".
Caption reads:
"FIG. 11. Schematic Diagram of High-Speedd Reproducer" ]

	There are two complete card feeding assemblies similar
to that of the tabulator. In one of these feeds the control
brushes have been replaced by eighty punching knives. Fig-
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ure 11 shows the paths of the cards through the two feeding
mechanisms. The two may be operated independently, or in
synchronism with a card in feed "B" in exactly the same
phase as the corresponding card in feed "A." In addition
to the brushes shown in the above diagram there are several
movable "x" brushes of the type used in the multiplier.
	A file of cards may be reproduced by passing it through
feed "A" II while new cards are passing through feed "B."
The reading brushes of "A" are wired to the magnets operattng
the punches of "B," and the comparing brushes of "A"
and "B" are wired to the "comparing relays." The purpose
of the comparing relays is to compare the holes in the two
cards which arc passing the brushes connected with the relay.
If the holes do not match, the machine stops.
	"Gang punching" gives a large number of duplicates
from one master card. Gang punching is done in feed "B."
The master card is put in front of the new cards and the
information from each card is passed back to the one following.
For this operation the comparing brushes of feed "B II are
wired to the punches.
	While feed "B" is being used for gang punching, feed
"A" may be used independently for checking.
	Where several master cards are interspersed among the
"detail" cards, an "x" punch in the master card will prevent
its being punched.
	The reproducer operates at the rate of 6,000 cards per hour.

The Interpreter
	The interpreter prints in type at the top of the card the
information punched in the card (see fig. 4-). A plugboard
permits the selection of the columns to be interpreted and the
position of the printed characters on the card. Interpreting
may be done on the front or the back of the card. Our
machine will interpret as many as 45 digits. The newer
models will interpret more columns, but have smaller type.
	The capacity of the machine is 6,000 cards per hour.
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The Calculation Control Switch

	The machines described above are all designed for computation
where each operation is done on many cards before
the next operation is begun. For the process of numerical
integration the different arthmetical operations must be done
in succession for the first step before proceeding to the next.
For this work it must be possible to change quickly from one
complicated setup to another.
	The calculation control switch controls the operatlOn of
the multiplier, tabulator, and summary punch so that a cycle
of different arithmetical operations can be performed in rapid
succession.

	4. OTHER MACHINES AND DEVICES NOT NOW IN USE
			AT THE BUREAU

Alphabetic Equipment

	For commercial accounting there is a complete line of
punches, tabulators, and interpreters, which handle literal

[Image of alphabetic tabulator record text includes:

"DISTRIBUTION OF PRODUCTIVE MATERIAL
					DATE May 31
Order  | MATERIAL 		      | Actual Amount | Standard Amount  | VARIANCE 
Number |Number| Name		      |		      |		   	 | Over | Under
3120   |   832| BRASS PLATE 	      | 416|50 	      |   400|00i 	 |161|50|
3120   |  1044| BRASS TUBING	      |  19|46 	      |    20|00	 |      |   |54
3120   |  1120| 6ⅹ9 INSULATING PANEL |  10|71 	      |    12|50 	 |      |  1|79
3120   |  1141| COILS 		      |  22|47 	      |    30|00	 |      |  7|53
       |      | 	 	      | 469|149•      |   462|50•	 |      |

8450   |   402| RH£OSTAT	      | 650|00 	      |   650|00	 |      |
8450   |   730| BOLTS AND NUTS	      | 210|00	      |   210|00	 |	|
8450   |   955| SO OT SWITCH 	      |  50|00        |    60|00	 |	| 10|00
8450   |  1201| NICKEL- SILVERWARE    | 125|40        |   120|00 	 |  5|40|
8450   |  1322| COPPER WIRE 	      |  12|25	      |    12|00	 |   |25|
       |      | 		      |l047|6S •      |  10S2|00•	 |      |

8462   |   832| BRASS PLATE 	      | 247|63        |   248|00 	 |	|   |37
8462   |  1044| BRASS TUBING          |  15|20	      |    16|00 	 |	|   |80
8462   |  1120| 6ⅹ9 INSULATING PANEL |   6|80	      |	    7|00	 |	|   |20
8462   |  1141| COILS		      | 125|14	      |   125|00 	 |   |14|
8462   |  1425| BUS BAR		      |   6|00	      |     6|00	 |	|
8462   |  1470| POTCHTIOMETER 	      |  19|20	      |    20|00	 |	|   |80 
84622  |  2055| PURCHASED PARTS       |  10|00	      |    10|00	 |	|
       |      |			      | 429|97•	      |	  432|00•	 |	|"

Caption: "FIG. 12. Alphabetic Tabulator Record" ]

as well as numerical data. Figure 12 shows information
printed from punched cards by means of an alphabetic tabulator.
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	This method of recording litera l in formation on cards
and printing partial lists from the cards is an efficient means
off furnishing current information from a file which is changing.
For instance in the U. S. Coast and Geodetic Survey the
descriptions of bench marks and their elevations could be
kept in such a file. Of course where the information is predominantly
numerical it is possible to code a small amount of
literal data and use the numerical machines.
	The alphabetic tabulator has many advantages over the
numerical tabulator for numerical work. For instance some
of the "carry-over" mechanisms of the counters are wired
through the plugboard so that the counter capacity may be
arranged to suit the problem. The chief disadvantage is the
increased cost.

The Collator
	The collator is a new type of machine which is important
for scientific work. It is really a mechanical filing clerk
which may be used to insert or remove cards from a file.
For example you have a file of 10,000 cards numbered consecutively
in columns 1-5. You have another set of "detail"
cards containing any 2,000 of the numbers from 1-10,000.
The collator will match each one of the detail cards with the
proper file card and remove it from the file. The collator
will later insert each of the file cards in its proper place.
	In the work of the Bureau the sorter has been used to
per form many of the functions of the collator, but of course
with less efficiency.

The Test Scoring Machine
	This machine is not a punched card machine, and it is not
a computing machine, but it is an interesting device which
might well have uses in science.
	It was designed to score the so-called "new type" examinations. 
Each question in the examination is answered by
placing a lead pencil mark in one of five positions. The completed
test sheet will thus have a pencil mark for each of the
100 questions. The test scoring machine compares these
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marks with similar ones the master sheet and instantaneously
gives the total score on a voltmeter. The score
may be computed in several ways wiht various weights assigned
to the different sections of the examination.
	The operation of the machine depends upon the electric
conductivity of a lead pencil mark rather than upon contact
made through a hole in a card.

Miscellaneous Devices
	With a high-speed listing device a standard tabulator will
list 6,000 cards per hour instead of 4,500.
	A tabulator may be equipped with special "fractional
wheels" for adding fractions or degrees, minutes, and seconds.
	We have seen above how the operation of the counters in
the tabulator may be controlled by an "x" punch in a given
column of the card. The "digit emitter" is a device which
will operate the counters according to the digits punched in
a given column. For example if our tabulator had a digit
emitter we could add in counter (1) information from all the
cards with a "2" in column 17, in counter (2) from those with
a "7" in column 17, etc. This device lends itself to a great
many useful tricks.
	Intermediate control provides an additional subclass in
the tabulation.
	The collating device on a sorter will deal cards for a bridge
game or similar purpose.
	The high-speed summary punch consists of a gang punch
attached to the tabulator instead of a key punch. It is
much faster, but of course more expensive.

[PAGE BREAK]

		PART II. MACHINE METHODS

			CHAPTER III

	    GENERAL CONSIDERATIONS OF THE
		PUNCHED CARD TECHNIQUE*

		1. INTRODUCTION

	The successful and efficient application of the machines
to a scientific problem requires careful consideration of human
and mechanical elements. The machines are not infallible,
so numerical checks must be provided. The machines are
not infinitely efficient, so it is necessary to take advantage of
all short cuts. The obvious direct method is frequently much
more costly than necessary. It is often possible to substitute
a little manipulation of the cards by the operator for long
machine processes. On the other hand the operator should
not be assigned tasks which are psychologically difficult to
perform without error. Part III of this book illustrates these
features, but a few of them may be summarized here in general
terms. The planning of an extensiYe program thus requires
a careful analysis of many factors in the light of all available
data and experience.
	It should be remembered that electric computing machines
of the keyboard variety are now so efficient that unless each
number required in the computation is to be used several
times the punched card method would not be efficient. The
advantage of the punched card method increases with the
number of uses to be made of each number concerned.

		2. NUMERICAL CHECKS

	The calculation must be so arranged that a check is provided
at each step. The machines are remarkable in their

* See also Practical Applications of the Punched Card Method in Colleges and
Universities, edited by G. W. Baehne. Columbia University Press, 1935.
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accuracy), but errors do occur. It must be remembered that
any failure in the entire machine may in a short time ruin
large quantities of work. It is not sufficient to have a check
at the end of a problem requiring weeks to complete, but each
operation must be checked before proceeding to the next.
	While it is not the purpose of this book to discuss tabulat-
ing technique which is part of the stock in trade of the average
operator, it is not out of place to mention some things which
are more important in scientific than in most other work with
the punched cards. In scientific computation a given card
may receive more varied treatment than in other work.
More different machines are used, the card may be run
through the machine with the back side forward, etc. For
these reasons the cards must be treated with great care and
the adjustment of the machines must be carefully checked.
Each new set of cards should be checked on the "card gauge"
to see that they are the exact size. The cards should be
carefully stored in files which hold them properly. Above
all the operator must insist that the machines be kept in
perfect adjustment so as to prevent "off punching." This
can be avoided if cards coming from all the machines are
frequently tested with the gauge. These are all standard
techniques but carelessness is fatal to successful work of the
kind described here.
	The first thing to be checked is naturally the punching of
the original cards from the manuscript. Three methods are
available for this: the use of the mechanical verifier, proofreading,
and the use of check sums; which is to be used
depends somewhat on the nature of the problem. If the
initial data are not extensive and the problem is long and
costly all three methods may be employed. If the data are
extensive and of a rough nature such as are used in many
statistical discussions there is no point in trying to find errors
which would not materially affect the result.
	The mechanical verifier in the hands of a conscientious
operator is probably the best means of checking ordinary
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numerical data. It docs not, however, detect double punching
in a given column.
	The cards may be listed on the tabulator and proofread
against the manuscript. This has two advantages. The
comparison may be made by someone who is familiar with
the data but is not a skilled key punch operator. Again, if
the original record is not written out in full, or if the information
about each case is not complete, the method of proofreading
is simpler.
	Checks which depend upon sums or upon sorting are
desirable because of the efficiency with which these operations
may be performed on the cards. Check sums detect misplaced
cards and errors of transposition. Another advantage
of check sums is that they may be carried on to the next step
in the work.
	Ordinary sorting is checked by looking through the holes
for small groups, and by passing a sorting needle through the
larger ones. This checks each sorting as it is made, but does
not insure against getting them out of order later. For
example in putting the cards in order on several digits, sorting
is done first on units, then on tens, etc. While sorting on the
tens it is essential that the units keep their proper place. If
the stacking mechanism of the sorter is out of adjustment one
card may stand on edge long enough for the next one to pass
under it.
	There is no direct infallible check on the listed numbers or
printed totals from the tabulator, but with proper precautions
the probability of such errors can be made very small
indeed. While listing with the tabulator, the numbers being
listed should also be added in the counters. The totals accumulated
in the counters may then be compared with similar
totals obtained by the tabulator. This insures that the numbers
were properly read from the cards. The failure of the
printing mechanism itself is rare and usually of such a nature
that it is easily detected.
	Tn tabulation it is desirable to punch the totals with the
summary punch. The addition of these summary cards then
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gives a perfect check on the tabulation. It is sometimes
desirable in checking work done on the tabulator to vary the
procedure so as to eliminate the possibility of the tabulator
making the same error twice. For instance listing and tabu-
lating are done at different speeds of the machine. The
counters may be interchanged. The counters may be wired
to change the signs of all the numbers. The cards may be
reversed from left to right so as to use different adding brushes.
The number of such precautions to be taken depends upon
how serious an undetected error would be.
	In checking summary cards. most totals may be checked
by straight tabulation. The parts, however, which come from
the group indicator require a special technique. The plugboard
is so wired that the tabulator does not break control
on each control group, but the counter which is adding the
group indication operates only from the impulses which in
the summary punching caused the break in the control. The
wiring for this operation is given in appendix B.
	In transferring information with the high speed reproducer,
the comparing brushes should be used to check the cards as
they are made. Tabulation of the cards is, however, a desirable
additional safeguard. It is possible for a card to feed
incorrectly into the reproducer, to be incorrectly punched, and
still to pass the comparing brushes.
	Ordinary multiplication is checked by running the cards
through again with the machine set for "check." As stated
in chapter II, section 3, this interchanges the numbers in the
multiplier and multiplicand counters, and accumulates the
sum of the products. This form of check is practically perfect
as far as the multiplication is concerned. It leaves the slight
possibility of an error in the transfer from right-hand components
to left. In the newer multiplying machines the
product formed in the check multiplication is compared with
the punched value on the card. This checks the individual
cards and locates the error at once, whereas in the older
method the fact that the sum of the products does not agree
simply tells that an error has been made.
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	When many products are formed with a constant multiplier
they should not be checked as above. The sum of the
products should equal the product of the constant multiplier
by the sum of the multiplicands. Usually not all the digits
of the product are punched on the card so the check on the
tabulator is not exact. However, when the multiplying is
being done additional digits of each product may be accumulated
in the summary products counter to check the last
required digit. Suppose for example you have a group of
500 cards with the multiplier and the multiplicand each of
eight digits and the products punched to ten. The summary
products counter has ten-digit capacity so we must drop the
three largest digits in order to leave room for the three additional
smaller ones. The larger digits in the sum of the
products are of course obtained by tabulating the punched
products.
	Where many small products not involving a constant
multiplier are to be checked we can form the sum of the
products by the method of section 5. When this method is
to be employed the cards should be sorted in order of the size
of the multipliers before multiplication. Small groups may
then be checked on the tabulator as they come from the
multiplier.
		3. SHORT CUTS IN MULTIPLICATION

	Multiplication is one of the most expensive operations, as
may be seen from the operation times given for each machine.
It is therefore desirable to find means of avoiding it whenever
possible. When the individual products are not required but
only their sum, the method of section 5 should be considered.
	When the multiplying punch is used it should be used to
capacity. For instance, it is possible to do more than one
multiplication at a time if the numbers are small. It is only
necessary to avoid one product overlapping the other. As an
example consider the formation from a, b of a^2, 2ab, b^2, where
a, b arc two digits each. Suppose for illustration a = 96,
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b = 84. The multiplication is as follows
[Equation roughly as follows:
		a      b
                /\    /\
		96000084  Multiplicand
		96000084  Multiplier
                --------
	       384000336
	      768000672
	 576000504
	864000756
	----------------
	9216016128007056  Product
	\  / \   /  \  /
	a^2   2ab    b^2  ]

Unfortunately this method will not work for two three-digit
numbers. It will work for a two and a three. Where a
great many such products are to be formed and a and b are
numbers each with not more than three digits it is possible
to sort out before multiplication those for which a is less than
three, then those for which a is 3 but b is less than three.
The remainder for which both a and b are three digits must
be done the long way. This technique may also be used as
follows: a + b in multiplicand and c + d in multiplier giving
ac, ad + bc, bd in the product.
	It should be remembered that when a product is being
formed the addition of a crossfooting operation usually
requires no additional machine cycles. The two operations
should therefore be combined whenever possible. It is possible
to perform for example a + b, c ⅹ d at the same time.
	Where large numbers of small products are required they
may be obtained with the sorter and gang punch. The cards
are sorted into groups according to the multiplicand and
multiplier. A small table gives the product for each group
and it is put on the cards with the gang punch. This method
pays only when most of the groups contain a considerable
number of cards.

[HEADER:	 	GENERAL CONSIDERATIONS 		31 ]

		4. SPECIAL SORTING TECHNIQUE

	For many calculations it is necessary to sort the cards
into groups with characteristics depending indirectly upon
the information punched in a given field. As an example of
such sorting consider the problem of sorting m cards, numbered
consecutively, into n groups according to the remainder left
after dividing the number by n. The general method requires
as many sortings as there are digits in m. Let us illustrate
with the case where there are 9,999 consecutive numbers and
the period n is 13. The four digits of the numbers are
punched in columns 1-4. Sort first on column one, forming ten
packs of cards. Then sort each pack on column two. The
first pack when sorted on column two gives 00, 01, 02, ··· 09.
The second gives 10, 11, ··· 19. When we come to pack 13
we combine it with pack 00, combine 14 with 01, and so on.
The sorting is now done on column three in the same manner,
and finally column four. The final result gives 13 packs of
cards. The first contains cards 00, 13, 26, ···, 9997, the
second 01, 14, 27, ···, 9998, and so on.
	The general method is not necessary for all periods. For
periods 2, 5, 10, only one sorting is required; for periods 4, 20,
25, only two; and for 8, only three.

5. THE FORMATION OF SUMS OF PRODUCTS WITHOUT
		THE INDIVIDUAL PRODUCTS

	In the method of least squares, and in the computation
of correlations we frequently require the sum of a great many
products, the individual products being unnecessary. Such
sums of products may be very efficiently formed with the aid
of the sorter and tabulator.*
	The process is known as "progressive digiting" and may
hest be explained by a simple example.

	* The Mendenhall-Warren-Hollerith Correlation Method. Columbia University
Statistical Bureau Document No. 1. Also The Use of Machine Factoring
in Multiple Correlation, A. E. Brandt, Journal of the American Statistical Asso-
ciation, September, 1928.
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	We require the sum of the following products
[equations
			1ⅹ2 =  2
			1ⅹ8 =  8
			2ⅹ3 =  6
			2ⅹ1 =  2
			2ⅹ5 = 10
			3ⅹ9 = 27

			5ⅹ6 = 30
			5ⅹ2 = 10
			6ⅹ7 = 42
			---------
			Sum = 137 ]

This sum may be found by addition only. Simply start at
the bottom and add the seven six times, the two and the six
five times, etc. The required additions are given below.
The last column contains in each line the sum of the numbers
in the adjacent column. In other words it is the "progressiYe
total" of thr numbers in the preceding column.

[ Equations
		  7=7		  7
		2+6=8		 15
		  0=0		 15
		  9=9		 24
	      5+1+3=9		 33
		8+2=10		 43
				---
				137  ]

The sum of the numbers in the last column gives the required
sum of the products by simple addition.
	The additions of the type outlined above arc done with
great efficiency on the tabulator. The cards are put in order
with the sorter, and the automatic control causes the printing
of a total after each group. Blank cards must of course be
inserted for the missing groups. If the summary punch is
available the progressive totals are punched on cards and
added with the tabulator; if not they must be added from the
printed tabulator record.
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	The great efficiency with the tabulator comes from the fact
that the product of one number by several others may be
treated simultaneously. For example, in the method of least
squares or for correlations we have several columns of numbers,
say a, b, c, d, ··· and we need the sum of products
[aa], [ab], [ac], [ad]. These quantities arc formed simultaneously
on the tabulator, [aa] in counter (1), [ab] in
counter (2), [ac] in counter (3), and [ad] in counter (4).
When the numbers in a, b, c, d are small, one counter could be
used for a and b, thus doubling the capacity of the machine.
	In the simple example above we have assumed a, b each
of one digit. Two digit numbers could be treated in the
same way with, at most, 99 groups; or the work could be done
in two steps with 9 groups each. For this latter method the
procedure is as above for the first step except that both digits
of the multiplicand arc added. When the first step is completed
the operation is repeated for the tens digit of the
multiplier. The results of the two steps are of course combined
by addition. This method may obviously be extended
to three-digit numbers.
	The following example illustrates the two-digit case.
		 1ⅹ 8 =   8
		 2ⅹ13 =  26
		 3ⅹ 7 =  21
		11ⅹ16 = 176
		12ⅹ 3 =  36
		23ⅹ 2 =  46
		------------
		   Sum = 313 ]

For the unit digits we have

[equations
		 2+ 7=  9		 9
		 3+13= 16		25
		16+ 8= 24		49
					--
					83  ]
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For the ten digits
	     	   2= 2		 2
		3+16=19		21
				--
				23 ]

Combining, 10ⅹ23+83=313.

		6. SPECIAL CARDS FOR BREAKING CONTROL

	In many problems it is desirable to form control groups
according to information not directly punched on the cards.
For instance we have on the detail cards a set of consecutive
numbers and wish to form control group of three each. This
may be accomplished by sorting into the detail cards a special
set of "control cards." These control cards must be punched
so they may be sorted into the detail cards and come at the
desired places. In addition to the punching required for sorting
the control cards into their proper places they must contain
an "x" in some given column. This "x" when used in connection
with a special wiring of the plugboard (see appendix B)
will cause only one break in control for each control card
instead of the normal two. An example of the use of control
card5 will he given in chapter V, section 3.

	7. SUGGESTIONS FOR CLERICAL WORK
	Where many computations are being made simultaneously,
some of which extend over long periods of time, general
schemes of identification and indexing are desirable.
	One method of identification is by means of colored cards.
In our laboratory manila cards are used for most of the computations
because they are the cheapest. Colored cards are
used in special cases. For instance when a card jams and
must be reproduced the reproduction is made on a salmon
colored card. The color aids in checking to see that the torn
card was properly reproduced and also shows up in the file
as a card which is more likely to contain error than the others.
	For a very extensive project it is desirable to have special
ruling and printing on the cards. Small files may be identified
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by ruling stripes on the top edges of the cards. The positions
of these stripes can indicate the number of the file. With
three stripes about five hundred files may be quickly identified.
The numbers should of course refer to a description of the
file. These stripes are easily applied to a trayful of cards
with a ruling pen. The marks are visible on the cards without
opening the tray, and also on the individual cards.

[Page Break]

		   CHAPTER IV

	THE CONSTRUCTION AND USE OF SPECIAL
	    TARLES OF TABULAR FUNCTIONS

		1. INTRODUCTION

	The tnbulator furnishes such an efficient means of con-
structing tabular functions that every computational
problem must be examined to see to what extent special tables
maybe used in its solution. In some cases the construction
of several manuscript tables to be used in the usual way may
furnish the most economical solution of the problem. In
others the special tables should be prepared in the form of a
punched card catalogue to be used in connection with other
machine work.
	The interval adopted for the tabulation depends upon the
number of times the table is to be used. In some extreme
cases it may pay to tabulate so that interpolation is unnecessary.
In this case "critical" tables are convenient. The
most common tables are those in which the interval of tabula-
tion is chosen just small enough that linear interpolat10n with
first differences is legitimate. Where comparatively few
values of the function are to be used, the interval may be
chosen so large that second and higher differences are required
in the interpolation. In view of the efficiency with which
the tables may be built up on the tabulator, smaller intervals
will in general be taken than if the tables had to be constructed
by more laborious means.

		2. SUMMATION WITH THE TABULATOR

	The tabulator will automatically build a table from the
differences, and print argument, function, and difference with
one run of the machine. The method is extremely efficient
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and is self-checking. Consider as an example the construction
of a table of logarithms beginning at 5000.

[Numerical table
Arg 		 Log 		Diff
5000 		69897
				  9
5001 		69906
				  8
5002 		69914
				  9
5003 		69923
				  9
5004 		69932
				  8
5005 		69940
				  9
5006 		69949  
]

	In building this table on the tabulator three counters are
used; the first for the argument, the second for the function
and the third for the difference. The first counter builds up
the argument as a progressive total by adding "1" from each
card, and the second builds the function by adding successive
values of the difference. The third counter is used only for
the printing of the individual differences and is cleared or
reset after each line is printed.
	The differences 9, 8, ··· are punched on a series of cards
by means of the key punch (unless they are the result of a
previous calculation with the punched card machine). With
the gang punch a "1" is punched in each of these cards.
We also prepare one card with the first value of the argument
and of the function, in this case 5000, 69897. The tabulator
is wired so that 5000 is added in the first counter and 69897
in the second. Further the "1" from the difference cards
must be wired to the unit position of counter (l), and the
field containing the difference must be wired to counter (2)
and also to counter (3). The automatic control is set to
break control on each card. Counters (1) and (2) are set for
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"total, non-reset," and counter (3) for "minor total and
reset." The tabulator record would appear as follows.

[table of figures
5000 	69897

5001 	69906 9

5002 	69914 8

5003 	69923 9 ]


The vaIue of the difference which is printed on any given line
e \11 uc 
is the one used in obtaining the value of the function on that
line from the preceding. It is possible to set the machine to
print the difference on the odd line of the table if desired.
	In order to check such a table it 1s only necessary to
check the first and the last values of the function and argu-
ment and the values of the differences in column three. Thus
each line contains ten digits, only one of which requires
personal attention. With the subtraction tabulator we can
take care of positiYe and negative values of the function and
difference.
	If the table is to be used in the form of a punched card
file we use the same process but punch summary cards to
match the above printed record. There is one important
difference, however, and that is that we want on a given card
the argument, the function, and the value of the difference
which is to be used in interpolating between the value of the
function on the card and the value corresponding to the next
vaJue of the argument. This may be accomplished by building
the table from the bottom rather than from the top. The
printed record is then inverted but the cards can be sorted
into the normal order.
	In the above example the differences were taken from a
table which had already been printed. In the construction
of a new table the differences are obtained without first
computing the values of the function. Methods of obtaining
the differences required to build up a table with a smaller
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tabular interval than that of a given table will be discussed
in chapter V.
	The process described above may be applied successively
to build a table from second or higher differences. In this
case if the summary punch is not available it is necessary to
punch the results of the first summation on cards in order to
use them as the differences for the second process. Where the
table is built from second differences counter (4) would be
used to print the second difference in the final table.
	In printing the higher differences of a table the arrangement
depends upon the interpolation formula to be used in
taking values of the function from the table. For instance it
may be desirable to use one of the well-known formulae
requiring horizontal instead of diagonal differences. Before
building the table of the function from the first difference it
is possible to move all the values in the second difference
column up or down the table by means of the high-speed
reproducer. For example, suppose we have a file of cards
each containing a first and a second difference and it is
desired to transfer to each card the value of the second difference
from the preceding card. For this operation feed "B"
of the high speed reproducer is used. The plugwires from
the brushes are so arranged that numbers read by the comparing
brushes are punched in a new field on the card.
Checking may be done on the tabulator by adding the numbers
from the two fields. If the cards are run through the
machine in reverse order we get on each card the following
instead of the preceding second difference. If the interpolation
formula requires the sum of two adjacent second differences,
this sum may be obtained by crossfooting with the
multiplying punch.

3.  THE USE OF TABLE IN THE FORM OF PUNCHED CARDS

	For many problems the ease of taking values of a function
from a table by the punched card method is a very important
factor. This is especially true where the previous and following
steps in the calculation are done by means of punched
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cards. This method is most efficient when a great many
values of the function are to be takcn from a given table.
	Let us consider a file of cards which constitutes thee table;
each card containing an argument, a function, and a difference.
There is also a large number of detail cards each containing
an argument for which the function is required, and 
an "x" punch in some given column. Suppose further that
the function is tabulated in the table with three digits in the
argument, and the values of the argument given on the detail
cards have six. The three digits of the argument on the
detail card which correspond to the three digits of the argument
of the table cards must be in corresponding columns.
	The first step is to sort the table cards and the detail
cards together in order of the three digits of the argument,
and with the table card preceding the detail cards in each
case. Thus for instance if the argument of the table runs
from 1 to 999 we have, after sorting, table card 001 followed
by all the detail cards 001; then table card 002 followed by
detail cards 002, and so on.
	The entire set of cards is now passed through feed "B" of
the high-speed reproducer. The wiring is arranged so that
information from each card is transferred to the next, pro-
vided the card being punched has an "x" in the prescribed
column. Thus the values of the function and of the difference
are punched on each detail card following it, but nothing is
punched on the table card. The machine may be wired to
punch only the cards without an "x," in which case the "x"
should be on the table cards. The "x" on the detail cards is
preferable because in case the machine fails to pick up an
"x" the table card will not be spoiled. The comparing
relays are used to check the punching. The tabulator may
of course be used for additional checking.
	In planning a given computation provision should be made
that the field on the detail cards containing the digits of the
argument which match those of the file should be in the
proper columns, either when the card is direct or reversed.
If the cards must be reproduced to obtain this result, efficiency
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demands that an attempt be made to have the previous card
completely used before reproducing.
	The collator is much more efficient than the sorter for
combining the file and detail cards. In addition, the question
of matching the fields of the argument does not arise, for the
plugboard permits sufficient flexibility in the wiring.
	One run on the sorter is required to separate the file from
the detail cards. The column with the II x" is suitable for this.
	The interpolation of the function for the last three digits
of the argument is now done on the detail cards. This is
performed on the multiplier as an operation of the form
A ± B ⅹ C. The special "x" brush on our machine makes
it possible to do all interpolation of this form without separating
the plus ones from the minus.
	Where interpolation with second differences is required
the procedure is as above except that two runs on the multiplier
are in general required. The particular values of the
first and second differences which are to go on the card depend
upon the particular interpolation formula to be used (see ch.
V). Stirling's formula is most convenient for this purpose
but requires smaller third differences for a given accuracy
than Bessel's. Stirling's formula for this purpose would be
written:

[formula
f_n = f_0 +  n [ (δ'_-1/2 + δ'_1/2)/2 + n (δ_o'')/2] ]

This requires two runs through the multiplier of the form
A ± B ⅹ C provided each card contains n; f_0, δ_o''/2, and
(δ'_-1/2 + δ'_1/2)/2. Bessel's formula

[formula
f_n = f_0 +nδ'_1/2 + n (n-1)/4 (δ_o'' + δ_1'') ]

is not as simple to apply but might be used in special cases,
for instance when the third differences are just too large for
Stirling's formula.
	In view of the methods here suggested for interpolation it
will be seen that the considerations leading to a choice of the
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tabular interval are different from those used in connection
with hand work. For example in interpolating with first
differences it is only necessary to secure an interval where it
is legitimate to do so, the size of difference being unimportant
to the machine. For example, Peters' table of natural functions
designed for ordinary use is given for each thousandth
of a degree, whereas for the card catalogue we need only
hundredths.
	In some cases the table should be tabulated for such a
narrow interval that part of the interpolation may be done
with the sorter and the gang punch. ( See chapter IX.)
	A list of the tables in the form of card files available at
the Bureau is given in appendix A.
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				CHAPTER V

		INTERPOLATION, MECHANICAL QUADRATURE,
			AND ALLIED SUBJECTS

			1. GENERAL REMARKS

	It is not the purpose of this chapter to discuss all the
available formulae for interpolation, numerical integration,
and numerical differentiation, but rather to indicate how the
punched card method may be applied to advantage. The
reader 1s referred to the many books on the subject to the
collections of formulae, and to the tabulated values of the
various coefficients required in the different formulae. The
following short list is suggested as a beginning:

(1) Steffensen, J. F. Interpolation. Williams and Wilkins.
	1927.
(2) Comrie, L. J. Interpolation and Allied Tables. Reprinted
	from the Nautical Almanac for 1937. His Majesty's Stationery Office.
(3) Comrie, L. J. On the Construction of Tables by Interpolation.
	Monthly Notices of the Royal Astronomical
	Society, volume 88, April 1928.
(4) Pearson, K. Tracts for Computers, No. II. On the
	Construction of Tables and on Interpolation, Part I,
	Uni-Variate Tables. London, Cambridge l'ni\'ersity
	Press, 1920.
(5) Bower, E. C. Systematic Subdivisions of Tables. Lick
	Observatory Bulletin 467. 1935.

	All the well-known methods of interpolation and mechanical
quadrature are equivalent to finding a polynomial of
degree n which represents n + 1 consecutive values of the
function, and then computing the required quantity from this
polynomial. In the formulae of the Lagrangian type the result
is expressed as a function of the n + 1 successive tabular
values of the function; in the difference formulae one value
of the function and the values of n successive differences are
employed.
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	The justification of this process may be established for
known functions which have continuous derivatives of a
certain order.* For other functions the computer must be
satisfied by the apparent convergence of the calculations
themselves or from a knowledge of the general behavior of
the functions considered.
	The choice between the two methods depends not only
upon the problem itself but upon the computation which
prercdrs and follows it. For example we might use a difference
method instead of the Lagrangian where the differences
had been formed in the previous work.
	The principal disadvantage of the difference method is
that the work of forming the differences may in some cases
be greater than the complete interpolation by the Lagrangian
method. On the other hand the differences form a check of
previous work which should not be overlooked. Another
advantage of the difference method is that the differences are
smaller than the functions and hence the necessary multiplication
is easier.
	The relative merits of the various methods have been
frequently discussed, for example see references (3) and (4)
at the beginning of this chapter. In addition to the classical
arguments should be included those which are introduced by
the punched card method. By way of illustration, consider
two examples from chapter XII where the decision was easy to
make.
	For the numerical integration of the differential equations
of planetary motion the difference method was chosen because
the check afforded by the differences was urgently needed,
and because two extra places in the computation were possible.
In this case a single error would invalidate weeks of
subsequent computation. The differences checked not only
the integration itself, but the perturbations which had been
previously computed. In the second place, the multiplying
punch has only eight digit capacity which would limit the
Lagrangian methods to eight figures. The difference method
gives the possibility of ten.

* See page 2 of reference (1).
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	In the interpolation of the perturbations for the integration,
however, the Lagrangian method is used. This interpolation
is usually to halves and quarters instead of to tenths
so that the efficiency of the summation process is not so great
as for tenths. The previous calculations did not require dif-
erences so it would have been necessary to compute them
Just for the interpolation. Finally the integration which is
the next step in the calculation furnishes an automatic check
not only of itself but of the perturbations as well.

			2. NOTATION
	In this and other chapters we shall use a uniform notation
for tabular functions and differences. A given table consists
of successive values of the argument and the corresponding
tabular values of the function. The quantities formed by
subtracting each successive value of the function from the
following one are called the differences of the first order.
The differences of the second order are formed by similar
operations on the first differences. Conversely we may think
of the first differences as the first summation of the second
differences. As an example tabulate f(x) = x^3 for values of
the argument from -2 to +4.

[Table
	  |	      |   Differences
	  |	      |______________________
Argument  | Function  | 1st | 2nd | 3rd | 4th
__________|___________|_____|_____|_____|____
  -2 		-8 
			+7
  -1 		-1 	      -6
			+1 	    +6
   0 		 0 	       0	   0
			+1 	    +6
  +1 		+1 	      +6 	   0
			+7 	    +6
  +2 		+8 	     +12 	   0
		       +19 	    +6
  +3 	       +27 	     +18
		       +37
  +4 	       +64   ]
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In order to avoid ambiguity it is necessary to assign a value
of the argument to the difference of each order as well as to

the function. This may be done by reading along a hori-
zontal line in the table, or along either of the two diagonals.
We shall use all three methods and use the symbols δ, Δ, ∇ in
designating them. They are referred to as horizontal, de-
scending, and ascending differences. Superscripts will be used
to indicate the order, and parentheses to indicate the argument.
Thus in the above table

[equations
f(+  2)=+ 8
δ'(+1½)=+ 7	 Δ'(+2)=+19	 ∇'(+2)=+7
δ''(+2)=+12	Δ''(+2)=+18	∇''(+2)=+6 ]

	In dealing with a small region of the table we can refer to
a few successive, values of the argument as -1, 0, +1, +2,
···; the idea being that we are not interested in the particular
yalue of the argument but only in the sequence of the values
of the function and differences. Interpolation formulae are
usually written in this form. In this case the argument is
written as a subscript rather than in parentheses. We list a
few of the most common of such formulae:

Bessel's,

[equation
f_n= f_0 + nδ'_½ + n (n-1) δ''_0 + δ''_1
		     _____ _____________
		      2! 	  2

		+ n (n-l)(n-½)
		    __________ δ''_½ + ···.
			3!			]



Stirling's,

[equation 
f_n = f_0 + n (δ'_-½ + δ'_½) + n^2
	      ______________   ___ δ''_0
		     2 		2
			+ n (n - 1)(n + 1) (δ''_-½ + δ''_+½)
			    ______________ _________________ + ···,
				   3! 		   2


Lagrange's four-point,

[equation
	f_n = L_-1 f_-1 + L_0 f_0 + L_1 f_1 + L_2 f_2  ]
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where

[equation
	L_-1 = -n (n-l)(n-2)
		  __________ ,
		       6

	L_0 = (n+l)(n-l)(n-2)
	      _______________ ,
		     2

	L_1 = - (n+1)n(n-2)
		___________ ,
		     2

	L_2 = (n+l)n(n-1)
	      ___________
		   6

	L_1 + L_0 + L_1 + L_2 = 1. ]

As special cases of Lagrange's four- and six-point formulae we
have for interpolation to halves:

[equation
f_½ = 1/16(-f_-1 + 9f_0 + 9f_1 - f_2),

f_½ = 1/256(f_-2 - 25f_-1 + 150f_0 + 150f_1 - 25f_2 + 3f_3). ]


Simpson's rule for mechanical quadrature:

[equation 
	∫^{2}_{0} f(x)dx = 1/3ω(f_0 + 4f_1 + f_2). ]

	The use of Stirling's and Bessel's formulae has been illustrated
in chapter IV. The use of formulae of the Lagrangian
type will be discussed in section 3. The question of systematic
subtabulation will be treated in section 5.

	3. INTERPOLATION BY THE LAGRANGIAN METHOD

	The adaptation of the punched card method to formulae
of the Lagrangian type is well illustrated by the interpolation
to halves by the four and six point formulae. The former is
particularly efficient. Interpolation to fourths is done by
two successive applications of the method. We now proceed
to describe interpolation to halves.
	Examination of the formulae of section 2 shows that in the
6-point formula only three multiplications are required for
each value of the function, and these arc easily checked by
addition since they involve a constant multiplier. The four-
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point formula requires only one multiplication, for the multi-
plication by 9 can be performed as plus ten minus one times
the function. In the six-point formula the factor 1/256 is
combined with the 3, 25, 150; but in the four-point the mult1-
plication by 1/16 is reserved till after the add1t1on of the
components. The remainder of the work consists of sorting
the cards into proper groups and adding on the tabulator.
	Consider the case where we have three different functions
of a similar nature tabulated for the same values of the argument.
We shall designate them as f(x), f(y), f(z), the subscript
indicating the argument being understood. An ex-
ample of three such functions would be the three components of
the force acting on a planet. For each value of the argument
there would be a card punched in five fields: (1) argument,
(2) f(x), (3) f(y), (4) f(z), (5) identification of the problem
and of the particular step in the computat10n. In the interpolation
of perturbations for a planet we reserve columns 1-10
for fields (1), (5); columns 11-20 for field (2); columns 21-30
for field (3); and columns 31-48 for field (4).
	The first step after the cards have been punched and
verified is to make one or two additional sets of cards with the
high-speed reproducer. For the six-point formula we require
two duplicate sets in addition to the original. For the four-
point formula we require one set like the original except that
f(x), f(y), f(z) have been moved one position on the card,
which amounts to multiplication by ten. While making the
cards on the reproducer, the new identification is applied by
the gang punch. If the original set was identified by a zero
in a given column we use "1" for the first set made on the
reproducer and "2" for the second. In addition an "x"
must be punched in those groups which are to be subtracted
according to the formula, which is in both cases those identified
as "1." In the four-point case this "x" will be used to
add and the "no x" to subtract.
	For the six-point formula f(x),f(y),f(z) are now multiplied
by 3/256 for the "0" set, 25/256 for the "1" set, and 150/256
for the" 2" set. This is of course done as rate card multiplica-
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tion, and checked on the tabulator (see ch. III, sec. 2). The
products should be so punched that when the card is reversed
the product will occupy the same columns as the original
number.
	The cards are now ready to be arranged in groups according
to the proper formula. For the six-point formula we
must make up a control group containing a card for each of
six consecutive values of the argument, the first and the last
being "0," the second and fifth being "1," and the third
and fourth being "2" cards. In addition the tabulator must
change the sign of the "l"s by means of the "x" punch.
Further the tabulator must give the argument of the result
as that of f_0 plus 1/2 the tabulator interval. The groups will
be made up with the sorter, and special "control" cards (see
ch. III, sec. 6) will be used to give the proper control group
and to give the proper argument of the result.
	The control cards constitute a set which serves for all
interpolations of this type. The arguments are punched in
the same columns as on the detail cards. In another field of
the control card is the argument which is to go with the
interpolated result. Thus for the six-point formula, where
detail cards in the group have arguments 4, 6, 8, 10, 12, 14 the
control card must sort into the deck as a three but must
print the resultant argument as 9. One set of control cards
with a range of one thousand in the sorting argument can carry
in different fields the proper interpolated argument for interpolation
with an interval of 2, 4, etc., by the four, six, or eight
point formula. The necessary stock of cards for this method
is thus very small compared to that of any other method.
	The sorting into control groups is done as follows for the
six-point. All the cards "0," "1," "2" are sorted together
according to the argument, and separated into packs according
to the period six (see ch. III , sec. 4). They are then
further subdivided according to the "0," "1," "2." There
are thus eighteen packs of cards which are placed on a table
in three columns of six rows each. The three columns
contain, respectively, the " 0," "1," and "2" cards, and the
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six rows contain czrds whose arguments, after division by 6,
have residuals 0, 1, 2 , 3, 4, 5. If now we take up from the
table all those whose position is 1ndicated by "a" in figure 13A
we have all the cards necessary for obtaining every sixth
interpolated value in the table. These cards are placed in the

[table
		  A		  B
		a b c		d f e
		b a c		d e f
		b c a		e d f
		c b a		e f d
		c a b		f e d
		a c b		f d e

FIG. 13 Scheme for Collecting Cards in Lagrangian Interpolation

sorter together with the set of control cards, and sorted in
order of the argument. This set is now ready for tabulation.
When this set has been tabulated and the results recorded
with the summary punch, we proceed to do the same for the
cards whose positions are indicated by "b"; and finally by "c."
It is noted that the scheme of figure 13A gives half the
required groups. The others must be gotten by replacing the
cards on the table and collecting according to figure 13B.
The above sorting may be done at one time for any number
of sets of functions having the same arguments, a final sorting
being required to separate them.
	For the four-point formula are required four consecutive
values of the function "0," and two values of "1." After
the summary cards have been obtained the interpolated values
must be multiplied by 1/16.

		4. THE FORMATION OF DIFFERENCES
	The formation of successive differences of a function is
the first step in most methods of interpolation and numerical
integration. We give here two alternative methods of forming
these differences with a tabulator. In the first method the
differences of all orders are derived for a given value of the
argument before proceeding to the next, while in the second
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the differences of the first order for all the arguments are
derived before proceeding to the second. The first method
is necessary in numerical integration, and may be desirable
where high-order differences are involved, especially if the
table contains errors which are to be found from the differences.
The second method is more efficient in that it does
not require the full attention of the operator.
	In either method all the counters of the tabulator and a
large part of the card should be used at each step. If there
are several functions of a similar nature for each value of the
argument, they should be taken together. For instance, in a
great deal of astronomical work we use the same card for f(x),
f(y), f(z). If there is only one function and the table is long,
it is possible to break the table into three or four parts so
that the values of the function for arguments separated by
100 or 1000 may be considered together. The first counter
of the tabulator is used for the argument, the order of the
difference, and the identification of the particular set of functions
being differenced. The second, third and fourth counters
carry the corresponding values of the differences of the
three functions. In some cases it might be desirable to use
the four counters for four functions and to apply the argument,
order of the difference, and identification by means of the
high-speed reproducer.
	In the first method we wish to find all the ascending
differences ∇^i for a given argument before proceeding to the
next. This is accomplished by finding first the difference of
the highest order and then obtaining the others by successive
additions. By definition we have for any diagonal in the
table:

[equations:
f_0 - f_-1 = ∇_{0}'
∇_{0}' - ∇_{-1}' = ∇_{0}"
· · · · · · · · · · · · ·
∇_{0}^{n} - ∇_{-1}^{n} = ∇_{0}^{n+1} ]

Adding, we find

[equations
∇_{0}^{n+1} = f_0 - (f_-1 + ∇_{-1}' + ∇_{-1}'' + ··· + ∇_{-1}^{n})  ]
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Thus at each step of the calculation we subtract all the num-
bers along the diagonal from the next value of the function.
We then to this result the same previous numbers one
after the other (progressive totals). This method is obviously
suitable for the punched card technique. It is only necessary
to work out a simple routine which will make the proper ad-
ditions and subtractions, and to punch summary cards and to
clear the counters at the proper times. The manipulations
of the cards and control switches must be made simple.
	The tabulator is wired so that when the cards arc tabulated
in the direct position the numbers are added (algebraically),
and in the reverse they are subtracted. For this
purpose an "x" must be punched in one column of all the
cards. The dial switches are set for "total, non-reset," and
the counters are cleared by tabulating the last summary
card reversed. Minor control is wired to the column which
indicates the order of the difference, and the corresponding
column on the reversed card must have a zero in each card.
Three colored cards, the first punched zero in the control
position, and the other two punched with the digit of the highest
order difference, facilitate the routine.
	The operator's routine is:

	(1) Turn switch "S. P. Minor" to "off" position.
	(2) Tabulate: f_-1, f_-1, ∇_{i}', ∇_{i}'', ···, ∇_{-1}^{n} reversed; f_0,
three control cards direct.
	(3) Place in tabulator hopper: ∇_{-1}^{n}, ∇_{-1}^{n-1}, ···, ∇_{-1}', f_-1.
	(4) Turn switch "S. P. Minor" to "on" position, and
start tabulator.

	In the above routine we have started with all the values of
f_0 punched on cards, and have in addition formed a duplicate
card for each value to simplify the clearing of the counters.
	In the process of numerical integration (see ch. VIII) we
form five or six differences and two summations for f(x),
f(y), f(z) by an adaptation of the above method. For that
work a special circuit in the calculation control switch provides
for the clearing of the counter, thus eliminating the
additional card f_0.
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	In the second method we wish to subtract (algebraically)
the numbers on each card from those on the following. This
requires that the cards be tabulated in pairs, and that the
first card of each pair be subtracted and the second added.
	The necessary wirings for both these operations are de-
scribed in appendix B. Section D of figure 21 refers to the
tabulation of the cards in pairs. Section A gives the wiring
for subtracting numbers on one set of cards from those on
another. The subtraction requires an "x" in alternate cards.
	Tabulation of the function cards with the above wiring
gives the first difference for every other value of the argument.
When this operation is completed we remove the first card of
the pack and run them through again to get the remaining
arguments. For the second run we must change the wiring
of the "x" so that in this case also we subtract the first
card of each pair from the second.
	For inspection of the differences it is necessary to sort
the two sets of differences together and list them.
	The "x" on the function cards may be put on by the key
punch operator or with the gang punch. The "x" on the
first differences, needed to get the second differences, may be
applied during the summary punch operation.

		5. THE SUBDIVISION OF TABLES
	In the construction of tables where each value of the
function involves long computation, the function is computed
for a few key values and the ones between are inserted by
interpolation. We have discussed the Lagrangian method
of interpolation to halves and quarters. Where more than
three values of the argument are to be inserted the difference
method has a great advantage. We have seen in chapter IV
how tables may be easily built by summation, once the differences
are known. The problem is then to find a method for
determining the differences in the subdivided table from those
of the original. To be suitable for the machine the method
must be continuous, i.e., it is not desirable to make any
numerical adjustment by "inspection." On the other hand
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a machine method could be efficient even if more arithmetic,
within reason were involved than for a hand method.
	Two methods have been suggested for use with punched
card machines and I shall refer to them as the "end figure" 
method, and the "superfluous difference" method. Both
methods are important contributions, but from the nature of
the problem, it seems that a still better method could be
devised.
	In the end figure method the end figure of each interpolated
value of the function is computed. These end figures
are differenced to an order high enough to give smooth dif-
ferences, say second or third. Where the differences run
smoothly the entire value of the difference can be filled in
when the end figure is known. With these values of the
highest order differences, those of the lower orders are built
up by summation. The end figures of the interpolated values
of the function are formed on the tabulator. Cards are selected
from the files depending upon the various order differences
of the original table, and these are added on the
tabulator. (See reference (3), p. 43.)
	The end figure method has the advantage of giving the
interpolated functions and differences just as they are to be
printed in the final table. It has the two disadvantages:
(1) a collection of thousands of special cards is required for
the interpolation; (2) a large amount of clerical work is necessary
in selecting the proper cards from the files, and in
differencing the end figures. This clerical work could of
course be adapted to the machines if sufficient interpolation
were contemplated to warrant it.
	The superfluous difference method builds the subdivided
table from differences of one order higher than that required
in the original table! In addition to the high order of the
differences, extra decimal places are required for the summation.
Another disadvantage is that, having used extra decimals
in the differences, they will not be the values of the
differences which should be printed in the table. After the
functions have been properly rounded, the end figures of the
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differences must be adjusted according to the rounded end
figures of the functions. From reference (5) at the beginning
of the chapter is given below a table showing the differences in
the subdivided table in terms of those in the original. It
must be remembered that in building tables from any given
difference by this method, the accuracy of the interpolation is
only that of one order less.

[Table
			Differences
_________________________________________________________________
Third	|Fourth	|     Fifth	|     Sixth	| Seventh 	|
________|_______|_______________|_______________|_______________|
  +.055	|-.0165	|  -.0086625 	|  +.00329175	| +.0017007375	|
  -.045 |+.034  |  +.0257125 	|  -.0134145 	| -.0085259625	|
    0   |-.0165 |  -.0247875 	|  +.0202555 	| +.016953475   |
    0   |  0    |  +.0078375    |  -.0134145    | -.016716525   |
    0   |  0    |    0 		|  +.00329175 	| +.0081802875  |
    0 	|  0 	|    0 		|    0 		| -.0015910125  |
    0 	|  0	|    0		|    0		|   0		|
    0	|  0 	|    0 		|    0 		|   0		|
    0 	|  0 	|    0 		|    0 		|   0		|
    0 	|  0 	|    0 		|    0		|   0		|
________|_______|_______________|_______________|_______________|
]

For comparison of the end figure method with the superfluous
differences method let us assume that interpolation
must include fifth differences, the sixth being negligible. In
the superfluous difference method we must build with sixth
differences using eight extra decimal places. In the end
figure method the third differences in the subdivided tables
are of course only about one thousandth of those in the
original, and no extra decimals are used in the summation.
With ten-digit counters in the tabulator and six-digit functions,
the superfluous difference method would work well for
fourth differences. The alphabetic machine has greater
counter capacity and could be used for this purpose. A
numerical tabulator with special carry-over from one counter
to another would be useful.
	In all these difference methods one is impressed by the
number of cards which must be punched by the machines and
used only once or twice. In England there is on the market a
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machine which will transfer from one counter to another.
This "rolling total" machine would be much more efficient for
forming differences and summations since it would not be
necessary to punch cards. This observation does not apply
to the process of numerical 1ntegration for then the tabulator
is used for other operations in the meantime and the numbers
could not remain in the counters.
	The newer models of the multiplying punch offer possibilities
in the building of tables by summation. In these
machines it is possible to transfer from one counter to another
and to punch progressive totals. With slight modification one
of these multipliers might become a remarkable "D1fference
Engine."
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			  CHAPTER VI

		NUMERICAL HARMONIC ANALYSIS
			AND SYNTHESIS

1. INTRODUCT!ON

	A great many functions occurring in science may be represented
with sufficient accuracy by harmonic series containing
a limited number of terms:

[equation
(1) F(x) = a_0 + a_1 cos x + C_2 cos 2x + ··· +a_n cos nx,
	       + b_1 sin x + b_2 sin 2x + ··· + b_n-1 sin (n - l)x.]

The problem of analysis is the determination of the a's and b's
from special values of F(x). The problem of synthesis is the
determination of special values of F(x) from a given series.
	The usefulness of the method depends to a great extent
upon the ease with which the analyses and syntheses can be
carried out. Both are done with great efficiency by the
punched card method where the series have many terms or
where many functions are to be analyzed or synthesized at
the same time.
	The numerical work in either case consists of the solution
of linear algebraic equations, and the special Yalues of x are
chosen in such a manner as to simplify the numerical work.
If we choose special values of x as follows:

[equation
x_i = ih, 	h = π/n ]

we shall have the well-known formulae:

[equation
(2) F_i = ∑_{k=0}^{n} a_k cos ikh + ∑_{k=1}^{n-1} b_k sin ikh

							(i = 0, 1, 2, ···, 2n - -1),
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n ε_j a_j = ∑_{i=0}^{2n-1} F_i cos ijh 	(j = 0, 1, ···, n),


ε_j = 2 for j = 0, n
    = 1 for j = 1, 2, ···, n - 1,
2ll-1

nb_j = ∑_{i=0}^{2n-1} F_i sin ijh(j = 1, 2, ···, n - 1).

	For a great many problems it is more convenient to separate
the sines from the cosines and consider ea ch separately.
	Let
			F = C + S,
where

	C(x) = a_0 + a_1 cos x + ··· + a_n cos nx,
	S(x) =       b_1 sin x + ··· + b_{n-1} sin (n - l)x.

Then for values of x chosen in pairs, a and 2π - α, we have

		F(α) + F(2π - α) = 2C(α),
		F(α) - F(2π - α) = 2S(α),

and the formulae of analysis and synthesis become:

(3)  	C_i = ∑_{k=O}^{n}a_k cos ikh		(i = 0, 1, ···, n),

	S_i = ∑_{k=1}^{n-1}b_k sin ikh		(i = 1, 2, ···, n - 1),

  n ε_j a_j = ∑_{i=0}^{n}ε_i C_i cos ijh	(j = 0, 1, ···, n),


       nb_j = ∑_{k=1}^{n-1}ε_i S_i sin ijh	(j = 1, 2, ···, n - 1),

	ε_i = 1 	for	i = 0, n
	    = 2		for	i = 1, 2, ···, n-1.


	The analysis or synthesis is easily executed with the help
of a table of sin (ijh) and cos (ijh). We give here such a
table for n = 6.
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			h = π/6

		d = sin h = ½ 	e = cos h = ½√3

			Sin
_________________________________________
i \ j | 0   1   2   3   4   5   6  | ε_i
______|____________________________|_____
  0   | 0    0   0   0   0   0  0  |
  1   | 0   +d  +e  +1  +e  +d  0  | 2
  2   | 0   +e  +e   0  -e  -e  0  | 2
  3   | 0   +1   0  -1   0  +1  0  | 2
  4   | 0   +e  -e   0  +e  -e  0  | 2
  5   | 0   +d  -e  +1  -e  +d  0  | 2
  6   | 0    0   0   0   0   0  0  | 
______|____________________________|_____
  n   |      6   6   6   6   6

		Cos
_________________________________________
i / j |  0   1   2   3   4   5   6  | ε_i
______|_____________________________|____
  0   | +1  +1  +1  +1  +1  +1  +1  | 1
  1   | +1  +e  +d   0  -d  -e  -1  | 2
  2   | +1  +d  -d  -1  -d  +d  +1  | 2
  3   | +1   0  -1   0  +1   0  -1  | 2
  4   | +1  -d  -d  +1  -d  -d  +1  | 2
  5   | +1  -e  +d   0  -d  +e  -1  | 2
  6   | +1  -1  +1  -1  +1  -1  +1  | 1
______|_____________________________|____
 nε_j | 12   6   6   6   6   6  12  |
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	Each computation consists of multiplying given numbers
by the numbers from these tables and adding the results.
A great deal of the work consists of adding the same numbers
in different combinations, which is done most efficiently by
punched cards.

		2. MACHINE METHODS FOR SMALL SERIES
	The machine technique is best described by means of an
example. Consider the analysis of a cosine series when n = 6
and where each Yalue of C_i is given to seven significant
figures. Suppose there are 99 such sets of functions to be
analyzed at the same time. For each set are given C(0°),
C(30°), C(60°), C(90°), C(120°), C(l50°), C(180°) or C_i where
i = 0, 1, ···, 6.
	Since our tabulator has four counters it will be possible
to add corresponding quantities for three different functions
in three of the counters and to use the fourth for identification
of the series and the operation. Nineteen cards will be required
for each three sets of analyses.
	The giyen functions are punched on the cards as follows:
Columns 1-2 contain the number of the series being analysed
i.e. 1, 4, 7, ···, 97. Column three contains i = 0 1, ···, 6.
Columns 5-11 contain C_i for the series whose number is
punched in columns 1-2, i.e. numbers 1, 4, 7, ···, 97. Column
12 is punched with an "x" if C_i is negative, and with a "y"
if positive. C_i for numbers 2, 5, 11, ···, 98 are similarly
punched in columns 14-21, and for numbers 3, 6, 12, ···, 99
in columns 23-30. A zero is punched in columns 4, 13, 22.
When the cards have been punched and checked they are
passed through the reproducer to give the necessary dupli-
cates and to punch quantities needed for the tabulation. It
is necessary to make duplicate sets for i = 1 - 5. To identify
the duplicates a "y" is punched in column three. All
cards with i = 1 and 5 (including duplicates) are now sorted
on columns 1, 2, and the same numbers gang punched in
columns 80, 79. Finally we gang punch on the cards the
following pattern:
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[Table
_________________________________________________________________
Value of    |			Column on Card
    i       |____________________________________________________
            |31| 32| 33| 34| 35| 36| 37|  |44| 45|46|47|48| 49|50
____________|__|___|___|___|___|___|___|__|__|___|__|__|__|___|__
    0 	    | 0| 1 | 2 | 3 | 4 | 5 | 6 |  |  |   |  |  |  |   |
    I 	    | 0|   | 2 |   |x_4|   |x_6|  |  |x_5|  |  |  | 1 | 
    1y 	    | 0|   |   |   |   |   |x_6|  |  |   |  |  |  |   |
    2       | 0| 1 |x_2|x_3|x_4| 5 | 6 |  |  |   |  |  |  |   |
    2y 	    | 0|   |   |X_3|   |   | 6 |  |  |   |  |  |  |   |
    3	    | 0|   |x_2|   | 4 |   |x_6|  |  |   |  |  |  |   |
    3y	    | 0|   |x_2|   | 4 |   |x_6|  |  |   |  |  |  |   |
    4	    | 0|x_1|x_2| 3 |x_4|x_5| 6 |  |  |   |  |  |  |   |
    4y	    | 0|   |   | 3 |   |   | 6 |  |  |   |  |  |  |   |
    5	    | 0|   | 2 |   |x_4|   |x_6|  |  | 5 |  |  |  |x_1|
    5y	    | 0|   |   |   |   |   |x_6|  |  |   |  |  |  |   |
    6	    | 0|x_1| 2 |x_3| 4 |x_5| 6 |  |  |   |  |  |  |   |
____________|__|___|___|___|___|___|___|__|__|___|__|__|__|___|___ ]

This pattern insures three things in the later tabulation: (a)
that the proper cards are included in each group, (b) that
they are tabulated in the required direct or reverse position,
and (c) that each has the proper sign.
	The numbers on cards 1 and 5 are now multiplied by 2e
as rate card multiplication, and checked on the tabulator
(see ch. III, sec. 2). The products are punched so that on
reversing the card from left to right the product will occupy
the same position as the multiplicand. No signs are indicated
in the products. Three wirings of the multiplier are necessary
because three functions occur on each card.
	The tabulator is wired to control on columns 1, 2, 31.
Counter (1) is wired as follows: positions 10, 9, 8 to columns
1, 2, 31 of card through group indication; position 7 to card
count. Counters (2), (3), (4) are wired to columns 4-12,
13-21, 22-30 respectively.
	Counter (1) always adds. Counters (2), (3), (4) add and
subtract from "x" in columns 12 or 69, 21 or 60, 30 or 51
when there is no x in 31, and subtract and add when there is
an x in 31. This is accomplished by the method described
in chapter II, section 3.
	The summary punch is wired to punch positions 10, 9, 8
from counter (1) in columns 1, 2, 3 on the card. The quan-
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tities from counters (2), (3), (4) are punched in columns 4-12,
13-21, 22-30 with " x" in columns 12, 21, 30 for negative numbers.
This summary card therefore resembles the original
cards.
	All the cards are now sorted in order according to columns
1, 2. Tabulation with summary punch on minor control gives
12a_0 for each function. The results are checked by tabulation.
	The cards are now sorted on i. The cards for i = 0, 1 y, 4,
6 are put in the sorter in the usual manner, but those with
i = 1, 5 are reversed. All are now sorted on columns 1, 2.
The plug wires on the tabulator going to column 31 of the
adding and control brushes are moved to column 32 and everything
is ready to give 6a_1.
	The above process is continued according to the schedule
until i = 6 has been completed.
	We have thus obtained 12a_0, 6a_1, 6a_2, ···, 12a_6. Rate
card multiplication by 1/12, 1/6 gives a_0, a_1, ···, a_6. These
results may be listed or used for further computation. If the
results are to be synthesized they may be run reversed through
the reproducer to obtain a set of cards of the same form as the
original.
	The modifications necessary for the synthesis or the ana)y-
sis of a sine series are obvious and need no further description.
	The checks outlined above are very rigorous, and the
general appearance of the results is usually a good test of the
work. However, in very important cases the results in analysis
may be checked by synthesis and vice versa.
	The time required for this example should be:

		Punching and checking. . . . . . . . . . . 3 hrs.
		Reproducing and gang punching . . . . . . ½ hr.
		Multiplying . . . . . . . . . . . . . . . 2 hrs.
		Tabulating and summary punching. . . . . 2 hrs.

		3. MACHINE METHODS FOR LARGE SERIES

	In the above example the machines were used advantageously
on series of only a few terms in those cases where
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many analyses or syntheses were to be done at the same time.
For series of many terms it may pay to use the punched card
method for only a few series at a time.
	Consider for example the cosine series with thirty harmonics
 where numerical values of the coefficients are given
and it is required to compute C(x) for the special values of
x  = nh, h = 6° , n = 0 , 1, 2, ···, 30. Four figure accuracy
will be assumed, though this can be anything up to eight.
	For this method a set of 228 prepunched cards are used
for each three series to be synthesized. They contain the
required cosines and also 30 columns with numbers for sorting
and for g1ving the proper algebraic sign of the products. The
a's are punched by hand onto other cards. Each "a" card
will contain corresponding a's for three different series. Negative
values are indicated by an "x" punch. The "a" cards
and the prepunched cards are sorted together and passed
through the multiplier to give the product of each a by the
required cosine; the "a" cards are the rate cards and the
prepunched cards are the detail cards. The multiplication
is checked as it is completed by means of the sum of the
products (without regard to sign). The cards must go
through the multiplier three times, once for each of the three
series. In addition to the above check on the absolute values
of the products, the sign values are checked by one run on
the tabulator. The cards are now sorted according to the
value of "n" and added on the tabulator. If the tabulator
has a summary punch the results may be checked by addition.
	We have set aside 30 columns in the prepunched cards for
sorting, control, and indication of the signs of the products.
These 30 columns correspond to the desired values of "n."
To add the products for a given value of "n" the cards are
sorted on that column and all the cards which are not to be
added are rejected by the sorter. An "x" in this column is
used to operate an "x-distributor" on the tabulator in such a
manner that terms with negative a's will be added, and those
with positive ones will be subtracted.
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	The time required to synthesize 48 series will be:

		Reproduce the prepunched cards (16
		  sets) . . . . . . . . . . .  . . . . . 1 hour
		Punch and check the a's   . . . . . . . 3 hours
		Multiply. . . . . . . . . . . . .  . . 15 hours
		Sort, tabulate, and summary punch . . 10 hours.

	The methods described aboYc in this and the previous
section require yery little in the way of prepunched cards,
and hence leave great freedom in the choice of the particular
scheme to be adopted. Should occasion arise where one
particular scheme is to be used day in and day out, the stock
of prepunched cards could be efficiently increased. For example
Dr. Linus Pauling of the California Institute of Technology
in his work on crystal structure has synthesized
thousands of series of thirty terms each. It has proved
efficient for his work to prepare a set of about 65,000 cards
containing all the possible combinations of coefficients. He
can then select the proper cards by hand and add the components
on the tabulator.
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			CHAPTER VII

		THE MULTIPLICATION OF SERIES

			1. THE PROBLEM

	In many physical problems we have two functions each
expressed as a harmonic series with a limited number of terms
and we require the harmonic series which represents the
product of the two functions. The formation of these product
series is efficiently done by the punched card method when
the series are long or when many product series can be formed
at one time. The application of the method here described
to power series is obvious.
	We shall discuss the case where the argument of each
term is a linear function of the independent variable involving
several fundamental frequencies. This includes of course the
ordinary Fourier series. Three kinds of series occur: sines,
cosines, and exponentials.
Let

(1)		S(x) = ∑b_i sin α_i,
		C(x) = ∑a_i cos α_i,
		E(x) = ∑c_i e^{iα_i},

where

	α_i = (jp + kq + lr + ···)x,
		j, k, l, ··· = 0, ± 1, ± 2, ···,
		P, q, r, ··· = given numbers.

The product of any one of these by it self or by any one of the
others may be required. In forming the product of two series
it is necessary to consider each term of the first in connection
with each term of the second. The two coefficients must be
multiplied and the harmonic functions combined according to
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the trigonometric formulae:

[equations
(2) 	2 sin x sin y = cos (x - y) - cos (x + y),
	2 sin x cos y = sin (x - y) + sin (x + y),
	2 cos x cos y = cos (x - y) + cos (x + Y),
		 e^x e^y = e^{x+y},
		 2 cos x = e^{ix} + e^{-ix},
		2i sin x = e^{ix} - e^{-ix},    ]


	Each pair of terms in the product of two trigonometric
series gives a term with the sum of the arguments and one
with the difference. Exponential series require only the sum,
but since exponential series are usually used for the computation
of functions where both the function and its conjugate
are required, the difference of the arguments is needed also
for that case.
	The multiplication of one harmonic series by another thus
consists of forming all the sensible products of the coefficients
and the corresponding sums and differences of the arguments,
and collecting these terms according to the resultant arguments.
	The machine technique will vary with the particular problem,
depending upon the number of terms in a series, the
number of fundamental periods in the arguments, and other
considerations. We shall describe the procedure for two problems
of widely different nature. One or the other of these
two methods should be adaptable in most cases.

		2. EXTENSIVE SERIES
	Let us consider the case where we have S(x) and C(x)
each with five hundred terms, four fundamental periods in
the argument, and ten significant figures in the largest coefficient.
Suppose for simplicity the largest coefficient is 1.0.
Assume further that extreme care is to be exercised to prevent
error. In such a series the integers j, k, l, m would not
exceed ± 10.
	Each term of each series is punched on a card showing
identification of the series, the argument, and the coefficient.
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Since Lhe argument has more than one component it would be
troublesome to examine each argument to see if it is positive
or negative. Instead we decide arbitrarily on the proper
order of the components of the argument and treat it as
positive or negative according to this scheme. Thus we may
consider an argument as positive if m is positive and negative
if m is negative. If m equals zero, l decides, and then k and
finally j. This leaves only the constant term. The coefficients
are punched as true figures with an "x" to indicate
negative numbers. The integers j, k, l, m however are made
positive by the addition of say 50. Thus the argument
(p - 2q + s) would be indicated on the card as 51485051.
When the cards have been punched and verified they should
be listed in the order of the argument for proofreading.
	The cards of each series are now sorted according to the
numerical value of the coefficients with 1 ⅹ 10^{-10} first and
1.0000000000 last. It may seem at first thought that this
would require ten complete sortings, but this is not the case.
The chance that a number is zero in the units and tens position
and not in the remaining positions is only one in one
hundred. Hence after the second sorting large groups may
be eliminated from further sorting by looking through the
holes. After the third sorting the total number may be still
further reduced. The final sortings involve only a few cards.
	We now print with the tabulator two lists for each series.
The first has four columns: (a) consecutive numbers from 1
to 500, (b) the argument, (c) the coefficient with sign, (d) the
coefficient without sign. The second list is just like the first
except that progressive totals are printed in columns b, c, d.
	It is now necessary to set a rejection limit for the products,
say 1 ⅹ 10^{-11}. In other words all products smaller than
1 ⅹ 10^{-11} will be rejected. We must examine each pair of
coefficients to sec if the product is within the rejection limit.
Consider the largest term of C(x); it will give a sensible
product for each term of S(x). The next largest term in C(x)
may permit the rejection of all terms in S(x) smaller than
say 2 ⅹ 10^{-10}. This rejection limit is computed for all terms
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in C(x) greater than √10 ⅹ I0^{-6}. It is more convenient to
examine the remaining terms in the reverse order, i.e. com-
pute the rejection limit for the largest. term in S(x) and so on
until we come to the term in S(x) whose rejection limit is
 √10 ⅹ I0^{-6}. Thus all the sensible products may be found
by considering about 100 terms.
	With the help of the computed rejection limit we form all
the required pairs of terms. A new card is made with the
high-speed reproducer for each pair of terms. These are
made as follows: the first card of C(x) must be paired with
each card of S(x), the second with all greater than the rejection
limit, etc. The set S(x) is placed in the feed A of the
reproducer; and new cards, preceded by the proper card from
C(x), are put in feed B. Thus by a combination of gang
punching and reproducing we obtain cards each of which has
two arguments and one coefficient. The coefficient which is
common to the whole group need not be punched on the
cards but may be punched on a special red rate card to go
with the group. By a continuation of this process we obtain
about 100 groups of cards each with its rate card. Since the
capacity of the multiplier is eight digits, those groups which
haYe 9 and 10 digit multipliers and eight digit multiplicand
are done in two groups. The few terms with large multiplier
and multiplicand are done by hand. In a series with reasonable
convergence there may be about 20,000 cards altogether.
The reproducer is of course used to check the punching.
	In addition to the check of the reproducer itself we may
now get an additional check from the tabulator. Each group
of cards is tabulated to give the sum of the arguments, the
card count, the sum of the coefficients, and the sum of the
absolute values of the coefficients. These totals are now com~
pared_ with the values obtained by subtracting the proper
Iines in list number 2 above. This checks everything connected
with the formation of the group.
	The multiplication is now performed for each group and
checked as rate card multiplication. Extra digits arc accumulated
m the summary products counter as described in
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chapter III, section 2. The check sum from the summary
products counter is compared at once. This sum checks the
multiplication but not the signs. Additional checks are obtained
during a later tabulation.
	The value of the sum and difference of the two arguments
on each card is evaluated by means of the sorter and high
speed reproducer. For those cases where the value of j, k, l
or m which is to be added and subtracted is zero, the desired
result should be applied at the time the duplicate cards are
made; it being only necessary to punch the same number for
the argument and for the sum and difference of the two
arguments. The largest groups of course have mostly zeros
in the argument, and hence most of the work is done at the
time the duplicates are prepared. For the others we have a
set of master cards giving the sum and difference of any pair
of integers from minus ten to plus ten. These cards are
sorted together with the detail cards for j of the two arguments,
and run together through the gang punch. The wires
are then moved on the reproducer and the cards are re-sorted
to give the corresponding results for k and then l and m.
	When all the sums and differences of the arguments have
been formed the cards are sorted to test for the proper sign
of the product, and are then sorted into the ongmal groups.
These groups are now tabulated to give the sum of the differences
of the arguments, the sum of the sums of the arguments,
the sum of the products of the coefficients, and the
sum of the absolute values of the products. The sums and
differences of the arguments are checked exactly. The other
two sums give fairly accurate checks on the punching of the
products.
	The cards are now sorted in order of the sum of the arguments
and tabulated with control on the sum of the arguments.
The results are recorded on the printed record and on tne
summary cards. The resultant argument of course comes
through the group indicator. The checking is done by tabu-
lation of the old cards and the summary cards. We have
outlined checks for such work in chapter III, section 2.
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Here every precaution should be taken with the checking,
such as reversing the cards in the tabult1tor, and interchanging
the counters.
	The cards are now sorted according to the difference of
the arguments, and then separated into two groups accord ing
to whether the cards were made with C(x) or S(x) as common
multiplier. In other words the difference of the arguments
may be (x -y) or (y - x) of equation (2). The summary
cards are obtained for the two groups separately, and by
means of the switch on the summary punch the signs of the
coefficients are changed for the set (y - x).
	We now have three sets of summary cards with arguments
(x + y), (x - y), (y - x). They must be combined, and
at the same time those having negative arguments may be
combined with those having positive. Since the resultant
terms are sines, the sign of the coefficient must be changed when
the sign of the argument is changed. We sort (x - y) and
(y - x) according to the sign of the argument. For negative
arguments we punch the complement of each component with
the aid of the sorter and reproducer. These new arguments
are punched so as to occupy the same position when the card
is reversed. The cards (x + y), (x - y), (y - x) are now
combined, those with positive argument direct and those with
negatie argument reversed, and sorted together according
to the argument. The cards are tabulated and the results
summary punched to giYe the resultant series 2 C(x) S(x) expressed
as a sine series.
	In forming the square of a series only one stack of cards
is involved, and fewer products are necessary. Each card
need be taken only with itself and those following. This
gives all the distinct products but not all the products of the
required series. Since there are many more cross products
than square terms the easiest thing to do is to divide all the
squares by two. This gives the square of the original series
for sines or cosines, and half the sHquare for exponentials.
	!he following identification has been used to identify the
various cards used in forming a product series. Two columns
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are required for identification. The first contains:

0 = first original series,
1 = second original series,
3 = :cd rate _cards made from first series and the correspondmg
detail cards made from the second series,
4 = blue rate cards made from second series and the corresponding
detail cards made from the first series,
5 = summary cards (x - y), (y - x),
6 = summary cards (x + y),
7 = product series.

The second contains a zero except:

	(1) Groups 3, 4 for the large digits of those groups whose
multipliers must be separated into two parts.
	(2) Group 5 for (y - x).

In these cases a "1 " is punched instead. In order to secure
uniformity those product terms which must be done by hand
are punched entirely on the (03), (04) cards, and no (13) or
(14) cards are used for these terms. Also, as soon as all the
products have been punched and checked the square cards
are removed and new ones made on green cards with half
the coefficients. If the sum of the products is formed first
for (x - y) these terms all come in one group with arguments
50505050.

		3. SMALL SERIES

	The case outlined in section 2 is an example of the largest
series that are apt to be encountered. One such series is
sufficient to justify the use of the punched card method. For
small series it would be necessary to do a considerable number
at one time, unless the previous and following computations
are also to be done on the punched cards.
	Consider the case where we haYe one hundred pairs of
series to be multiplied. Each series has one fundamental
frequency and nineteen Yalues of i from plus 9 to minus 9.
The coefficients have eight digits, and products less than
1 ⅹ 10^{-8} arc to be neglected.
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	The cards are punched so that each card con tains_ the
argument, the number of the series, and the coefficient.
Usually in such a case as this some of the series will each be
used to multiply several series, but only one set of cards
will be punched by hand.
	When the cards have been punched they should be listed
with the tabulator and totals obtained for the sum of the
arguments, the sum of the coefficients, and the sum of the
absolute values of the coefficients.
	Duplicate sets should be made with the high-speed reproducer
so that for each product series required there will be a
multiplier series and a multiplicand series, each carrying a
number to indicate the number of the product series. Colored
cards should be used for the multiplier series. The multiplicand
cards should have an "x" in a given column. In
case only a few duplicate cards are required this step could
have been avoided by modifying the initial hand punching.
All the sets of cards should now be listed as for extensive
series, and the totals punched on summary cards.
	For these small series the constant multiplier groups are
so small that the hand method of rejecting the negligible
products would not be efficient. It must be done by succes-
sive sortings, on all the series at one time.
	The colored multiplier cards for all the series are now
sorted according to the numerical value of the coefficient, and
the corresponding "rejection limit" for terms of the multiplicand
is applied with the gang punch. This latter quantity
need be accurate to only one or two significant figures, but
should be rounded so that additional terms will be included
in each group rather than that significant terms will be
omitted. It should be applied so that on reversing t he cards
it will occupy the same field as the original coefficient. The
argument and number of the series should also be transferred
to the other end of the card so that when the card is reversed
they will occupy the same columns as before.
	The multiplier and multiplicand cards are now sorted
together according to the numerical values of the multi-
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plicand cards. The multiplier cards are placed in the sorter
in the reverse position followed by the multiplicand cards
in the direct position. When the cards are all in order according
to the coefficients, they are then sorted into 100 groups
according to the number of the series. In each of these series
the multiplier cards (red) are followed by the necessary multiplicand
cards, and the corners of the multiplier cards stick
out so the top group in each case may be easily lifted from
the pile. These top groups are collected and the red cards
separated from the manila ones by sorting on the column with
the "x." The manila cards are now put through the reproducer
to form duplicates, and the duplicates are then placed
with the red ones. The cards which were used to form the
duplicates are replaced in the reproducer and to them are
added all those groups obtained by lifting the second colored
cards from the original series. These second colored cards
are attached to the second set of duplicates as before. This
procedure is continued until all the necessary duplicates
have been formed. Finally there will be nineteen packs of
cards containing red and manila cards. These are all placed
in the sorter, in each case red cards followed by the corresponding
manila ones. They are all sorted in order of the
numbers of the series. By means of the gang punch we now
transfer the argument and the coefficient of each red card
to the following manila ones. The "x" punch is used to
prevent the transfer to the red cards. The red cards are then
separated from the manila ones by means of the sorter. We
now have all the cards of a given series ready for the multiplier.
There should be perhaps forty or fifty in each series.
	The multiplication for each series is performed, and the
products punched with nine digits. The sum of the products
in the summary products counter may be checked as each
group is completed by comparing it with the product of the
sum of the absolute values of the multipliers by the sum of
the absolute values of the multiplicands for the series.
	With the tabulator the sum of the punched products may
be formed. Each of these should be compared with the
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product of the sum of the multipliers by the sum of the
multiplicands. If the check products.were formed on puncheu
cards the comparison can be made with the tabulator.
	All the cards are now sorted on the two arguments and
the sum and difference applied with the gang punch.
	From here on the procedure is the same as for large series.
	The machine time for 100 series is as follows:

		Original punching and verifying. . . . . 2 days
		Listing and summary punching check
		   sums. . . . . . . . . . . . . . . . . ¼ day
		Sorting and gang punching duplicates. . 1 day
		Multiplying . . . . . . . . . . . . . . . 1 day
		Sorting, tabulating, and summary
		   punching. . . . . . . . . . . . . . . 1 day

	In this example the checks would have been exact had
we included a card for each product regardless of size, but
this would haYe required about forty thousand cards instead
of four or five thousand, and the multiplying time would
have been over a week instead of one day.
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			CHAPTER VIII

	THE NUMERICAL SOLUTION OF DIFFERENTIAL
			EQUATIONS

	The method described here for t he numerical solution of
differential equations will work for the vast majori ty of cases
arising in physical science. It includes all simultaneous systems
of ordinary differential equations where t he accelerations
are functions of the time and of the coordinates and t heir
derivatives.

	 Let x, y, z = the dependent variables,
		t = the independent variable,

		X = d^{2}x, and similarly for Y , Z.
		    ______
		    dt^2

The differential equations are

(1) 		X = X (x, y, z; x', y', z'; t)

and similarly for Y, Z. Higher derivatives may also be
included. The solution is given as a table of values of x, y, z
for equidistant intervals w of t. It is giYen by the well-known
expression

(2)   X = ''X + 1/12 X - 1/240X'' + 1/1951 X^{Ⅳ}  + ···,

and simila rly for y, z. Here ''X , X , X '', X^{Ⅳ}, ··· are numbers
taken from an ordinary difference and summation table. ''X
is a column of numbers of which X is the second difference,
X'' is the second difference of X, and so on. The six arbitrary
constants of the solution may be taken as the values of ''X,
''Y, ''Z for any two successiYe tabular valucs of t.
	All the numbers in equation (2) refer to a given Yalue of
t, i.e., the values of X'', X, ''X, ··· lie on a horizontal line in
the difference ta ble. ln the notation of chapter V they a re
''X, δ^2 ''X, δ^4 ''X,
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	The solution is built up step by step as follows. As soon
as a new value of X is computed another diagonal line in the
difference and summation table may be found. This gives
the next value of ''X, required by equation (2), but not the
required values of X, X'', X^{Ⅳ}, ···. These latter must be
extrapolated in the table with sufficient accuracy for this
purpose. X must be extrapolated one tabular interval, X''
two, X^{Ⅳ} three, etc. If the interval of tabulation is small
enough this extrapolation may always be done correctly. The
difference between the extrapolated value and that computed
later is evident from the difference table.
	The values of x, y, z just deriYed from equation (2) are now
inserted in equation (1) to compute the next value of X. The
values of x', y', z' for the same Yalue of the argument are also
required, and are computed from

(3) 		x' = ('X) - 1/12(X') + ···,

where ('X) and (X') are the means of the numbers lying just
below and above the horizontal line of ''X, X, ···.
	The mechanical procedure is as follows. Special tables of
functions required for the computation of X (x, y, z; x', y', z'; t)
are prepared in advance in the form of files of punched cards,
each card containing argument, function, and difference. All
the multiplication and interpolation of tabular functions required
for equations (1), (2), (3) are done on the multiplying
punch, and all the addition and formation of differences and
summations are done on the tabulator. In all this work with
the tabulator we assign the first counter to the argument and
to the identification of the particular step in the calculation, and
counters (2), (3), (4) to the corresponding numbers pertaining
respectively to x, y, z. The work of the operator consists of
passing cards from one machine to the other, the selection of
a card from the files of tabular functions, and the inspection
of the highest order differences printed. The calculations
are always repeated in a given sequence which is quickly
learned. All switching of the machines to change from one
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operation lo another is performed automatically by means of
the calculation control switch. Thus the machines are always
ready to perform the next computation when the last has
been completed.
	The extrapolation of X, X'', ··· for equation (2) requires
only a few seconds, and the same is true of equation (3).
Take for example the extrapolation of X one step in advance.
It is only necessary to add together all the differences down
the diagonal from the right to left. This gives the extrapolated
quantity with an error of the lowest difference neglected.
In practice where fifth and sixth differences are employed
the highest should always be less than one hundred.
The factor 1/12 as well as the small factor contained in the
function X should be sufficient to give 1/12 X correctly.
	The calculation control switch contains a row of electric
contacts each of which is operated by a rotating cam. The
cam is a circular fiber disk which is notched at various points
around the circumference. A series of about twenty of these
disks are attached to a common shaft to form a sort of player
piano roll. When this roll is rotated from one position to
the next the various contacts open and close according to the
notches in the disks. The circuits from the contacts are used
to operate the various control switches on the tabulator and
multiplier, and also a number of multicontact relays which
effectively change the wiring of the plugboards. Each step
in the integration consists of a certain number of distinct
machine operations which always come in the same order.
Hence in order to have the machines ready for each operation
it is only necessary to rotate the roll from one position to the
next, one complete revolution corresponding to a complete
step in the integration. One roll serves for all equations of a
given form, and a new one could be prepared in a few hours.
	As an example of the method consider the integration of
the equations of planetary motion, which has been performed
for many different planets. A detailed discussion of this
problem will be given in chapter XII, but we shall describe
it here in general terms.
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	For planetary motion equation (1) becomes
		X =  -k^{2}x + P,
		      ______
		        r^3

where 		P = k^2 ∑m_i (x_i - x) 
			     _________
			     ρ_{i}^{3}

		x_i = a function of t alone,

		p_{i}^{2} = (x_i - x)^2 + (y_i - y )^2 + (z_i - z)^2,

		m_i = small constants,

		r^2 = x^2 + y^2 + z^2 .

	Since m_i are small we use approximate values of x, y, z and
compute P in advance. P, for our purposes, is thus a f unction
of t alone which is computed and punched on a card for
each value of the argument. The equation then becomes

X = - k^{2}x + P (t).
      ______
       r^3

A table with argument r^2 and giving 1/r^3 is also made up in
advance. This works for all planets. During the process of
integration, r^2 is computed and 1/r^3 is found by interpolation in
the table.
	For this integration ten significant figures are used in the
calculations, about 60 points are used for each revolution of
the planet, and fifth and sixth differences employed. The
mtegration proceeds at the rate of about 12 points per hour.
	It might be added that even if the tabulator and summary
punch only were available, it might in some cases be efficient
to evaluate X with an ordinary electrical computing machine,
and use the tabulator for all the differences, summations, and
estimat1ons.
	The Punched Card Method may be used in connection
with other methods of solving differential equations. For
example, chapter XI illustrates solution by the substitution
of harmonic series.

[PAGE BREAK]

	PART III. ASTRONOMICAL APPLICATIONS

		CHAPTER IX

	THE CALCULATIONS INVOLVED IN THE
	CONSTRUCTION OF A CATALOGUE OF
	PHOTOGRAPHIC STAR POSITIONS

		1. INTRODUCTION

	The most important catalogue of this type is that being
constructed by Dr. Frank Schlesinger and Dr. Ida Barney at
the Yale Observatory, and known as the Yale Zone Catalogue.
It is described here because of the part played by the punched
card method in performing the laborious calculations which
are inherent in the problem. The punched card method was
not adopted until the program had been under way for several
years, and for this reason the methods previously employed
were changed only where the advantage of the new method
was very great. Had the punched card method been aYailable
at the outset it would no doubt have been used to fuller
advantage. As it is, a saving of fifty per cent has been introduced
in the computations which are done on the cards, and
as these constitute about half the complete computations, a
saving of twenty-five per cent has been introduced in the total
cost of the calculations. The calculations now performed on
the cards illustrate many of the methods of part II of this book.
	The Yale Zone program consists of the reobsenation of
the relative positions of all the stars in the so-called A. G.
Catalogue. This catalogue was completed about forty years
ago, so the comparison of the old and new positions gives very
accurate values for the relatiYe motions. The positions in
the A. G. Catalogue were obtained by visual measurements
with a meridian circle and constitute the combined work of
several observatories for many years. The photographic tech-
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nique which is being used in the new program is so efficient
that the calculations now become an important part of the
total cost of the project.
	In the photographic method a series of photographs of the
sky are taken and the positions of the stars on the photo-
graphic plates are measured with a measuring engine. From
these measures of rectangular coordinates (X, Y) on the plate
must be computed the spherical coordinates on the sky.
These spherical coordinates are known as right ascension and
declination (α, δ). The computation must be made for each
star as it appears on two overlapping plates, and with seven
significant figures. In order to facilitate the identification of
the stars in the measuring engine the process is reversed, i.e.,
the spherical coordinates given in the A. G. Catalogue are
converted into the rectangular coordinates as they appear on
the plates. Then, since the relative motions of the stars in
fifty years are minute, the small differences between the observed
and computed X and Y may easily be converted into
the small differences in α and δ. The application of these
small corrections to the old α and δ, give the new α and δ, and
the long conversion does not have to be made again. The
most costly calculation in the process is thus the conversion
from spherical to rectangular coordinates. The well-known
formulae are:

(1)		 X = K           tan ∆α
		       ____________________________
 		       tan δ' sec ∆α sin δ + cos δ'

		Y = K  tan δ' sec ∆α cos δ - sin δ
		       ___________________________
                       tan δ' sec ∆α sin δ + cos δ

X, Y are the rectangular coordinates of a star on the plate,
δ' is the declination of the star,
δ is the declination of the center of the plate,
∆α is the difference in right ascension between the star and
   the center of the plate, and
K  is a constant for the entire program.

The average number of stars on a plate is about 1,000 and the
plates overlap so that each star appears on two plates.
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	Special tables have been printed for effecting this transformation 
but when Lhe plates arc large, and accurate positions
are required, the use of the tables becomes more expensive
than the direct calculation.
	The second computation which is now done on the cards
is the computation of the precession for each star. Right
ascension and declination are measured with respect to a
moving coordinate system, so that the changes in these
coordinates consist not only of the small relative motions of
the stars themselves but also of the motion of the axes which
is much larger. In order to find the so-called "proper motions"
it is therefore necessary to allow for the "precession."
For the purposes of the present catalogue the precession in α
and δ for a given star may be computed from the following
formulae:

(2) α = α_0 + p_{α}t + s_{α}t^2,
    δ = δ_0 + p_{δ}t + s_{δ}t^2,
    α_0, δ_0 = right ascension and declination to be corrected,
    α, δ = corrected right ascension and declination,
    t = number of years for which precession is to be applied.
      In this case t is the same for a very large number of
      stars. p_α, p_δ, s_α, s_δ, are functions of α_0, δ_0. They have
      been computed for each star and are given in the A.
      G. Catalogue.

The computation is simple but must be made a great many
times. It is complicated by the fact that astronomers still
use the sexagesimal system for star positions. α, α_0 are given
in hours, minutes and seconds while δ, δ_0 are given in degrees,
minutes and seconds.
	In order to make the results of such an observing program
immediately available the observations are made in zones of
declination. In the A.G. Catalogue, for instance, the various
zones were constructed at different observatories. The new
photographic catalogue is all being prepared at Yale, but the
various zones are published as soon as they are completed.
The following shows the work done in this program and the
part played by the punched cards:
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[TABLE:
_____________________________________________________________
	  |Plates|	  |	    |	
 Zone	  |Taken |Measured|Published|Punched Card Computation 
__________|______|________|_________|________________________
+55 to +60| 1916 |  1926  |  1930   |
+50 to +55|1916  |1916,22 |  1925   |
+25 to +30| 1928 |  1928  |  1933   |
+20 to +25| 1928 |  1928  |  1934   |
-2 to + 1 |1914  |  1914  |   1926  |
	  |	 |	  |	    |
+ 9 to +20|1940  |   1940 |	    | Coordinates, Precession
- 2 to + 9|1936  |   1939 |	    | Coordinates, Precession
- 2 to -10| 1933 |  1938  |	    | Coordinates, Precession
-10 to -14| 1933 |   1934 |   1939  | Coordinates
-14 to -10| 1933 |   1934 |  1940   | Coordinates
-20 to -30| 1933 |   1934 |         | Coordinates, Precession
__________|______|________|_________|________________________
]

		2. THE PUNCHING OF THE CARDS

	For each star in the A. G. Catalogue a card is punched
containing all the pertinent information: the star number;
the right ascension in hours, minutes, seconds and hundredths;
the declination in degrees, minutes, seconds and tenths; the
precession and secular variation in α and δ; the mean epoch,
and the magnitude. The last two items are for a purpose
other than the computations here described. In addition to
the customary verification of the cards, totals have been
formed for some of the columns of the catalogues and these
totals used for checking.

		3. THE COMPUTATION OF THE PRECESSION

	The precession for an interval of fifty years is to be applied
according to formulae (2). The necessary calculations for
each star are:

			α = α_0 + 50 P_α ± 12.5s_α,
			δ = δ_0 ± 50 p_δ ± 12.5s_δ,

where α_0, δ_0, P_α, p_δ, s_α, s_δ are given on each card and

			  2^{S}.2 < p_α < 3^{S}.8,
			           p_δ < 20''.1,
				   S_α < 0^{S}.0250,
				   S_δ  < 0''.600.

The results are required to 0^{S}.001 in α and 0''.01 in δ.

[HEADER: CALCULATIONS OF PHOTOGRAPHIC STAR POSITIONS	 83]

	The necessary multiplications, additions, and subtractions
could be made by successive operations on the multiplier,
but for such a large number of cards it is much cheaper to
perform part of the work on the sorter and gang punch. By
Lhis method about 98 per cent of the reductions in α, and 80
per cent of the reductions in δ require only one run of the
multiplier. The secular terms are computed with the aid
of the sorter, and punched with the gang punch as true figures
or as complements. The work on the multiplier is always 0£
the form

			A ± B ⅹ C ± D,

where A is the uncorrected coordinate, B ⅹ C is the precession
and D is the secular term which has been applied with the
gang punch. D also includes such additional digits as are
necessary to convert excess minutes and seconds in the
sexagesimal system. These factors must be determined before
the operation on the multiplier. For most of the stars
this is easily done with the aid of the sorter, and for the
comparatively small remainder, a trial run on the multiplier
is necessary.
	In the reduction of the right ascension the precession is
always positive. Hence if positive secular terms are punched
as true figures and negative ones as complements the operation
on the multiplier will always be the standard one

			A + (BⅹC) + D.

	The secular variations range from - 0.0250 to + 0.0250,
which means that the cards must be divided into 500 groups
of about 60 cards each. This requires three sortings. The
gang punching is most easily done by takmg every eighth
group in order that 12.5 s_α will change by one unit.
	Before starting the multiplier the cards must be divided
into groups according to whether 60, 120, or 180 seconds are
to be converted into minutes, and O or 60 minutes are to be
converted into hours. Let us consider first the question of
the seconds. The secular term is never greater than about
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0^{S}.3, and hence ran affect only those stars which are near the
critical values. In making the separation to the nearesL 1^{S}.0,
or for about 98 per cent of the stars, we can ignore this
quantity. In order to save time at the sorter we make this
separation in two steps, first making the separation to the nearest
10^S and then to 1^{S}. In each case we first divide the cards
according to the number of seconds in ao and then for each
group find the critical valuc of p_α.

[Graph:
		p_α ---> [Indicates X-axis is increasing p_α to the right]

α_0 [Indicates Y-axis is increasing α_0 as you go down]
|
|
V

[X-axis along the top 0^S   ^S.20  ^S.40    ^S.60     ^S.80     1^{S}.00  1^{S}.20 ]

[Y-axis on the left 10^S  20^S  30^S    40^S   50^S   60^S ]

[A line from the bottom left 60^S Y-axis to the top right 1^{S}.20 X-axis ]

[Text in top left corner "Convert 0^S " ]

[Text in bottom right "Convert 60^5" ]

"FIG. 14. Computation of Precession in Right Ascension"]

	The sorting and gang punching is simplified by the use of
the diagram shown in figure 14. We could plot each star
on this diagram according to the seconds of α_0 and the value
of P_{α}. Those lying below and to the right of the shaded
diagonal must convert 60^S into 1^{m}, while those lying to the
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left and above must not. For those within the shaded area
iL cannot be decided without considering the value of
12.5 s_{α}. If they are separated according to the sign of s_{α}
half of these can be eliminated. There are also two tiny
regions of uncertainty in the corners due to S_{α}, It is obvious
that for the majority of the stars the decision can be made by
sorting only in the tens of seconds in α_0 and the tenths of
seconds in p_{α}. A large portion of the remainder requires
sorting on the seconds of α_0 and the hundredths of p_{α}. The
process is continued to the point where it is cheaper to make
a trial run on the multiplier than to continue the sorting
process with small groups.
	The trial run on the multiplier is necessary only for a small
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percentage of the cards. ln order to avoid punching these
few cards and thereby making them different from the others,
the operator simply watches the numbers as they appear in
the summary products counter and separates the cards as they
come out of the machine.
	As soon as the number of seconds to be converted is determined
we can sort the cards on the minutes of α_0 and
determine which ones are to have a carry-over into the hours.
	Since we now know in advance where all the sexagesimal
conversions are to be made, we can incorporate the necessary
quantities into D of the machine operation A + (B ⅹ C) + D.
The punched result is thus the right ascension corrected for
p_α and s_α and expressed in the customary system of hours,
minutes, seconds, and decimals.
	For the precession in declination the problem is complicated
by the fact that δ_0 and p_ δ, may be positive or negative
and by the fact that s_ δ, is larger than s_. The cards are therefore
sorted first according to the signs, and separate diagrams
are used for the different cases. The diagram in figure 15 is
the one used for the case  δ_0 +, p_ δ +, s_ δ -, or for the case
 δ_0 - , p_ δ - , s_ δ +. In this diagram there are the added diagonal
lines indicating values of 12.Ss_ δ = 2", 4", 6", 8". In
this case also the test run on the multiplier is used for the
doubtful cards.

4. THE COMPUTATION OF THE RECTANGULAR COORDINATES

	This requires the evaluation of formulae (1) for each star
and for two overlapping plates, with an accuracy of seven
significant figures. This calculation illustrates the extensive
use of special tables in the form of files of punched cards
(see ch. IV). The special tables used are : (1) a table of
tangents with the tabular interval 10''. This table was made
from the values in Brandenburg's table. (2) A special table
of tangents and secants for angles less than about 6° and with
a tabular interval of O^{S}.1. This table was constructed for the
purpose with the aid of Brandenburg's table of cosines.
(3) A table of reciprocals multiplied by K, for all values of
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the denominator in formula (1). This table was constructed
for the purpose.
	The features of the machine work consist mostly of avoiding
multiplication where possible, and of going from one card
to another at the proper time.
	The first card which was punched from the A.G. Catalogue
is used also for the computation of the precession. There is
still some space left on card (1) and this is used for tan  δ' and
the tabular difference. There is not room for the interpolation
 on this card, but having found the tangent and difference
it will not be necessary to carry all of  δ' itself to the next card.
The application of the tangent and difference is of course
accomplished automatically by the method of chapter IV,
section 3. The cards are tabulated for a check.
	A second card is now punched with the high-speed reproducer,
transferring from the first card the star number and
the minutes and seconds of the right ascension. At the same
time the gang punch applies the complement of the right
ascension of the plate center, with an additional factor of
3600 for those stars which require a carry over of one hour.
It should be remarked that the right ascensions of the plate
centers have been rounded to the nearest whole second so the
decimals will not have to be subtracted.
	One run on the multiplier of card (2) now gives as a crossfooting
operation ∆α_1 and ∆α_2 in seconds of time. Checking
is done on the tabulator.
	The values of the tangent and of the secant of∆α_1 and of
∆α_2 together with the differences of the secant are now applied
and checked in the usual manner. The differences of the
tangent are not necessary, as they may be inferred from those
of the secant.
	The card files of tan ∆α, sec ∆α have been made with a
small interval of tabulation to facilitate interpolation. This
was done because the range of argument is small and there are
four interpolations for each star. The tabular interval is 0^{S}.1,
so in the seventh place the difference in the tangent is about
seventy-five and varies very little, while the difference in the
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secant is never more than eight but varies considerably. The
interpolation of both is done by one crossfooting operation.
A set of master cards is sorted into the detail cards according
to the value of the difference of the secant and according to
the hundredth and thousandth of a second of ∆α. In addition
to the quanlities used for sorting, these cards contain the
number to be applied to the tangent and secant for the
interpolation. The operation on the multiplier is done as rate
card crossfooting. Thus with two sortings, one crossfooting
operation, and one check run on the tabulator we are able to
perform and check two interpolations. The tangents were
made by subtabulation from Brandenburg; and to avoid
accumulation of error, were taken to eight places, and are
interpolated with eight. In order to save space on the card
the first three digits of the secant which are constant (1.00)
are not punched. The above is done for ∆α_1 and then for ∆α_{2}.
	The third card is now required. From card (2) are transferred
the star number, sec ∆α_1, and sec ∆α_1; and with the
gang punch is added the plate number and the sequence of
digits 1 2 3 4 5 6 7 8 9. From the first card is transferred
tan  δ' and its difference.
	The tangent of  δ' is interpolated on the multiplier as an
operation of the form A + B ⅹ C. Before this interpolation,
the cards are sorted according to the difference so as to form
constant multiplier groups which can be checked on the
tabulator.
	The product tan  δ' ·sec∆α  is now punched and checked on
the multiplier. The first digit, "1," of sec ∆α which was not
punched on the cards at the time of the interpolation is now
obtained from the sequence of digits mentioned above.
	The denominator of formula (1) is computed on the multiplier
as an operation A ± B ⅹ C, and checked on the tabulator.
Sin  δ is taken from a rate card and cos  δ, is made up
from the sequence of digits. This is done for ∆α_1 and then
for ∆α_{2}.
	Card (4) is now prepared. From card (3) comes the star
number, the denominator, and the quantity tan  δ' sec ∆α_{1}.
Card (5) contains the same quantities for ∆α_{2}.
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	The numerator of the Y equation is next to be computed.
The value of sin  δ is gang punched on the card and the cosine
comes from a rate card. The numerator is more troublesome
than the denominator for zones near the equator since cos  δ
is larger than sin  δ. We therefore have to differentiate between
cases where tan  δ ≷ tan  δ' sec ∆α. Checking is done
with the tabulator.
	From the special table we find K times the reciprocal of
the denominator and interpolate the result. Checking is done
by progressive digiting.
	We now reproduce from card (2) on to cards (4) and (5)
the values of tan ∆α_{1}, tan ∆α_{2}, and we are ready to form X
and Y by straight multiplication. These are eight-digit
multiplications and require eight runs for each star, including
the checking.
	Finally the results are listed with the tabulator.

[PAGE BREAK]

			CHAPTER X

		   STELLAR PHOTOMETRY

		    1. INTRODUCTION

	The Ruthrrfurd Observatory of Columbia University is
engaged in the determination of the relative brightness
of stars by means of photography. The present program
which should be completed in a few yea rs includes about
150,000 stars. While this work was started before the application
of the punched card, it is very doubtful if the program
in its extended form could have been under taken had it not
been possible to perform the necessary calculations by the
punched card method. The process is now so automatic
that the only personal attention given the individual star is
the actual reading of the photometer.
	In the photographic method of determining the brightness
or magnitude of stars two exposures are taken on each
plate: one of a standard region containing a selected sequence
of "comparison stars" and the other containing the stars
whose magnitudes are to be determined, the so-called "program
stars." The plate is then placed in a photometer which
measures the amount of light obscured by each of the star
images, the amount of obscuration being indicated by the
deflection of a galvanometer. The readings of the galvanometer
for the comparison stars of known magnitudes give a
calibration curve for converting the galvanometer readings
of the program stars into magnitudes. The photographs are
taken with a grating over the lens of the camera so t hat each
star gives a central image and two diffraction images, and
each star is photographed on two overlapping plates. For
each star on each plate are recorded the reading of the galvanometer
on the general background in the neighborhood,
the reading on the central image, and the mean of the readings
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of the Lwo diffraction images. For the brightest stars the
diffraction images only are used, and for the faintest ones only
the central image.
	The readings on the star images must be corrected for the
transparency of the background. The reduced reading for
the star image (central or diffraction) is given by subtracting
the reading on the star from the reading on the background
and multiplying the difference by the ratio of the standard
background reading to the observed. In the program here
described the galvanometer readings range from 0 to 200,
and 200 is the standard background reading. Hence if S is
Lhc reading on the star and G the reading on the background
near the star the corrected reading is (G - S) 200 .
					___
					 G

	The corrected readings for the comparison stars may be
plotted against the corresponding magnitudes to give a calibration
curve for the plate. From this curve may then be
read the magnitudes corresponding to the corrected reading
for each program star. However, a study of a large number
of these calibration curves has shown that within the experimental
error they may all be reduced to a standard curve by
a linear reduction. Thus if m is the magnitude read from the
correct curve for a given plate and m' is the magnitude read
from the standard curve, then

[equation		m =a + bm',	]

where a and b are constants for the plate in question. This
fact makes it convenient to substitute a numerical process
for the plotting and reading of the graphs. The standard
graph is made into a table from which m' is read. The
forty or so comparison stars are used to determine the constants
a, b and their probable errors. The magnitudes of the
program stars are taken from the table and reduced with the
a, b obtained from the comparison stars.
	In addition to the above reduction of the observations,
there is the problem of identifying the stars on the plate for
measurement. This identification is made by means of rec-
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tangular coordinates which must be computed from the known
spherical coordinates.
	When the magnitudes have been determined for a large
number of stars they are examined for systematic errors de-
pending upon the magnitude, the position on the plate, etc.

		2. THE USE OF THE PUNCHED CARDS

	The first step is naturally the punching of the cards from
catalogues such as the B.D. The data for each star are
punched on a card and verified. Each card contains the
number of the zone in question, the number of the star in the
zone, the right ascension and declination, and the visual
magnitude. With the gang punch is now added the number
of the plate on which the star will appear.
	The computation of the rectangular coordinates on the
plate is the same problem as that of chapter IX except that
for the present purpose very little accuracy is required. For
this low accuracy the series expansions in terms of ∆α,  ∆δ are
used and terms beyond the second order are neglected. To
simplify the computations, the centers of the plates are so
chosen that for a gffen zone they all have the same declination,
and the difference in right ascension between centers is
a constant number of whole minutes. Further, the scales
on the photometer in x and y were ruled in minutes of arc.
The formulae required are:

		X = k_{1}∆α^s - k_{2}∆α^{s}∆δ',
		Y = ∆δ' + k_3(∆α^s)^{2},

k_{1}, k_{2}, k_3 arc constants for the zone,

∆α^S and ∆δ' are the differences in α and δ between the star and
	the center of the plate, the first expressed in seconds of
	time and the latter in minutes of arc.

	The computation of X_{1}, Y_1 and X_{2}, Y_{2}, the coordinates of
a star on two different plates, is accomplished by a small
amount of sorting and gang punching, the pulling of one card
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from a file, and one summary punch operation. The first
operation wilh the sorter and gang punch gives ∆δ'. The
second gives the minutes of ∆α. The third gives the small
term k_{2}∆α^{s}∆δ' for one of the plate centers, and then for the
other. A card is now drawn from the file with argument
∆α^{m,s} and giving k_{1}∆α^s and k_3(∆α^{s})^{2}. This card is still taken
from the file by hand because there arc very few of them and
because the question of duplicates is troublesome when the
sorter is used instead of the collator. The tabulation of the
original card and that pulled from the file gives at one run
the coordinates of the star as it will appear on two plates.
The summary punch gives a new card containing the zone, the
plate number, the star number, the visual magnitude, X_{1}, Y_{1},
X_{2}, Y_{2}.
	The cards just obtained are to be used at the photometer.
To facilitate the reading of the cards they are passed through
the interpreter which prints at the top of the card the numbers
which are punched on the card. Thus, in type, at the top of
the cards appear the number, the visual magnitude, and the
four coordinates.
	The observer uses the printed numbers to locate the star,
and then places the card in a card punch to record the readings
of the galvanometer directly on the cards. By the time the
galvanometer has come to rest the card is in position for
punching, and the scale reading for background, central image,
and mean of the diffraction images are recorded. Each card
is used at the photometer in connection with two plates.
	For the machine application of the background correction
we rearrange the formula:

					  200
					  ___
		Corrected reading = 200 -  G ·S.

The instrument is adjusted so that 200 varies little. The varies little. The
				   ___
				    G
cards are sorted on G, and 200 is used as a rate card on the
			   ___
			    G
multiplier to compute 200·S. Since G differs very little from
		      ___
		       G
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200 we can form this quantity for both the central and the
diffraction images with one run of the machine. This is done
as an operation A ± B ⅹ C where A is the reading on the
central image and the reading on the diffraction image, B is
200 - 1, and C is the two largest digits of each of the two
___
 G
parts of A. As an example take the case where the reading
on the background is 195, the reading on the central image is
136, and the mean of the diffraction images is 183. Then

			A = 13600183
			B = 	    .026
			C = 13000018

The answer is 139.38, 187.68 which are rounded to 139.4 and
187.7. These quantities may be in error by as much as 0.26
because of the omission of the end digits of C, and by greater
amounts for larger values of 200 - G.
	The table giving the provisional magnitude m' is made in
the form of a file of punched cards which may be used in the
usual manner. To avoid the subtraction of 200·S from 200,
				          ___
				           G
the 200 is absorbed in the argument of the table. The opera-
tion of "looking up" the number in the table requires about
five hours for 20,000 entries.
	The values of m' for the comparison stars must now be
used in conjunction with the corresponding values of m for
the determination of the plate constants a and b by the
method of least squares. There are n (about forty) observation
equations of the form

		a + m'b - m = 0.

The normal equations are

		n·a + [m']·b - [m] = O,
		[m']·a + [m'^2}·b - [m m'] = o,
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and the solution is given by
		   [m][m'^2] - [m'][m m'J
		   ______________________
		a = n[m'^2] - [ m']^2

		b =n[m m'] - [m][m']
		   _________________
		    n[m']^2 - [m']^2


The formation of the coefficients of the normal equations is
clone very simply on the machines. The card for each comparison
star contains m and m', so m m' and m'2 can be
formed with two runs of the multiplier. (Two runs are required
because of the wide variation in a and b.) The tabulator
then gives with one run of the cards n, [m'], [m], [m'^2],
[m m']. If relatively few plates are being reduced at one
time the normal equations can be solved by hand from the
quantities given on the tabulator record, and summary cards
are not required. If say five hundred plates are reduced
at once it might be desirable to eliminate this handwork.
The summary cards are run through the multiplier to give
[m][m'^2], [m'][m m'J, n[m'^2], [m']^2, n[m m'], [m][m'].
The crossfooting device then gives the denominator and the
two numerators. A small table of reciprocals is required to
avoid the division. One run on the multiplier for the interpolation
of the reciprocal and two more for the multiplication
are then required to give a and b.
	Having found a and b for each plate, we can evaluate m
for the program stars by a run on the multiplier of the form
			A + B ⅹ C,
where A = a, B = b, C = m'. This quantity should be
formed at the same time for the comparison stars in order to
determine the weight of the solution for a and b. This ':eight
is determined from the sum of the squares of the reS1<luals
where the residual v = m - a - b m'.
	The values of m may now be used to form the mean of the
observations for each star and to form the difference for each
star of the values obtained from the two plates. These are
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crossfooting operations except for the multiplication by 1
						     _
						     n
in forming the mean. n is here 2, 3, or 4 depending upon
whether or not central and diffraction images were measured
on both plates. The ,values of the difference for the two
plates are used to discuss the accuracy of the results. The
values of the mean of this quantity for a given plate is subtracted
from the individual values, and the outstanding discrepancies
are reobserved. These outstanding differences
sometimes show that the wrong star was observed on one of
the plates or that the star is variable. Finally the differences
are examined for systematic effects depending upon t he magnitude,
the X and Y coordinates on the plate, and the color
(difference between photographic and visual magnitude). All
these quantities are on the card and the necessary computations
are done most effectively on the tabulator and counting
sorter.
	Finally the definitive magnitudes are listed with the tabulator
for publication.

[PAGE BREAK]

			   CHAPTER XI

		NUMERICAL LUNAR THEORY

			1. INTRODUCTION

	In this ch pter the solut ion of a set of differential equations
is determined in the form of harmonic series with
numerical coefficients and with arguments which are linear
functions of ·the independent variable. An approximate solution
is substituted in the equations, and from the residuals
corrections are obtained for the coefficients of the assumed
solution. The project here described is the most extensive
which has been done, and it was feasible because of the availability
of the punched card method.
	It is well known that the solution of the main problem of
the lunar theory may be expressed with great accuracy as
harmonic series giving the coordinates in terms of the time.
The argument of each term is a linear function of the time
with four fundamental frequencies. The coefficients in these
series and two of the frequencies, the motions of the perigee
and node, are determined from the differential equations of
motion. In the well-known theories of Delaunay and Brown
these coefficients and motions were determined by algebraic
expansions in powers of the arbitrary constants, numerical
values of the constants being inserted only at the end. These
developments are built from the ground up, and no use is
made of previous developments except as a general guide.
Hansen's theory, though numerical. followed the general plan
of the algebraic theories and was independent of other solutions.
	The basic principle involved in the present development
is essentially different. It consists of substituting in the
differential equations of motion the numerical values of the
best available previous solution. The residuals thus obtained
are then used with the equations of variation to determine
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small corrections to each coefficient and to the motions of the
perigee and node. This method gives at the same time a
precise verification of the previous work, and increased accuracy
where necessary. It utilizes the previous theory, but the
final results are independent of it. The coefficients involved
are ordered solely according to numerical value and not
according to the algebraic order of the parameters involved.
The only assumption is that the square of the errors may be
neglected.
	The numerical method by itself would leave much to be
desired, but when used in connection with a good literal
theory it gives results which could be obtained by the literal
method alone only at the cost of tremendous effort. Thus
in the present case the general accuracy of the theory of
E. W. Brown was ample for observational purposes, but a
few coefficients, particularly those giving the motions of the
perigee and node, had not been obtained with sufficient accuracy.
To obtain the additional accuracy by an extension of
the original method would require years of highly specia lized
work. Then too there was the possibility of error in the
original work which had required twenty years to complete.
By the numerical method it has been possible in about two
years to verify completely Professor Brown's results and to
obtain the required additional accuracy.
	The general method here employed was suggested by Sir
George Airy about half a century ago. He spent fourteen
years and considerable sums of money on the project but
was unable to complete it successfully.* The numerical work
involved is so extensive that with the methods then available
numerical errors were almost unavoidable. In addition, as
pointed out by Radau,† he had made one or two mistakes in
theory.

			2. THE PROBLEM

	The problem investigated here, known as the main problem
of the lunar theory, is the determination of the motion of

* Monthly Notices of the Royal Astronomical Society, 1888, vol. 49, p. 2.

† Bulletin Astronomique, 1887, vol. 4, p. 274
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the moon under the gravitational action of the earth and the
sun, with the assumptions that the three bodies are spherical
and that the sun moves in a fixed ellipse about the center of
gravity of the earth-moon system.
	The position of the moon is referred to rectangular axes
through the earth's center with x and y in the plane of the
sun's orbit, and rotating about the z axis with the mean angular
velocity of the moon.
	The equations of motion may be written in the homogeneous
form.

(1) xDy - yDx + (x^2 + y^2) = D^-1 (x δΩ -  y δΩ) + C_{1},
			     ____     ___
			       δ       δ

(2) ½D^2(r^2) + ½(x^2 + y^2) - ½[(Dx)^2 + (Dy)^2 + (Dz)^2]

= Ω + [r δΩ] ] +D^{-1} [ 2 ( x δΩ - y δΩ) -1 δΩ] + C_{2},
	__		 __     __ _ __
	δr		 δy	δx n δt


(3) xD^{2}z + z [- D2x + x + 2Dy + δΩ - x δΩ] = 0
			     __ _ __
			     δx z δz

where x = 1 + x_1 and x_1, y, z are the deviations from circular
motion in the plane of the orbit.

		Ω = ∑_{j}∑_{h}k_{jh}ω_{jh}

where k_jh  contains successive powers of a small constant and
each ω_jh is the product of a homogeneous functions of x, y, z
of order j by a harmonic series with argument t and  with
known numerical coefficients.
	The adopted solution is of the form
			x = ∑a_i cos α_{i},
			y = ∑b_i sin α_{i},
			z = ∑c_i sin α_{i},
		  α_{i} = (jp + kq + lr + ms)t
			j, k, l, m = 0, ± 1, ±2···
			p, q, r, s = given numbers.

The derivative series are formed by multiplying the coefficient
of each term by the coefficient of t in α_{i}, and changing the
signs where the harmonics are cosines. For each term, the
numencal value of D is the coefficient of t in α_{i}. To substitute
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this solution in the equations ( l ), (2), (3) it is necessary to form
all the product series of the second order for the parts independent
of Ω, and for the terms of second and higher order
involved in Ω. Since the individual parts of Ω contain small
factors, these latter are not required to the same accuracy as
the former. By far the largest part of the work consists of
forming the product series for the second order terms. The
coefficients in the assumed solution were taken with ten
significant figures and there were about 500 terms for each
coordinate.
	For comparison with observation the rectangular coordinates
must be connected into polars. The equations for
this are:

	tan(V-nt) = y,
		    _
		    x
	    sin φ = s
		    _,
		    r

		    1
		    _
		π = r.

Here again the transformation may be checked by substituting
the numerical coefficients for the series representing x, y, z;
V ,  φ, π. We find from the above equations

		(x^2 + y^2)Dv = xDy - yDx,
		φ = zπ + 1/6(zπ)^3 + 3/40(zπ)^{5},
			r^{2}π^2 = 1.

		3. THE METHOD WITH THE PUNCHED CARDS

	The method described in chapter VII for large series was
employed. here. The multiplications were done as trigonometic
series in some cases and as exponential in others, so that
all the cases of that chapter occurred. Conjugate imaginaries
u = x + iy and s = x - iy were employed to facilitate the
evaluation of Ω and its derivatives. There were about eight
hundred terms in u, and in forming u^2 and u s about twenty
thousand distinct products greater than 1 X 10^{-11} had to be
considered for the sum and for the diffcerence of the arguments.

[PAGE BREAK]

			CHAPTER XII

		THE COMPUTATION OF PLANETARY
			PERTURBATIONS

			1. THE PROBLEM

	The computation of planetary perturbations furnishes
excellent examples of the methods of chapters VI, VII, and
VIII. The methods of chapters VI or VII or both would be
used in the so-called method of general perturbations while
the application of chapter VITI is known as the method of
special perturbations.
	The differential equations of planetary motion are

		d^{2}x + k^{2}x = δΩ = k^{2}∑m_i x_i - x 
		______   ______   __		 _________
		dt^{2}     r^3    δx		 ρ_{i}^{3}

with similar equations for y and z.
		   m_i
		   ___,
	Ω = k^{2}∑ ρ_i

	ρ_{i}^{2} = (x_i - x)^2 + (y_i - y)^2 + (z_i - z)^{2},

	m_i = ratio of the mass of the disturbing planet to that of
		the sun.

	The integration of these equations is always based on the
fact that m_i is small, never greater than 1 ⅹ 10^{-3}. When
m_i = 0 the solution is an ellipse with the well known law of
areas. The elliptic solution may be substituted in the right
hand side of the equation and a second approx1mat1on to the
solution made which gives the so-called first order perturbations.
The substitution of this then gives the second
order, etc.
	The solution of the complete equations may be given in
the form of a table of numerical values of x. y and z for equidistant
intervals of time through a certain range, or in the
form of harmonic senes with terms whose arguments are linear
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functions of the time. The expression of the values of the
coordinates as a table is the nwthod of special perturbations,
while that employing the harmonic series constitutes the
method of general perturbations. The advantage of the
method of general perturbations is that it gives explicit expressions
for the coordinates in terms of the time. They may
be used for computing the position of the planet at wide!y
separated times without computing all the intermediate
positions. The determination of these expressions with high
accuracy is howeYer a very laborious process. The tab!e of
numerical values ts useful only for the range of the table, but
it has the advantage that for limited intervals of time the
table may be constructed with great accuracy and comparatively 
little work.

	2. THE COMPUTATION OF SPECIAL PERTURBATIONS*

	This problem furnishes an ideal example of the method of
chapter VIII. There are nine planets with m_i large enough
to be cons1dered and about 1,500 planets and about 500 comets
to which the process could be applied without changing the
machme setup. In other words the necessary initial cost of
the application of the punched card method can be distributed
over a great many cases.
	The card files required are: (1) a table of quantities to
the minus three halves power and (2) tables of the coordinates
of the nine perturbing planets. The first of these tables was
made from Comne's table and is accurate with linear interpolation
to eight places. The coordinates of the perturbing
planets were prepared from Comrie's table and other sources.
The coordinates are tabulated at ten-day intervals with five
decimals in x, Y, z. In order to use the center of the sun as
origin instead of the center of gravity of the solar system we
have on the same card the attractions of the perturbing
planets on the sun.
	The quantity k^{2}∑m_i  x_i - x  is called the perturbation in x,
			       _________
			       ρ_{i}^{3}


	*The Computation or Special Perturbations by the Punched Card Method.
W. J. Eckert. Astronomical Journal No. 1034, 1935.

[HEADER:	THE COMPUTATION OF PLANETARY PERTURUATIONS 	103 ]

and there are similar expressions in y and z. These perturbations
are computed before the main process of integration is
started. Since m_i are small quantities we can use approximate
values of x for this purpose. Obviously if "elliptic"
values are used for x, y, z the integration will give coordinates
which include the first-order perturbations. If first-order
perturbations are included in the approximate coordinates,
the results will contain second-order perturbations.
	The computation of the perturbations is a straightforward
process on the machines. It is most efficiently done as mass
production. The more disturbed and disturbing planets and
the more dates included the more efficient it becomes. The
cards containing the approximate values of x, y, z are sorted
together with the cards containing the coordinates x_{i}, y_{i}, z_{i},
and the sorting is done according to the date. The tabulator
and summary punch then give a new card for each date containing
(x_i - x), (y_i - y), (z_i - z). Repetition of this process
gives successively the corresponding quantities for i = 2,
3,···, 9, or as many of them as are required. The card must
carry in addition to the date and the differences in the coordinates,
numbers to represent the values of i for the disturbing
planet and of j for the disturbed planet.
	When all the differences of the coordinates have been
formed they are passed successively through the multiplier
to give (x_i - x)^2 , (y_i - y)^2, (z_i - z)^2 and ρ_{i}^{2}. At this point
it is necessary to make sure that no accuracy is lost due to
close approaches of the disturbing and disturbed planets,
since the force increases as the distance decreases. A simple
method for doing this is given in the last reference.
	The values of     1    corresponding to the values of ρ_{i}^{2} are
		      _________
		      ρ_{i}^{3}
taken from the table in the usual manner with the sorter and
high-speed reproducer. Interpolation is done with the multiplier.
The multiplication by m_i is done on the multiplier.
In addition there is a factor (10^h)^{-3/2} which must be applied
because our table of quantities to the minus three halves
power was constructed with arguments from 2 to 20 whereas
ρ_{i}^{2} is frequently outside this range.
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	Final multiplilication by (x_i - x), (y_i - y), (z_i - z) gives
the· perturbations due to each planet.
	Since the computation of perturbations is complicated,
they are computed for fewer dates when possibk and the missing
ones found by interpolation. The subtabulation is usually
to halves and quarters and is done automatically by the
Lagrangian method described in chapter V section 3.
	Finally all the perturbations for a given disturbed planet
for a given date are collected and added with the tabulator.
The summary card is made in red lo facilitate its identification
during the process of integration. These perturbations may
now be considered as a function of the time for the process
of integration.
	The equations for the integration proper thus reduce to
			X = - k^2  x + P(t)
				  ___
				  r^3

the solution of which is given by

	x = ''X + l /12X - 1/240 X'' + 1/1951^Ⅳ + ···
	  = ''X + the "reduction,"
with similar expressions for y and z.
	The following table shows the various cards used in the
computation of the perturbations and in the integration.

Card				Information
1   0	approximate coordinates, x, y, z
x_1 0	same with signs changed
1   i	x_i - X, y_i - y, z_i - z; ρ_{i}^{2}
2   i	computation of k^{2}m_{i}/ρ_{i}^{3}
3   i	multiplication of x_i - x, y_i - y, z_i - z, by k^{2}m_{i}/ρ_{i}^{3}
3   y	sum of perturbations
4   i	ith difference of the second summation of the total
	   attraction

5   i	extrapolated value of (4i)
6   0	integrated values of x, y, z and the values of the solar
	   attraction computed from them.
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Additional cards used in the interpolation of the perturbations
are not shown. On each card appears the date in the form of
the Julian day, the identificalion of the particular card in the
general scheme, and the quantities required. Six columns
are used for the date, and two for the identification. Cards
1, 2, and 3 are used in the computation of the perturbations
while 4, 5, and 6 are used in the integration.
	The integration is best described by means of an example.
Figure 16 is a portion of the tabulator record and shows all
the totals printed during three steps of integrat1on with
seventh differences. Each line is printed at one time. The
first group of ten columns comes from the first counter of the
tabulator and contains the date and the identification of the
operation; the second, third, and fourth groups come from
the second, third, and fourth counters and contain information
pertaining to x, y, z respectively. In the first group of columns,
digits 2 to 7 indicate the Julian day while columns 1 and 8
indicate the quantity contained in a given line. Thus from
the seventh line (46) in the figure we find that for J. D.
[2]427980, X^Ⅳ = + 1344, Y^Ⅳ = + 1585, and Z^Ⅳ = + 444
and find also from the line (40) that ''Y = - 1.846294849.
	Summary cards are punched for each line except those
identified as (49) and the second (40). The first two lines (49)
are steps in the computation of the last (49) which contains
the true seventh differences of the attraction.
	There is in addition a card with the same identification
as (60) except for a "y" punched in columns 77-80. This
card is used in computing the solar attraction.
	In the machine process we consider a step completed when
we have just finished computing the differences and summa-
tions. Thus on the tabulator record have just been printed
lines (49), (48), ···, (40), and the corresponding summary
cards have been punched. The new step starts with the
difference and summation cards and proceeds to the estima-
tion of the quantities needed in the reduction. The value of
(42) must be estimated one step in advance with considerable 
accuracy. This is easily done on the tabulator by simply
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[table/ photo image of tabulator output

4|427980|900		237232-		8S96969-	862302-
4|427980|900		237232-		8596969-	862302-
4|427980|900		     7-		     24-	    54
4|427980|900		     7-		     24-	    54
_______________________________________________________________
4|427980|800		    95 -	    119-	     4-
4|427980|700		   346		    379-	    52
4|427980|600		  1344		   1585		   444
4|427980|500		  9098-		   6165		  1356-
4|427980|400		 27471-		  67496-	 12239-
4|427980|300		648272		 150201-	124235
4|427980|200		23732S		8596945		862356
4|427980|100	     136441552-		5211857-      29629728-
4|427980|000	     187384632-	     1846294849-     214362843-
4|427980|000	     187384632-	     1846294849-     214362843-
5|427970|200		850523-		8386500-	973488-
5|427970|200		850523-		8386500-	973488-
5|427970|200		850523-		8386500-	973488-

6|427980|000	     187313551-	     1845596017-     214281644-
6|427980|000	     187313551-	     1845596017-     214281644-
6|427980|000	     187313551-	     1845596017-     214281644-
6|427980|000	     187313551-	     1845596017-     214281644-

4|427970|900		850523-		8386500-	973488-
4|427970|900		850523-		8386500-	973488-
4|427970|900		    20-		     82		    54-
4|427970|900		    20-		     82		    54-
_______________________________________________________________
4|427970|800		   115-		     37-	    58-
4|427970|700		   231		    416-	     6-
4|427970|600		  1575		   1169		   438
4|427970|500		  7523-		   7334		   918-
4|427970|400		 34994-		  60162-	 13157-
4|427970|300		613278		 210363-	111078
4|427970|200		850503		8386582		973434
4|427970|100	     135591049-		3174725	      28656294-
4|427970|000	     322975681-	     1843120124-     243019137-
4|427970|000	     322975681-	     1843120124-     243019137-
5|427960|200	       1422955-		8124107-       1070811-
5|427960|200	       1422955-		8124107-       1070811-
5|427960|200	       1422955-		8124107-       1070811-

6|427970|000	     322856863-	     1842443175-     242929821-
6|427970|000	     322856863-	     1842443175-     242929821-
6|427970|000	     322856863-	     1842443175-     242929821-
6|427970|000	     322856863-	     1842443175-     242929821-

4|427960|900	       1422955-		8124107-       1070811-
4|427960|900	       1422955-		8124107-       1070811-
4|427960|900		    17-		     31		    64
4|427960|900		    17-		     31		    64
_______________________________________________________________
4|427960|800		   132-		      6-	     6
4|427960|700		    99		    422-
4|427960|600		  1674		    747		   438
4|427960|500		  5849-		   8081		   480-
4|427960|400		 40843-		  52081-	 13637-
4|427960|300		572435		 262444-	 97441
4|427960|200	       1422938		8124138	       1070875
4|427960|100	     134168111-	       11298863	      27585419-
4|427960|000	     457143792-	     1831821261-     270604556-
4|427960|000	     457143792-	     1831821261-     270604556-
5|427960|200	       1950322-		7818013-       1154643-
5|427960|200	       1950322-		7818013-       1154643-
5|427960|200	       1950322-		7818013-       1154643-

6|427960|000	     456981061-	     1831169887-     270508299-
6|427960|000	     456981061-	     1831169887-     270508299-
6|427960|000	     456981061-	     1831169887-     270508299-
6|427960|000	     456981061-	     1831169887-     270508299-

		FIG. 16. Tabulator Record for Numerical Integration  ]
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adding the cards (42) , (43), · · ·, (48) just obtained. This
process obviously gives the desired result with an error equal
to the value of (49) which will appear when the step is completed.
The results will be (52) for the next date, and in
order to simplify the multiplication by 1/12 the results are
punched on three separate cards for x, y, z. The line (52) is
thus printed three times, once when the x card is punched,
once when the y and once when the z. In order to secure the
proper change in date and identification a colored master card
is tabulated with the group.
	The contribution to the reduction from the first term is
computed accurately with the multiplier, but a simpler process
is usually adopted for the second term. It has been found
that in most cases the quantities in (44) may be estimated by
inspection with sufficient accuracy. A small file of cards has
been made up with argument equal to the estimated values
of (44) and with the corresponding contribution to the reduction.
These cards are pulled from the file while the (42) term
is being multiplied by the machine. The third term (46) is
usually negligible.
	The new values of x, y, z are now evaluated from their
component parts by the tabulator. The totals (60) are
printed four times as three cards are required for x, y, z and
an additional card with the same identification is required for
		  k^2
		 _____
the computation of r^{3} .

	The cards containing x, y, z are passed through the multiplier
to give r^2 in the summary products counter. This
					  1
					 ___
number is used to draw a card from the file of r^3 . This file
card is copied onto the fourth (60) card by the reproducing
mechanism of the punch. Four of the eight digits of r^2 arc
used in selecting the card from the file and the remaining four
are punched by hand for use in the interpolation. The interpolation
which is linear is done automatically on the
machine. The interpolated value of   1   is then multiplied by
				___
				r^3
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the gravitational constant, and the result is ready for multiplication
hy x, y, s. These quantities, if added to the perturbations, 
would give the new ,values of (42) to he entered
in the difference table.
	For simplicity of the machine operation we first find the
highrst-order difference and build the others including the
two summations by progressive totals (See ch. V,sec. 4).
The highhest-order difference is ohtained as follows.

	  ∇_1^{i} = f_1 - [f_0 + ∇_{o}' + · · · + ∇_{0} ^{i-1}]

where f is the attraction. In this case we add the solar attraction
and perturbation and subtract the differences for the last
step. This leaves  ∇_1^{i} in the counters of the tabulator, and
by adding successively  ∇_{0} ^{i-1},  ∇_{0} ^{i-2}, · · ·,  ∇_{0} ^{-1},  ∇_{0} ^{-2} we obtain
the values of all the required differences and summations of
X, Y, Z. These quantities are printed on the record as
(49), (48), · · ·, (40). One step of the integration is thus
completed.
	The various operations in the above cycle are repeated
for each step in the integration and always in the same order.
All the necessary switching of the machine has been taken
care of in the calculation control switch as will be described
presently. The calculation control switch has twelve distinct
settings and these are run through in connection with
the various steps outlined above. For six of the settings the
switch moves automatically to the next, but for the others
the operator must press a conveniently located button.
	We now describe the activities of the operator during one
cycle.
	The calculation control switch is set on position 1. The
operator takes cards (4i) from summary punch, removes
(40), (41) and adds a yellow master card and two blank cards.
These are tabulated in the direct position. Three summary
cards (52) appear while the switch automatically moves to
positions 2, 3, 4.
	Operator takes a blank card, a salmon rate card, and
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three cards (52) from the summary punch, and passes them
through the multiplier. While multiplying ·is being done the
operztor reads fourth differences from tabulator record to
nearest 10,000 and writes them on a computing sheet. From
this sheet  he estimates these quantities two steps in advance
and selects cards from the file for this part of the reduction.
 (Part of this may be done as soon as the differences appear.)
He now  tabulates the following cards: (40) direct, three cards
 (52) from  multi lier reversed, the above cards from the file,
a blue master, and three blanks direct. The summary punch
gives four cards (60) while the switch moves automatically
to 5, 6, 7.
	Operator places first three (60) cards in multiplier. While
multiplier is in operation he lifts lever on summary punch to
prevent reset, and presses ekect lever. He substitutes in the
punch the fourth (60) car reversed for the blank. By this 
time r^2 appears in the summary products counter of the 
multiplier, and the operator pulls the proper card from the
file of 1  . This is placed in the duplicating rack and copied
       ___
       r^3

on (60) and the remaining digits of r ^2 are punched by hand. 
	The button is pressed to move the switch to 8 and the
card from the punch is passed through the multiplier direct.
	The switch is moved to  9, then blank, green rate card
direct, and (60) card of last operation reversed are passed 
through the multiplier.   At this stage the card from the file
is replaced and the lever on the punch released.
	The switch is moved to 10, and into multiplier is placed:  
blank, above (60) reversed, and the other three (60) direct.
While multiplier operates, tabulate (42), (43), · · ·, (48) and
pink master. 
	The switch is moved to 11 and the following tabulated:
red perturbation card and three cards (60) reversed from
multiplier.
	Inspect numbers in counters 2, 3, 4. These are the values
of the highest differences. If they are satisfactory, proceed.
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	Move switch to 12, tabulate all (4i) cards including (40).
This completes one step in the integration.
	We have not indicated all the details of picking up and
laying down the cards and putting them away. The efficiency
of the process requires that these be done without lost motion.
	The following is the punching on the various cards: The
ordinary detail cards in the hopper of the summary punch are
prepunched with a "0" in column 71 and an "x" in column 76.
The yellow master contains a "0" in column 10 and 5400012
in columns 51-57 for forward integration and 5399992 for
backward. The salmon rate card has 33333338 in columns
41-48 and an "x" in column 46. The pink master contains
"0" in column 10 and a "4" in column 51 and numbers in
columns 52 57 which depend upon the backward or forward
integration and upon the order of the highest difference. If
i is the order of the highest difference we have 4000i for
forward and 3999i for backward.
	The summary cards (4i) have four ten-digit fields: argument
and identification, functions of x, functions of y, func~
ions of z. The summary cards (52), (60) have two ten-digit
fields: argument and identification; function of x, or function
of y, or function of z. The card which refers to x has "x" in
punches in columns 9-10, that which refers to y has "x" in
column 9, and that which refers to z has no "x". The "x's"
in these two columns operate two class selectors on the tabul~
tor so that functions of x always go into counter (2), functions
of y into (3), and functions of z into (4). In these
operations with the summary punch there are multicontact
relays which punch the successive cards from the proper
counters and indicate the results with the proper combination
of "x's."
	The tabulator and multiplier have had minor modifications
to permit remote control of their operation through the
calculation control switch. These modifications become effective
only when the cables from the switch are plugged into
the sockets on the machines, so that when the switch is un-
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attached the machine operates in the normal manner. The
modifications described here are in addition to those in chap-
ter II.
	The special features of the tabulator are three in number:
a special circuit for clearing the counters by a major break in
control even though the dial switches are set for progressive
totals; a third class selector operated by the switch rather
than by x punches on the cards; a switch to control summary
unchmg of minor totals. The summary punch has an add1-
tional electromagnet to eject cards automatically at more 
than one column of the card.
	The multiplier has special provisions for the following
automatic changes:

	(1) from punch to nonpunch,
	(2) clearing the summary products counter,
	(3) from rate card to standard multiplication.
	(4) to special crossfooting.

	In addit ion to thcse changes in the machines themselves,
the calculation control switch has six multicontact relays
ough which the w1~:!::m~ary punch and multith:
plier are plugged. These relays give a very quick change in
wiring, and are controlled by the switch. Relay number l
selects the multiplicand from one of two fields on the card.
Number 2 selects the multiplier from one of two fields.
Numbers 3 and 4 cause the product punched on the card to
come from one of the three positions in the summary products
counter. Relays 5 and 6 go to the summary punch and
punch in the second field of the card numbers from the x, y,
or z counters of the tabulator. Additional relays would be 
convenient, but by careful planning of the cards it has been 
possible to carry out all the wiring arrangements.
	The upper half of figure 17 shows the operation cycle of the
cams of the calculation control switch. An ⅹ means the
contact is closed. The lower half shows the specific relays
and switches operated by each cam.
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[Table 

Cam No.			Switch Position

	1   2  3  4  5  6  7 8 9 10 11 12 | 13 14 15 16 17 18 19 20
1	ⅹⅹⅹⅹⅹⅹⅹ   ⅹ ⅹⅹ | 
2		   ⅹ	       |
3		    ⅹ	       |
4	ⅹⅹⅹⅹ	    ⅹ  ⅹ ⅹⅹ |
5	    ⅹ   ⅹ	       |
6	  ⅹⅹ  ⅹⅹ	       |
7	ⅹⅹⅹ		ⅹ     |
8	ⅹⅹⅹⅹⅹⅹⅹⅹⅹⅹ     |	Not Used
9	ⅹⅹⅹⅹⅹⅹⅹⅹⅹⅹⅹ   |
10	ⅹⅹⅹⅹⅹⅹⅹⅹⅹ       |
11	ⅹⅹⅹⅹⅹⅹⅹⅹⅹ       |
12	ⅹⅹⅹⅹⅹⅹ    ⅹⅹⅹⅹ |
13		     ⅹ	       |
14		    ⅹ	       |
15		       ⅹ      |
_______________________________|
16	
17	Not used
18	


Cam No.
	1-	multicontact relay no. 1 to select multiplicand
	2-	multicontact relay no. 2 to select multiplier
	3-	multicontact relay no. 3 to control punching of products
	4-	multicontact relay no. 4 to control punching of products
	5-	multicontact relay no. 5 to control summary punching
	6-	multicontact relay no. 6 to control summary punching
	7-	multicontact relay no. 1A in tabulator
	8-	multicontact relay no. 1B in tabulator
	9-	multicontact relay no. 2A in tabulator
	10--	counter clear on major control
	11-	summary punch on minor control
	12-	rate card multiplication
	13-	A - B ⅹ C
	14-	punch-check
	15-	"x" punch in column 44.



FIG. 17. Calculation Control Switch for Computat1ion of Planetary Motion ]


	3. THE COMPUTATION OF GENERAL PERTURBATIONS
	The method of general perturbations requires the expansion
of the disturbing function Ω = k^{2} m_i and its derivatives
				 ____
				 ρ_{i}
as harmonic series with a rguments which are linear functions .
of the independent variable There is a wide variety in the
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choice of variables and of methods of expansion, and the
choice of the particular one to be used is dictated by the
nature of the problem. From the computational standpoint
however they are very similar, and the methods outlined here
can be readily applied to any case in point.
	The expansions are made in powers of four parameters:
The ratio of the semimajor axes of the disturbed and disturbing
planet α, the eccentricities e, e', and the mutual inclination
η. The exact definitions of the constants vary with the
method, but we shall employ the above symbols for the
general case. The parameters e, e', η are usually small and
the series involving them converge very rapidly, but the
ratio α is frequently very large and gives rise to slow convergence.
For this reason the effects of α are considered
first, and then the parameters e, e', η are introduced by
multiplication of series or by harmonic analysis and synthesis.
Efficient methods for both have already been given and the
application to this problem is obvious.
	The fundamental expansions required at the outset are
series giving t e fu nction

(1 + α^2 - 2α cos s)^{-1/2, - 3/2, -5/2, -7/2}

as a Fourier series in S . These series are laborious to compute
but since they involve only one parameter they have been
tabulated.*
	The expansions are:

(1 + α^2 - 2α cos S)^{-S} = ½β_{S}^(o) + ∑_{i=1}^{∞} β_{S}^{(i)} cos iS

		      = (1 - α^2)^{-S}(½G_{S}^(0) + ∑_{i=1}^{∞}G_{S}^{(i)}α^{i}cos iS).

The quantities in the tables are log G_{S}^{(i)}. They are given to
as many as eight decimals and must in general be interpolated
with second differences. They are tabulated for i = 0, 1,

* Tables for the Development of the Disturbing Function, by Ernest W.-
Brown and Dirk Brouwer. Transactions of the Yale Observatory, 1932, vol. 6
pt. V.
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⋯, 11, and s = 1/2, 3/2, 5/2, and 7/2. It must be remem-
bered that there is a different value of α for each disturbing
planet. The required derivatives with respect to α are obtained
from the same table, for

α d (1 + α^2 - 2α cos S)^{-S}
 __
 dα

		         =   ∑_{i=0}^{∞}  Sα^i  (G_{s+1}^{(i)} - G_{S}^{(i)}) cos iS,
				         ________
				        (l _ α^2)^S

the factor ½ to be inserted for i = 0.
	We thus see that even with the tables available the work
of obtaining the necessary value of  β_{S}^{(i)} and their derivatives
is quite considerable. In order to perform this part of the
calculation by the machines the tables have been punched on
cards, each card containing a value of log G_{S}^(i) together with
the first and second differences. The selection of the proper
cards from the file is easy, for there are only a few values of  α
required and all the cards for a given value of  α should be
kept together. The interpolation is done on t he multiplier,
and the factors (1 -  α^{2})^{S},  α^{i}, s are applied with the tabulator
or multiplier. No provision has yet been made for finding
the antilogarithms by machinery. This is a simple matter
by hand with the ordinary tables, but it could easily be done
automatically if sufficient use could be found for such a card
catalogue.
	Once the series involving a have been found, the other
expansions can be carried out by the methods of chapters VI
and VII. Thus the elaborate and extensive calculations involved
in the planetary theory may be done almost entirely
without hand work.
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			APPENDIX A
	LIST OF CARD CATALOGUES AT THE BUREAU

The following list contains the card files at the Bureau
which are most likely to be of use to others. It does not
include the files which constitute the various steps of calcula-
tions that have been performed. The existence of such files
may be inferred fromfprobiems described in part III of
this book.

[TABLE:
			Argument			Described
File 	   Function     	Interval  Range	  Accuracy   Source 	in Chapter

1	sin x, cos x, tan        0°.01	 0 to 90°    7	     Peters	
2	   tan x 		  10'' 	 0 to 30°    7	   Brandenburg	    IX
3	tan x, sec x	  	  1^S	 0 to 20^m   7			    IX
4	tan x, sec x		0^{S}.1	 0 to 20^m   8			    IX
5	   x^{-3/2}		 0.001	 2,0 to 7.5  8	     Comrie	   XII
				 0,01	 7.5 to 20.0
6	    K/x 		0.0001	0.99 to 1.17 8			    IX
7 	Tables of the disturbing function				   XII
8 	Boss General Catalogue of 33342 Stars
9 	A.G. Catalogue from -10° to + 20°				    IX
10 	B.D. Catalogue							     X
11 	Yale Zone Catalogues
12	Kohlschiitter's Catalogue
13 	Harmonic series representing the rectangular and polar co-	    XI
		ordinates of the moon.
14	Rectangular coordinates of the major planets.			   XII
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			  APPENDIX B

	ADDITIONAL DESCRIPTION OF THE MACHINES
			AND PLUGBOARDS

	For those who may wish a more detailed knowledge of the
machines we give here a few of the principles of their operation.
This is in no sense a complete description but is in tended
merely to remoYe some of the mystery.
	The regular tabulator service men are familiar with all the
circuits of the machines and are glad to answer questions, but
usuall) this is little help as the user of the machine is not in
a position even to ask questions. It is the purpose of these
few paragraphs to give a basis for asking questions.
	The tabulator contains about 200 electromagnets. The
passage of a current through one of these magnets causes an
armature to be pulled toward the magnet, and the armature
in turn is used to operate mechanical parts of the machine

[Diagram of the function of the solonoid/relay in a punched card machine
Text with arrows indicates:
"PRIMARY
CIRCUIT" 
"ARMATURE"
"MAGNET" 
"X" 
"NX" 
"C"

Caption:
FIG. 18. Relay ]

or to open and close the electric contacts of a second circuit.
Examples of the mechanical form are the armatures which
operate the clutches of the adding wheels in the counters
and the ones which operate the type bars in the printing
mechanism. The electrical form is called a relay and is
shown schematically in figure 18.
	When there is no current flowing through the magnet
(primary circuit) a spring holds the armature so that the
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"C" and "NX" are connected, but when the armature 
is drawn toward the magnet the "C-NX" circuit is
opened and "C-X" one closed. Sometimes a single
magnet is used to operate several sets of contacts multi-
contact relay) .

[Diagram of Tabulator Plugboard and Switches

Usually 20 columns of plugboard holes inside various rectangular lines. There are 34 rows.

Text over plugboard includes
"NO. 1 BANK"	"No. 2 BANK"
Over 10 column 2 rows of plugboard each.

Remaining banks of same size numbered at side "3", "4", "5", "6", "7"

Bank of three rows (10 columns) labeled "CLASS SELECTION - 1" Rows labeled "A" "B" "C"

Bank of 1 row (7 columns) 1st column is "CC" rest is "CLASS SEL PICK UP"  pairs of holes labeled "1" "2" "3". This is next to the 1st CLASS SELECTOR 3rd row.

Another bank of three rows "CLASS SELECTION - 3" and "CLASS SELECTION - 4[2?]" otherwise as first class selector bank.

4 rows of 20 columns labeled "ADDING BRUSHES" with number over the 5th, 10th, 15th and 20th columns in each row (except for the first row for the 10th because of ADDING BRUSHES). "5", "15", "20". "25" and so on to "80".

A section of 20 columns by 4 rows but with 8 plug holes missing at bottom right (4 columns from bottom 2 rows).
First Row labeled "To ADD BR	GROUP IND. & AUTO. CTRL		TO ADD BR"
Second ROw labeled "TO COUNTER"
Third Row labeled "TO SHUNT PLUGS"
Fourth Row labeled "TO CTRL BRUSHES"

Another 4 by 20 section numbered as ADDING BRUSHES but labeled "CONTROL BRUSHES"

A somewhat irregular heavily labeled section of 3 rows of 17, 19 and 19 columns, overall label "X-DISTRIBUTOR"

1st row: 1st column "TCC", 2nd column "UCL", 3-7 columns "PLUG TO C" 8-9 unlabeled because of "X-DISTRIBUTOR" then 10 "X" 11 "NO-X" 12 "C" 13 "X" 14 "NO-X" 15 "C" and 16 and 17 MAJOR
2nd row: 1st column "M(U)" 2-3 "X-1", then 4-6 are labeled "X", "NO-X", "C" as in previous row columns 11-12 and 13-16 as are 7-9 in this row. Columns 10-11 labeled "x-4", " 12-14 and 15-17 all fall under the "X", "NO-X", "C" in the 1st row (columns 10-15). Column 18-19 labeled "INTER".
3rd row: 1st column "M(L)" 2-3 "X-2", 4-6 and 7-9 as in row 2, 10-11 "x-5", columns 12-17 fall under columns 10-15 in the 1st row and are labeled accordingly. Columns 18-19 labeled "MINOR".

Last 7 rows all labeled under the same headings, only 19 columns of holes.

1st column unlabeled. 2-9 under "COUNTER ADD OR SUBTRACT" no holes but a column of text headed "CTR
MODS" putting a number in the middle of each row except the 1st 2-7, columns 10-12 "BAL. CONT.", columns 13-14 "LIST", columns 15-19 "SYMBOL"
At bottom column 1 labeled "ADD", columns 2-5 "+", columns 6-9 "-" most of the rest is illegible but by checking later reference under column 10 "+ & -", 11 "C",  12 "+", 14 "ONLY" 15 "CTR" 16 "L" 17 "*", 19 "-" (18 may also be "*"?)

Tables down the right side include

List of two categories
"
__________________________
| TABULATE |  |  LIST |  |
|__________|__|_______|__|"

8 column Bank table

"
_______________________________
| BANK 		|1|2|3|4|5|6|7|
|_______________|_|_|_|_|_|_|_|
|LIST |LIST|OFF	| | | | | | | |
|     |____|____|_|_|_|_|_|_|_|
|     |IND | ON | | | | | | | |
|_____|____|____|_|_|_|_|_|_|_|
|        |  OFF | | | | | | | |
|        |______|_|_|_|_|_|_|_|
|        | PROG | | | | | | | |
|        |______|_|_|_|_|_|_|_|
| TOTAL  |  MIN | | | | | | | |
|        |______|_|_|_|_|_|_|_|
|        |INTER.| | | | | | | |
|        |______|_|_|_|_|_|_|_|
|        | MAJ. | | | | | | | |
|________|______|_|_|_|_|_|_|_|
|          |OFF | | | | | | | |
|  BALANCE |____|_|_|_|_|_|_|_|
|          | ON | | | | | | | |
|__________|____|_|_|_|_|_|_|_|"

Tables with three columns:
"
____________________________________
|MISC SWITCHES	|	OFF |	ON |
|_______________|___________|______|
|AUTO. START	|	    |	   |
|_______________|___________|______|	
|AUTO. RESET	|	    |	   |
|_______________|___________|______|
|AUTO. CONTROL	|	    |	   |
|_______________|___________|______|
|	-INTER.	|	    |	   |
|_______________|___________|______|
|	-MINOR	|	    |	   |
|_______________|___________|______|
|GROUP INDICATOR|	    |	   |
|_______________|___________|______|
|SUMMART PUNCH	|	    |	   |
|_______________|___________|______|
|S.P.AUTO.RESET	|	    |	   |
|_______________|___________|______|
|S. P. MAJOR	|	    |	   |
|_______________|___________|______|
|S. P. INTER	|	    |	   |
|_______________|___________|______|
|S. F. MINOR	|	    |	   |
|_______________|___________|______|
|SUCCESSIVERESET|	    |	   |
|_______________|___________|______|"

"FIG, 19. Diagram of Tabulator Plugboard and Switches
 
. . 
Figure 19 is a d1agram of the plugboar d us cd for all 
numerical tabulators. Not all the holes indicated are used
on all types of tabulators. For instance the standard type 
3S has only five counters and five print banks whereas the
diagram shows seven of each. Again, the tabulator of the
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Bureau has only four counters and five banks. In these cases
we simply remember that some of the indica ted holes have no
meaning.
	We proceed now to describe the plugboard in terms of
magnets and relays in order to make clear some features not
easily grasped otherwise. In the tabulator most of the magnets
and relays have one end of the primary circuit permanently
grounded to the machine so that only the other end
need be wired on the plugboard. In many cases two or more
holes in the plugboard are wired in "common" to the same
circuit. This eliminates the need of double and triple plugwires.
	At the top of the diagram are seven enclosures with ten
pairs of holes in each. The two holes of a pair are common
and are connected to the clutch magnets of the adding wheels
of the counters and to the printing magnets. By wiring into
these holes from the adding brushes it is possible to list and
add at the same time.
	Directly under "NO. 7 BANK" is a hole marked "CC." This
hole, wired to a position of the counter. gives the card count.
	A class selector is a multicontact relay with ten sets of
three points. The diagram shows three class selectors. The
primary circuits are each wired to a pair of common holes
labeled "CLASS SEL. PICK-UP" and the thirty secondary contacts
are wired to the groups of 30 holes labeled "CLASS
SELECTION." In using a class selector, say number one, the
adding brushes from one field are connected to the "NO-X"
holes and those from the second are wired to the "X" holes.
The "c" holes are then connected to the counter. The impulse
to operate the primary comes from the "CONTROL
BRUSHES."
	The eighty holes labeled "ADDING BRUSHES" correspond
to the eighty positions of the card.
	The next four rows of holes (" GROUP IND. & AUTO. CTRL.")
are used for the control and group indication. For purposes
of group indication there are multicontact relays with only
two contacts for each circuit. When these points are closed
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each hole in the first row is connected with the one just below.
Under these circumstances we may wire from the adding
 brushes to the top row (" TO ADD BR.") and from the second
row ("TO COUNTER") to the counter, and the result will be
the same as if the brushes were wired directly to the counter.
During the process of group indication these relay points are
normally open, but they close during the passage of the first card
of each control group. In order to be able to list individual
cards through these circuits a switch is provided which oper-
ates the relays and keeps the points closed. The sixteen
pairs of holes tot  serve the same function as the first
four on the left, and in addition are associated with the control
operation.
	For the control operation there are sixteen independent
single circuit relays. of the primary of each is connected
to one of the holes marked "TO ADD BR." and the
other end to the corresponding hole in the last row marked
"TO CTRL. BRUSHES." It is thus possible to control on any
 lS cl'
number of columns on the card not exceeding 16. The next
to the last row of holes ("TO SHUNT PLUGS") are connected to
the contact points of the control relays. The first four on the
right of the plugboard are shown in figure 20.

[Diagram illustrating shunt to plugs

arrow pointing to plug hole  
"HOLES IN PLUGBOARD
("TO SHUNT PLUGS ")"

Arrow pointing to ground symbol
"GROUNDED TO
MACHINE"

Caption:
"FIG. 20. 	Control Relay Contacts"

	The first contact point is wired permanently to the ma-
chine and the others are connected in series as shown. Thus
when the·holes in the columns of the card which are wired for
the control match, the points are all closed and we get a
circuit from the machine through all the points and out the
last hole. From here a plugwire conducts the impulse to
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another hole in the plugboard which actuates the control
mechanism. The failure of any one of the control relays to
operate thus causes a "break in control." By plugging into
three of these positions we may secure the necessary impulses
for major, intermediate, and minor control.
	These control relay points, if we skip the first, may be
used for other purposes than control. When we require
several relays in series which are operated by the comparison
of holes in two cards, we simply wire in one of the holes and
out another. (See part B of figure 21.)
	The next four rows on the plugboard are connected to the
control brushes.
	Below these, on the right, are three pairs of common holes
marked "MAJOR," "INTER," and "MINOR." These are to be
connected by plugwires to the holes above ("TO SHUNT PLUGS")
as described for the operation of the control.
	The holes "TCC," "M(U)," "M(L) " are for special circuits
used in listing cards with consecutive numbers to show that
one of the sequence is missing. The hole "UCL" is for another
special circuit. Our tabulator has neither of these features.
The other holes in the section labeled "X-DISTRIBUTOR" are
all associated with the operation of the counters. There are
five double circuit relays foF this purpose. These relays are
operated by "x" punches on the cards. For each relay there
are eight holes in the plugboard, the first two of which are
common and go to the magnet of the relay. If one of these
holes is wired to the control brush of a given column an "x"
in that column of the card will cause the relay to operate.
The next three holes are wired to the three contacts of the
first secondary circuit and the succeeding three to the other
secondary circuit. When there is no "x" on the card the
two circuits "C-NO-X" are closed, but when there is an "x"
the circuits "C-X" are closed. The five holes "PLUG TO c"
supply electric impulses which arc usually wired through the
secondary circuits of the relays and used to operate the
counters.
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	The last seven rows are used for the operation of the
counters, and list banks. The first column "ADD " is merely
a source of impulses similar to those above "PLUG TO C."
The remaining eighteen columns are associated in rows with
particular counters and list banks. Thus all the holes in the
top row refer to counter (l) and list bank (1).
	In any row the first four holes are common, so an 1mpulse
in any one of them will cause the counter to add as a given
card passes the adding brushes. Similarly an impulse in one
of the next four holes will cause it to subtract. This impulse
may come directly from "ADD," or from "C" through one of
the "x"relays.
	The three columns "BAL. CONT." control the pnnttng of
the totals in the counters. A wire from "C" to "+ & -"
will cause all totals to be printed, while a wire from "C" to
"+" will cause negative totals to be omitted from the record.
	Similarly, in listing a wire connecting the two holes
"LIST" will cause the machines to list only the cards whose
numbers would be added in the counter.
	The last five columns are used for the printing of the
asterisk and minus s1gns. In each printing unit there are
eleven type bars. The first from the left has ten digits and
an astensk, the next nine have each ten digits, and the last
one has on Yan ast ·sk and a minus sign. The asterisks are
used to indicate the totals on the record, and the minus sign
to indicate negative values of the listed items or totals. If it
is desired to print an asterisk with type bar on the left,
the numbers in that column  must be sacrificed. The choice
is made by connecting "L" with "CTR" or with "*." Simi-
larly with the bar on the right it is necessary to choose the
"*" or the "-" by means of a plugwire.
	Most of the operation of the machine is controlled by the
plugboard, but in addition there are a number of ordinary
switches and one gear shift lever. There are shown in the
right-hand part of figure 19.
	The top line shows the gear shift lever. This is used to
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change from the high tabulating speed to the lower listing
speed.
	The first group of switches shown on the chart is associated
with pnrticular counters and list banks. The first seven
switches are simple two-way switches. They have two de-
scriptions in the chart "LIST-IND.," and "OFF-ON." In case
there is a counter associated with the print bank the "OFF-ON"
applies, for in this case you can list or not list. If the bank
is simply a list bank without a counter, provision is made to
use this for additional group indication. In that case the
switch causes it to list or indicate.
	The printing of totals from a given counter is controlled
by a dial switch with he settings. The counter will print
no totals; print totals on any break in control but not reset;
print and reset on minor, intermediate, or major breaks in
control.
	Negative numbers are added as complements in the
counters. If the total to be printed is negative it may be
printed as a complement, or as a true figure with a minus sign.
The "BA.LANCE" switches are provided for this purpose.
	The "MISC. SWITCHES" control the operation of the machine
as a whole. The "AUTO. START" causes the machine to
start feeding cards after each reset. The "AUTO. RESET"
gives a reset and total print operation on each break in control.
The "AUTO. CONTROL" switches must be turned off when the
control is not wired or there will be a break in control after
each card.
	The "GROUP INDICATOR" switch operates the multicontact
relay described in the explanation of group indication. This
switch should be labeled "LIST-IND." instead of "OFF-ON."
	The first summary punch switch, if "ON," prevents the
tabulator from operating unless the summary punch is operating
properly. For instance the tabulator will not operate
if the summary punch runs out of cards.
	The next three switches determine whether summary
cards are to be punched on major, intermediate, or minor
breaks in control.
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[Partial images of plugboards with connections.

There are 4 partial plugboards shown labeled "A", "B", "C", "D"

Plugboard A

Connection from "CONTROL BRUSHES" Row 1, Column 4 to "X-1" column 2 in "X-DISTRIBUTER"
Connection from "CONTROL BRUSHES" Row 2, Column 3 to "X-2" Column 1
Connection from "PLUG TO C" column 3 to left "C" row 1 directly below
Connection from left "X" row 1 to left "C" row 2 directly across (again below "PLUG TO C")
Connection from left "NO-X" row 1 to right "C" row 2.
Connection from left "X" row 2 to row 2, column 3 of "COUNTER ADD OR SUBTRACT" directly below
Connection from left "NO-X" row 2 to row 2 column 5 "COUNTER ADD OR SUBTRACT"
Connection from right "X" row 2 to row 2 column 6 "COUNTER ADD OR SUBTRACT"
Connection from right "NO-X" to row 2 column 4 "COUNTER ADD OR SUBTRACT"

Plugboard B

Connection from "ADDING BRUSHES" row 4 column 14 (74) to GROUP IND & AUTO CTRL row 1 column 14
Connection from "ADDING BRUSHES" row 4 column 15 (75) to GROUP IND & AUTO CTRL row 1 column 15
Connection from "ADDING BRUSHES" row 4 column 16 (74) to GROUP IND & AUTO CTRL row 1 column 16

Connection from "GROUP IND & AUTO CTRL" row 3 column 14 ("TO SHUNT PLUGS" column 10) to "X-1" 	column 1.
Connection from "GROUP IND & AUTO CTRL" row 3 column 17 ("TO SHUNT PLUGS" column 13) to "PLUG 	TO C" Column 4 (X-DISTRIBUTOR COLUMN 4)
Connection from "GROUP IND & AUTO CTRL" row 4 column 14 ("TO CTRL. BRUSHES" column 10) to 	"CONTROL BRUSHES" row 4, column 14 (74).

Connection from "GROUP IND & AUTO CTRL" row 4 column 14 ("TO CTRL. BRUSHES" column 10) to 	"CONTROL BRUSHES" row 4, column 14 (74).
Connection from "GROUP IND & AUTO CTRL" row 4 column 15 ("TO CTRL. BRUSHES" column 11) to 	"CONTROL BRUSHES" row 4, column 15 (75).
Connection from "GROUP IND & AUTO CTRL" row 4 column 16 ("TO CTRL. BRUSHES" column 12) to 	"CONTROL BRUSHES" row 4, column 16 (76).

Connection from "PLUG TO C" column 3 to left "C" directly below.
Connection from left "NO-X" below "PLUG TO C" to "COUNTER ADD OR SUBTRACT" row 1, column 3.

PLUGBOARD C

Connection from "CONTROL BRUSHES" row 4 column 13 (73) to "X-DISTRIBUTOR" row 1 column 7 ("X-3" column 2)
Connection from to the right of X-DISTRIBUTOR right "NO-X" to column 1 of MINOR.
Connection from to the right of X-DISTRIBUTOR right "C" to column 1 of "MAJOR" (directly to right).

PLUGBOARD D

Connection from "CONTROL BRUSHES" row 2 column 12 (32) to arrow to text indicating "to C.C."
Connection from "CONTROL BRUSHES" row 3 column 12 (52) to row 1, column 1 of "ADD" column to 	the left of "COUNTER ADD OR SUBTRACT".
Connection from "CONTROL BRUSHES" row 3 column 13 (53) to X-DISTRIBUTOR column 7 (X-3 column 	2)
Connection from "CONTROL BRUSHES" row 4 column 12 (72) to row 1, column 1 of "ADD" column to 	the left of "COUNTER ADD OR SUBTRACT".
Connection from to the right of "X-DISTRIBUTOR" right "X" to column 1 "MINOR"
Connection from to the right of X-DISTRIBUTOR right "C" to column 1 of "MAJOR" (directly to right).
 
Caption
"FIG. 21. Special Wiring Diagram for Tabulator"]
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	The last switch is used mostly in bookkeeping work. It
causes a total to be printed twice if the break in control is
both minor and major.
	We conclude the discussion of the tabulator with some of
the special wiring arrangements already mentioned in the
text (see fig. 21). Only the unusual parts of the wiring are
shown, and it is assumed for instance that the adding brushes
are wired to the counters.
	Part "A" of the diagram shows the operation where
counter (2) is to subtract algebraically numbers from the
cards of one set from those of another. The set which is to
be subtracted has an "x" in column 4 while both sets have
an "x" in column 23 for negatiYe numbers.
	Part "B" of the diagram shows part of the wiring for
adding in counter (1) only the numbers from the first card of
each control group, but without stopping to break control
(see ch. III, sec. 2). Columns 74-76 are used to determine
the control groups. The wire shown on the diagram as leading
from "TO SHUNT PLUGS" to "PLUG TO c" should lead
instead to column 17 of the" CONTROL BRUSHES." The cards
should contain an" x" in column 17.
	Wiring "C" is for the case where special control cards are
interspersed between the detail cards and we want one break
in control with each control card (see ch. III, sec. 6). The
control cards must have an "x" in column 73 and the detail
cards not. The major control switch must be turned on and
the minor off.
	Wiring" D" will cause a break in control after every other
card (see ch. V, sec. 4). For this operation the group indication
and major control switches must be on and the minor
off. The wire labeled "to C.C." is to be connected to the card
count on the plugboard.

The Reproducer

	Figure 22 shows the arrangement of the plugboard of the
reproducer. This diagram should be compared with figure 11.
As one would expect from figure 11, reproducing requires
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[Diagram of Reproducer Plugboard

"READING BRUSHES"
heads section of 4 rows by 20 columns plug holes "5", "10" and so on are numbered up to 80 (left to right from top row).

Next section left and right half of boards separated are the holes.
On left half:
"CLASS SELECTION-3"
4 rows of 10 columns of holes except for row 4 which has only 2 columns
To the right of column 1 the 1st 3 rows are labeled "NO X" "X" and "C" respectively.

On the right half:
"X-ELIMINATOR" 
3 rows of 5 columns (rest of half has no holes)
row 1 headed by "PUNCH 11 & 12"
row 2 headed by "PUNCH 0 TO 9"
row 3 headed by "COMMON"

Below X-ELIMINATOR (to the right of row 4 of "CLASS SELECTION-3") are 2 pairs of holes labeled "R-X" and "P-X" respectively.

Next row is split into one group of 5 columns and one group of 2 columns of holes.

Left group of 5 labeled "ZERO", right GROUP of 2 labeled "GP-CH. CL. SET UP" left "P-D" right "R-D".

Next section:
"PUNCHED MAGNETS"
heads section of 4 rows by 20 columns plug holes "5", "10" and so on are numbered up to 80 (left to right from top row). Except that "10" is occluded by heading.

Next row divided in half 10 columns each labeled "READING X BRUSHES" and "PUNCH X BRUSHES" in both halves the holes "1", "2", "3", "8", "9" and "10" are labeled.

Next section:
"PUNCH BRUSHES"
heads section of 4 rows by 20 columns plug holes "5", "10" and so on are numbered up to 80 (left to right from top row). Except that "10" is occluded by heading.

Next section is divided into two halves one is "CLASS SELECTION-1"  the other "CLASS SELECTION 2" these are as "CLASS SELECTION-3" except that "CLASS SELECTION-2" has the labels of its rows to the left of the right most colunm and the two holes in row 4 are on the right for it (and on the left for the other 2).

Next section:
Has a heading on the left "COMPARING RELAYS" and on the right "TO PUNCH BRUSHES".
heads section of 4 rows by 20 columns plug holes "5", "10" and so on are numbered up to 80 (left to right from top row). Except that "5" and "15" is occluded by the two heading.

Next section:
Has a heading on the left "COMPARING RELAYS" and on the right "TO COMPARING BRUSHES".
heads section of 4 rows by 20 columns plug holes "5", "10" and so on are numbered up to 80 (left to right from top row). Except that "5" and "15" is occluded by the two heading.

Next section:
"COMPARING BRUSHES"
heads section of 4 rows by 20 columns plug holes "5", "10" and so on are numbered up to 80 (left to right from top row). Except that "10" is occluded by heading.

Caption for diagram reads:
"FIG. 21.	 Diagram of Reproducer Plugboard"]
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wires from "READING BRUSHES" to "PUNCH MAGM£TS," while
gang punching requires wires from "PUNCH BRUSHES" to
"PUNCH MAGNETS."
	There are eighty "COMPARING RELAYS" each with two
holes in the plugboard. In checking the work of the reproducer
one of these holes must be connected to a comparing
brush and the other to the corresponding punch brush. In
checking the work of the gang punch one is connected to the
reading brush and the other to a comparing brush.
	The class selectors are the same as on the tabulator. The
two holes at the bottom are common and are for the wire to
the "x" brush.
	There may be from 6 to 10 "x" brushes which precede the
reading brushes, and from 6 to 10 more which precede the
punch brushes. The holes in the plugboard for them are just
above thm:e for the punch brushes. They are used to operate
the class selectors, and to prevent the punching of a card with
an "x," as described below.
	The "X-ELIMINATOR" is used to separate "x" and "y"
in a given column from the digits 0 to 9. For instance if we
have x and y and the digits 0 to 9 in the same column we can
form a new card with the x and y in one column and the 0 to
9 in another. This can be done for five different columns at
the same time.
	The holes marked "ZERO" furnish a source of "zero"
impulses just as the card count of the tabulator gave a "1."
	The pair of common holes "P-X" when wired to a "PUNCH
x BRUSH" will supress all punching for cards with an "x"
(see ch. IV, sec. 3). A special switch makes it possible to
punch only when there is or only when there is not an "x."
The holes "R-X" serve a similar purpose when connected to
a "READING X BRUSH."

The Multiplier

	The plugboard of a new type standard machine is shown in
figure 23. This should be studied in connection with the
description in chapter II.
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[Diagram of a Multiplier plugboard


First section:
"BRUSHES"
heads section of 4 rows by 20 columns plug holes "5", "10" and so on are numbered up to 80 (left to right from top row). Except that "10" is occluded by heading.


Next section split into left and right sections each 2 rows by 8 columns separated by a gap in the middle. The plug holes in the top row are connected to the one directly below by a line.
On the left a heading reads "MULTIPLIER".
On the right a heading reads "MULTIPLICAND".

Next section a row of 12 plug holes centred in the diagram with the heading:
"CROSS-FOOTING R H C - TO BRUSHES"

Next section a row of 12 plug holes centred in the diagram with the heading:
"CROSS-FOOTING L H C - TO 8RUSHtS"

Next section a row of 10 plug holes centred in the diagram with the heading:
"SUMMARY COUNTER-T0 BRUSHES"

Next section a row of 10 plug holes centred in the diagram with the heading:
"SUMMARY CTR TO PROD. CTR"

Next section a row of 16 plug holes centred in the diagram with the heading:
"PRODUCTS COUNTER-TO SUMMARY COUNTER"

Next section a row of 16 plug holes centred in the diagram with the heading:
"PRODUCTS COUNTER TO PUNCH"

Next section a row of 10 plug holes centred in the diagram with the heading:
SUMMARY COUNTER TO PUNCH


Next section
"PUNCH"
heads section of 4 rows by 20 columns plug holes "5", "10" and so on are numbered up to 80 (left to right from top row). Except that "10" is occluded by heading.

Next section
Seperated into left and right:

The narrow left section headed "½ PICK·UP" consists of 5 irregular rows.
1st row blank.
2nd row two plugholes one on extreme left labeled "IND" the other on the extreme right labeled "GROUP"
3rd row fills the space with 4 plug holes and the heading "X-BRUSHES" plugholes are labeled "I"  "4" in the appropriate columns, the heading occludes 2 and 3.
4th row has two plugholes on the left labeled "5", "6" respectively.
5th row has one plughole only on the extreme right labeled "C".

The broad right section headed "½ ENTRY" consists of 5 rows:
1st row 16 columns bordered above, below and to the left by a line
2nd-5th rows are all 12 columns 
Between 2nd and 3rd row a line splits to two arrows pointing at row 2, column 4 and row 3, 	column 4 with the text "NORMAL PLUGGING"
Between 3rd and 4th row a line splits to two arrows pointing at row 3, column 3 and row 4, 	column 3 with the text "A+B-C PLUGGINGS"
Between 4th and 5th row a line splits to two arrows pointing at row 4, column 2 and row 5, 	column 2 with the text "NORMAL PLUGGING"


The next section is split into six sets of columns, all of 4 rows.

On the left is a section headed "FIXED" with 1 column of holes, between row 1 and 2 the text 	"OFF", between 2 and 3 the text "MULTIPLIER" between 3 and 4 the text "ON"

2nd section headed "GROUP" with 2 columns, between rows 1 and 2 the text "OFF", between 2 and 	3 the text "MULTIPLIER" between 3 and 4 the text "ON".

3rd section headed "MULTIPLY CROSS-FOOT" with 4 columns, between rows 1 and 2 the text 	"MULTIPLY", between 3 and 4 "CROSS-FOOT".

4th section headed "UNIT" 1 column, betweeen rows 1 and 2 the text "OFF", between 2 and 3 the 	text "MULTIPLIER", between 3 and 4 the text "ON".

5th section headed "RHC TO LHC" 1 column, between rows 1 and 2 the text "ADD", between 3 and 4 	the text "SUBT."

6th section headed "DEC ACC" 3 columns, between rows 1 and 2 the text "OFF", between 2 and 3 	the text "DEVICE", between 3 and 4 the text "ON".

Next section two sections of plug holes on the left and some text on the right.

1st section on the left headed "CR0SS·FOOT" is 2 columns, between rows 1 and 2 the text "OFF", 	between rows 2 and 3 the text "SUM CTR", between rows 3 and 4 the text "ON".

2nd section on the left headed "PROD TO" is 2 columns, between row 1 and 2 the text "ADD", 	between rows 2 and 3 the text "SUM CTR", between rows 3 adn 4 the text "SUBT."

Text on the right:
"	PLUGBOARD DIAGRAM
	      ______

    AUTOMATIC MULTIPLYING PUNCH

	TYPE 600 - TYPE 601"


Caption of Diagram:
"FIG. 23. Diagram of Plugboard for Multiplying Punch"  ]
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	The eighty "BRUSHES"are used to read numbers from the
card as it is fed into the machine. These numbers may be
read into (a) counter (1) "MULTIPLIER," (b) counter (2)
"MULTIPLICAND," (c) counter (4) and four positions of counter
(3) "R. H. C.," (d) counter (6) and four positions of counter
(5) ''L. H. C.," and (e) counter (7) "SUMMARY COUNTER."
	The computed products are read into counter (7) through
the tenth row of holes" SUMMARY CTR." from counters (4)-(5)
here labeled "PRODUCT COUNTER." They are also punched
on the card from the same counters through the twelfth row
of holes. It is also possible to punch from counter (7) through
the thirteenth row of holes. The wires from rows 12 and 13
go, of course, to the eighty holes labeled "PUNCH."
	In order to "round" a computed quantity a "5" may be
added in any one of the sixteen positions "½ ENTRY," either
for individual cards" IND." or for each "GROUP." The group
would be determined by an "x" punch on a card. For con-
venience there are 6 "X-BRUSHES" and we can plug from any
one of these to "C."
	The four rows of twelve holes labeled "NORMAL PLUGGING"
and "A + B - C" are associated with multiplication and
crossf ooting. For normal multiplication these holes are
wired to connect the multiplying mechanism with the right-hand
and left-hand components counters. Since this wiring is
through the plugboard it is possible to change it for unusual
operations such as forming two products at one time.
	The remainder of the holes in the plugboard are used instead
of switches for controlling the machine. For instance
if we want a fixed multiplier we wire together the two bottom
holes in that box. If we do not want a fixed multiplier we
must connect the top pair. Some of these switches require
two circuits or even four.
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		APPENDIX C

	THE ORGANIZATION OF THE THOMAS J.
	   WATSON ASTRONOMICAL COM-
	      PUTING BUREAU

	Thc Bureau is a Joint enterprise of the American Astro- 
nomical Society, The International Business Machines Cor-
poration, and the Department of Astronomy of Columbia
University. Its operation is in charge of a Board of Managers,
the majority of whom are appointed by the American Astro-
nomical Society.*
	The operation of the Bureau is governed by the following
regulations: 
	1.  The Bureau shall be under the direction of a Board of
Managers which shall have charge of its operation. There
shall also be an Advisory Council of five which shall have
duties set forth in paragraph 8 .
	2. The purpose of the Bureau is to perform astronomical
calculations by means of the machines of the International
Business Machines Corporat1on.
	3. The machines and all persons employed by the Bureau
shall be under the immediate charge of the Board of Managers
or of such person as the Board may appoint.
	4. The Board of M~nagers shall make all decisions as to
 the acceptance or rejection of applications to perform compu-
tations. It may ca on one or more members of the Advisory
Council for advice on any application and shall do so in
the case of an application from  a member of the Board of
Managers.
	5. The Bureau shall make a charge for all computations
performed by its machines or operatives. The Board of

	* See also Pub. Ast. Soc of the Parific, vol. XLIX No. 291, Oct. 1937.
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Managers shall furnish a maximum estimate of cost to an
applicant whose project is deemed acceptable.
	6. T he Board of fivc Managers shall be appointed as
follows: one by the International Business Machines Corporation,
one by the Department of Astronomy of Columbia
University, and three by the Advisory Council. The three
appointed by the Advisory Council shall each serve a term of
three years and shall not be immediately eligible for reappointment
after serving two terms or portions of terms. The
terms shall be so arranged that one of the three is appointed
every year. A vacancy occurring before the expiration of any
term shall be filled by the Advisory Council for the unexpired
portion of that term. A member of the Advisory Council
shall be eligible for membership on this Board, but if so
appointed shall vacate his seat on the Council.
	The Board of Managers shall appoint a Chairman, Secretary,
and Treasurer from its own members. The presence of
three members shall constitute a quorum for the transaction
of business. Majority ruling of those present shall prevail,
the Chairman having an extra vote in the case of a tie.
	7. The President of the American Astronomical Society
shall be requested to appoint the five members of the Advisory
Council, each to serye for a term of five years and not to be
eligible for immediate reappointment . The terms shall be so
arranged that one member shall retire each year. A vacancy
before the end of any term is to be filled by the Advisory
Council for its unexpired portion. The Council shall elect its
own Chairman.
	8. The duties of the Advisory Council shall be: (1) to appoint
three Managers as set forth in paragraph 6; (2) to give
advice when requested by the Chairman of the Board of
Managers; (3) to meet at least once a year with the Board
of Managers at which meeting the Chairman of the latter shall
report on the work of the Bureau; (4) to act as trustees of any
funds or property acquired by the Bureau, decision as to the
use of such property, however, to rest with the Board of
Managers.
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9. Changes in these regulations may be made by corre-
5μondencc Y .1 b a unanimous vote of the Board of Managers and
the Advisory Council or by a favorable vote of six at a com-
bined meeting of the Board and Council provided  notice of
such change has been sent to every member of both at least
one month in advance of the meeting.
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			INDEX

Accounting machines, 7, 12		Cards, master, 9	
Adding brushes 12, 118				used for recording observations, 93
	wheel, 12, 116				used for perturbations, 104
Alphabetic machine,, 6, 22			punched for star catalogues, 82, 92
Analysis, harmonic, 57, 62			punching positions on, 6
Antilogarithms, 114				rate, 19
Armatures, 116					ruling on, 7, 34
Astronomical applications, 2, 79		summary, 17
Automatic control. 13, 118			tabulating, 4, 6
	plugboard, 13, 20			"x," 16
	spacing on key punch, 9		Carry-over mechanisms, 23	
Averages, 2				Catalogues, card, 115	
						star, construction, 79, 92
Background correction, 93		Check sum;, 26	
Balance control, 121			Checking with reproducer, 2I	
Banks, printing, 16				with mechanical verifier, 9
Bessel's formula, 41, 46			tables with differences, 38
Break control, 15, 34, 49, 120		Checks, 27, 28	
Brushes, adding, 12, 118		Class selector, 16, 118	
	comparing, 21			Clerical work, 34	
	control, 13, 118		Clutches, 116	
	movable, 10, 17			Collating device on sorter, 24	
	reading, 4, 21, 126		Collator, 23, 41	
	sorting, 10			Common multiplier, 19	
	"x", 18, 21			Comparing brushes, 21	
	wire, 4					relays, 21, 126
					Components counters, 18	
Calculation control switch,8, 22, 77,	Computing methods, 1	
	108				Constant multiplier, 19, 29	
Card count, 16, 118			Contact roller, 4	
	file at Bureau, 115		Control, automatic, 13, 118	
	gauge, 26				balance, 121
	punch, 8				break, 15, 34, 120
	punching speed, 9			intermediate, 24
Card for breaking control, 34, 49		major, 15, 120
	colored, 34				minor, 15, 120
	control, 34				pass, 15
	detail, 21, 23				"x," 16
	fields on, 7			Control brushes, 13, 118	
	used in intergration. 104		card, 34, 49
	interpretation on, 7			cards for interpolation, 49
						group, 15, 16
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Control relays, 119				Gearshift lever, 121	
	switch, calculation, 8. 22, 77, 108	General perturbations, 101, 112	
Conversion of coordinates, 81			Group indication. 16, 118, 122	
Correlation coefficients, 2, 31, 33					
Cos (ijh) table, 59				Half pick up, 19	
Cosine series, 58, 60, 6S			Harmonic analysis, 57, 62	
Counters, 12, 18					series, 2, 57, 65, 102
Counting attachment, 11				High-speed listing dedce, 24	
	sorter, 10					reproducer, 20
Critical tables, 36					summary punch, 24
Crossfooting, l7, 19, 20, 30			Hollerith method. 1, 4	
					
Dealing cards, 24				Integration, numerical, 22, 43, 52, 75,	
Detail cards, 21, 23					105
Difference iormulae, 43				Interpolation, 41, 43	
Differences, 2, 39, 50, 55				bibliography, 43
Differential equations, 75, 97, 1O1			formulae, 41, 46
Differentiation, numerical, 43				by Lagrangian method, 47
Digit emitter, 24				Interpreter, 7, 21	
Digiting, progTessive, 31					
Disturbing function, 113			Key punch, duplicating, 8, 16	
Double punching, 10				Keyboard computing machines, 25	
Duplicates, 21					
Duplicating key punch, 8			Lagrangian formulae, 43, 46	
	rack, 9						method of interpolation, 47
Duplicator, 8					Least squares, 31, 33, 94	
						List bank, 16	
Eject key on card punch. 9			Listing, 16	
Electrical multiplication table, 18			check, 27
Electromagnets, 4, 12, 116			Logarithms, 1, 37	
End figure method, 54				Lunary theory, 97	
Equations, differential, 75, 97, 101					
	normal, 2, 95				Machine methods, 25	
Exponential series, 65				Machines, accounting, 7, 12	
							not used at Bureau, 22
Feeding mechanism on tabulator, 12,		Magnitude of stars, determination, 90	
		15				Major control, 15	
	on multiplier, 17			Master card, 9, 17, 21, 108	
	operation on reproducer, 21		Mechanical quadrature, 43	
Fields on cards, 7					verifier, 9, 26
Filing, automatic, 23				Minor control, 15	
Fixed multiplier, 19				Minus signs, 121	
Fourier series, 65				Movable brushes, 10, 17	
Fractional wheels, 24				Multiplicand counters, 18	
Functions, tabular, 2, 36			Multiplication, check, 28	
							short cuts, 29
Gang punching, 21				Multiplication of series, 65, 74	
Gauge, card, 26					Multiplier, 17	
							counters, 18
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Multiplier, operation, 18		Punching knives, 20
	plugboard, 126				positions on cards, 6
	special features, 19, 11
	speed, 20			Quadrature, mechanical, 43
Multiplying punch, 17, 29, 56
					Rate card, 19
Negative numbers, 6			Reading brushes on tabulator, 4
Normal equations, 2, 95				on multiplier, 17
Numerical checks, 25				on reproducer, 21, 126
	differentiation, 43		Rectangular coordinates for star cata-
	harmonic analysis, 57			logues, 86
	integration, 22, 43, 52, 75, 105Reject on sorter, 11
					Rejection limit for products, 67, 72
Observations, recordings, 3, 93		Relay, 116
Off punching, 26				contacts, control, 119
						multicontact, 117
Partial products, 18			Relays, comparing, 21
Pass control, 15			Reproducer, check, 28
Periods, sorting for, 31, 49			high-speed, 20
Perturbations, planetary, 101, 112		speed, 21
Photometry, stellar, 90			Reproducer plugboard, 20, 124
Planetary motion, 2, 44, 78		Reset time, 16
	perturbations, 101		Rolling total machine, 56
Plate constatns, 94			Rounding, 19
Plugboard, automatic, 13, 17		Ruling on cards, 7, 34
	multiplier, 17, 126
	reproducer, 20, 124		Secular terms in precession, 82
	tabulator, 12, 117		Selector, class, 16, 118
Plugwires, 12				Series, cosine, 58, 60, 65
Power series, multiplication, 65		exponential, 65
Precession, computation, 82			Fourier, 65
Printing banks, 16				harmonic, 2, 57, 65, 102
	mechanism, 12, 16			multiplication, 65
Progressive digiting, 31			power, 65
	totals, 32, 52			Shunt plugs, 119
Proofreading to check punching, 26,	Simpson's rule, 47
	27				Sin (ijh) table, 59
Punch, card, 8				Sine series, 65
	key duplicating, 8, 16		Sorter, 10, 23
	multiplying, 17, 29, 56		Sorting brush, 10
	summary, 16, 24, 76			check, 27
Punch brushes, 126				needle, 27
	magnets, 126				technique, 31
Punched card method, 1, 25		Spacer key on card punch, 9
Punching, double, 10			Special perturbations, 102
	gang, 21, 30			Speed, interpreter, 21
	off, 26					listing, 16
Punching speed, 9				multiplier, 20
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Speed, punching, 9			Tabular functions, tables of, 36	
	reproducer, 21			Tabulating card, 4, 6	
	reset, 16				operation, 16
	sorting, 12			Tabulation checks, 27	
	summary punch, 17		Tabulator, alphabetic, 6, 22	
	tabulator, 16				Hollerith. 1, 12
Star catalogue, construction, 79, 92	Test scoring machine, 23	
Stellar photometry, 90			Type bars, 116	
Stirling's formula, 41, 46					
Subtraction tabulator, 12, 13					
Summary card, 17, 28			Verifier, mechanical, 9, 26	
	products counters, 18				
	punch, 16				
	punch, high-speed, 24		Watson Computing Burcau, iii, 7, 129	
Summation with tabulator, 36		Wire brushe, 4	
Superfluous difference method of in-	Wiring, 117, 124	
	terpolation, 54				
Switch, control, calculation, 8, 22, 77,"x" brushes, 18, 21	
	108					cards, 16
Switches, miscellaneous, on tabulator,		control, 16
	122					distributor, 120
Symbol printing, 121				eliminator, 126
Synthesis, harmonic, 57, 64			punch, 6, 19
					
Table, critical, 36					
	subdivision of, 53		"y" punch, 6	
	use in punched card form, 39	Yale zone catalogue of stars, 79
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